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We present a review of experimental and theoretical studies of the anomalous Hall eflfect ( AHE) , 
focusing on recent developments that have provided a more complete framework for understanding 
this subtle phenomenon and have, in many instances, replaced controversy by clarity. Synergy 
between experimental and theoretical work, both playing a crucial role, has been at the heart 
of these advances. On the theoretical front, the adoption of Berry-phase concepts has estab- 
lished a link between the AHE and the topological nature of the Hall currents which originate 
from spin-orbit coupling. On the experimental front, new experimental studies of the AHE in 
transition metals, transition-metal oxides, spinels, pyrochlores, and metallic dilute magnetic semi- 
conductors, have more clearly established systematic trends. These two developments in concert 
with first-principles electronic structure calculations, strongly favor the dominance of an intrinsic 
Berry-phase-related AHE mechanism in metallic ferromagnets with moderate conductivity. The 
intrinsic AHE can be expressed in terms of Berry-phase curvatures and it is therefore an intrin- 
sic quantum mechanical property of a perfect cyrstal. An extrinsic mechanism, skew scattering 
from disorder, tends to dominate the AHE in highly conductive ferromagnets. We review the 
full modern semiclassical treatment of the AHE which incorporates an anomalous contribution 
to wavepacket group velocity due to momentum-space Berry curvatures and correctly combines 
the roles of intrinsic and extrinsic (skew scattering and side-jump) scattering-related mechanisms. 
In addition, we review more rigorous quantum-mechanical treatments based on the Kubo and 
Keldysh formalisms, taking into account multiband effects, and demonstrate the equivalence of all 
three linear response theories in the metallic regime. Building on results from recent experiment 
and theory, we propose a tentative global view of the AHE which summarizes the roles played by 
intrinsic and extrinsic contributions in the disorder-strength vs. temperature plane. Finally we 
discuss outstanding issues and avenues for future investigation. 
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naturally in the first microscopic theory of the AHE by 



Karplus and Luttinger (Karplus and Luttinger 19541. 
However, because understanding of these concepts, not 
to mention the odd intrinsic dissipationless Hall current 
they seemed to imply, would not be achieved for another 
40 years, the AHE problem was quickly mired in a contro- 
versy of unusual endurance. Moreover, the AHE seems 
to be a rare example of a pure, charge-transport problem 
whose elucidation has not - to date - benefited from the 
application of complementary spectroscopic and thermo- 
dynamic probes. 



anomalous Hall mechanisms 






B. Kubo formalism 


EH 




1. Kubo technique for the AHE 


im 


1 


2. Relation between the Kubo and the semiclassical 




V 

O 


formalisms 






C. Keldysh formalism 


|43] 




D. Two-dimensional ferromagnetic Rashba model - a 




it 


minimal model 


ES] 


■o 


Conclusions: Future problems and perspectives on 




« o 

U. 5 
O — : 


AHE 


Ez] 


u > 

LU 


Acknowledgements 


Eo] 


"5 
z 


References 


Eq] 






5 10 15 

Magnetic Induction B kilogauss 



I. INTRODUCTION 

A. A brief history of the AHE and new perspectives 

The anomalous Hall effect has deep roots in the his- 
tory of electricity and magnetism. In 1879 Edwin H. 
Hall (Hall 18791 made the momentous discovery that. 



when a current-carrying conductor is placed in a mag- 
netic field, the Lorentz force "presses" its electrons 
against one side of the conductor. One year later, he 
reported that his "pressing electricity" effect was ten 
times larger in ferromagnetic iron (Hall 1881 1 than in 
non-magnetic conductors. Both discoveries were remark- 
able, given how little was known at the time about how 
charge moves through conductors. The first discovery 
provided a simple, elegant tool to measure carrier con- 
centration more accurately in non-magnetic conductors, 
and played a midwife's role in easing the birth of semi- 
conductor physics and solid-state electronics in the late 
1940's. For this role, the Hall effect was frequently called 
the queen of solid-state transport experiments. 

The stronger effect that Hall discovered in ferromag- 
netic conductors came to be known as the anomalous Hall 
effect (AHE). The AHE has been an enigmatic problem 
that has resisted theoretical and experimental assault for 
almost a century. The main reason seems to be that, 
at its core, the AHE problem involves concepts based on 
topology and geometry that have been formulated only in 
recent times. The early investigators grappled with no- 
tions that would not become clear and well defined until 



much later, such as the concept of Berry-phase (Berry 
1984). What is now viewed as Berry phase curvature. 



FIG. 1 The Hall effect in Ni [ data from A. W. Smith, Ph ys. 
Rev. 30, 1 (1910)]. [From Ref. [Pugh and Rostoker[ [l953l] 



Very early on, experimental investigators learned that 
the dependence of the Hall resistivity p^y on applied per- 
pendicular field Hz is qualitatively different in ferromag- 
netic and non-magnetic conductors. In the latter, p^y 
increases linearly with Hz, as expected from the Lorentz 
force. In ferromagnets, however, p^y initially increases 
steeply in weak Hz, but saturates at a large value that 
is nearly iJ^-independent (Fig. [ij. Kundt noted that, 
in Fe, Co, and Ni, the saturation va lue is roughly pro- 
portional to the magnetization Mz ( [Kundt 18931 and 
has a weak anisotropy when the field {z) direction is ro- 
tated with respect to the cyrstal, corresponding to the 
weak magnetic anisotropy of Fe, Co, and Ni (|Web- 



ster 



1925). Shortly thereafter, experiments by Pugh and 



coworkers (Pugh 1930 Pugh and Lippert 1932) estab- 
lished that an empirical relation between p^y, Hz, and 



Pxy 



R„Hz + RsMz 



(1.1) 



applies to many materials over a broad range of external 
magnetic fields. The second term represents the Hall ef- 
fect contribution due to the spontaneous magnetization. 
This AHE is the subject of this review. Unlike i?o, which 
was already understood to depend mainly on the density 
of carriers, Rs was found to depend subtly on a variety 
of material specific parameters and, in particular, on the 
longitudinal resistivity pxx — P- 



later dubbed "anomalous velocity" by Luttinger, arose 



In 1954, Kar plus and Luttinger (KL) (Karplus and 
Luttinger 1954) proposed a theory for the AHE that, in 



3 



hindsight, provided a crucial step in unraveling the AHE 
problem. KL showed that when an external electric field 
is applied to a solid, electrons acquire an additional con- 
tribution to their group velocity. KL's anomalous velocity 
was perpendicular to the electric field and therefore could 
contribute to Hall effects. In the case of ferromagnetic 
conductors, the sum of the anomalous velocity over all 
occupied band states can be non-zero, implying a contri- 
bution to the Hall conductivity a^y. Because this contri- 
bution depends only on the band structure and is largely 
independent of scattering, it has recently been referred to 
as the the intrinsic contribution to the AHE. When the 
conductivity tensor is inverted, the intrinsic AHE yields 
a contribution to pxy ~ axy/cr^x ™d therefore it is pro- 
portional to p^. The anomalous velocity is dependent 
only on the perfect crystal Hamiltonian and can be re- 
lated to changes in the phase of Bloch state wavepackets 
when an electric field causes them to evolve in crystal 



momentum space (Bohm et al. 2003 


Chang and Niu 


1996 


Sundaram and Niu 1999 Xiao and Niu 


2009 


). As 



mentioned, the KL theory anticipated by several decades 
the modern interest in Berry phase and Berry curvature 
effects, particularly in momentum-space. 
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FIG. 2 Extraordinary Hall constant as a function of re- 
sistivity. The shown fit has the relation Rs ~ p^"- [From 
Ref. [KooI] p54l] 



Early experiments to measure the relationship between 
Pxy and p generally assumed to be of the power law form, 
i.e., Pxy ~ p^, mostly involved plotting pxy (or Rg) vs. 
p, measured in a single sample over a broad interval of T 



(typically 77 to 300 K). As we explain below, competing 
theories in metals suggested either that /3 = 1 or /? = 2. 



A compiled set of results was published by Kooi (Kooi 
19541 (Fig. [2]). The subsequent consensus was that such 
plots do not settle the debate. At finite T, the carriers 
are strongly scattered by phonons and spin waves. These 
inelastic processes - difficult to treat microscopically even 
today - lie far outside the purview of the early theories. 
Smit suggested that, in the skew-scattering theory (see 
below), phonon scattering increases the value P from 1 
to values approaching 2. This was also found by other 
investigators. A lengthy calculation by Lyo ( Lyo[ 19731 
showed that skew-scattering e^tT^Qu (the Debye tem- 
perature) leads to the relationship pxy ~ (p^ + ap) , with 
a a constant. In an early theory by Kondo considering 
skew scattering from spin excitations ( Kondo| 
may be seen that pxy 



also varies as p^ at finite T, 



1962), it 



The proper test of the scaling relation in comparison 
with present theories involves measuring pxy and p in a 
set of samples at 4 K or lower (where impurity scattering 
dominates). By adjusting the impurity concentration rii, 
one may hope to change both quantities sufficiently to 
determine accurately the exponent /3 and use this iden- 
tification to tease out the underlying physics. 

The main criticism of the KL theory centered on the 
complete absence of scattering from disorder in the de- 
rived Hall response contribution. The semi-classical AHE 
theories by Smit and Berger focused instead on the in- 
fiuence of disorder scattering in imperfect crystals. Smit 
argued that the main source of the AHE currents was 



(|Smit 


1955 


1958 


Pxx 


p = l}. Berj 



AHE picture predicted that Rs 
the other hand, argued that the main source of the AHE 
current was the side-jump experienced by quasiparticles 
upon scattering from spin-orbit coupled impurities. The 
side-jump mechanism could (confusingly) be viewed as a 
consequence of a KL anomalous velocity mechanism act- 
ing while a quasiparticle was under the influence of the 
electric field due to an impurity. The side-jump AHE cur- 
rent was viewed as the product of the side-jump pe r scat- 

One 



1970). 



tering event and the scattering rate ( Berger 
puzzling aspect of this semiclassical theory was that all 
dependence on the impurity density and strength seem- 
ingly dropped out. As a result, it predicted Rg ^ Pxx 
with an exponent /3 identical to that of the KL mech- 
anism. The side-jump mechanism therefore yielded a 
contribution to the Hall conductivity which was seem- 
ingly independent of the density or strength of scatter- 
ers. In the decade 1970-80, a lively AHE debate was 
waged largely between the proponents of these two ex- 
trinsic theories. The three main mechanisms considered 
in this early history are shown schematically in Fig. [3] 

Some of the confusion in experimental studies stemmed 
from a hazy distinction between the KL mechanism and 
the side-jump mechanism, a poor understanding of how 
the effects competed at a microscopic level, and a lack of 
systematic experimental studies in a diverse set of mate- 
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a) Intrinsic deflection 

Interband coherence induced by an 
external electric field gives rise to a 
velocity contnbution perpendicular to 
the field direction. These currents do 
not sum to zero in ferromagnets. 




Electrons have an anomalous velocity perpendicular to 
the electhc field related to their Berry's phase curvature 



b) Side jump 




The electron velocity is deflected in opposite directions by the opposite 
electhc fields experienced upon approaching and leaving an Impuhty. 
The time-integrated velocity deflection is the side jump. 



c) Skew scattering 

Asymmetric scattehng due to 
the effective spin-orbit coupling 
of the electron or the impuhty. 




FIG. 3 Illustration of the three main mechanisms that can 
give rise to an AHE. In any real material all of these mecha- 
nisms act to influence electron motion. 



rials. 

One aspect of the confusion may be illustrated by con- 
trasting the case of a high-purity mono-domain ferromag- 
net, which produces a spontaneous AHE current propor- 
tional to Mz, with the case of a material containing mag- 
netic impurities (e.g. Mn) embedded in a non-magnetic 
host such as Cu (the dilute Kondo system). In a field H, 
the latter also displays an AHE current proportional to 



on the AHE problem starting around 1998. Theoretical 
interest in the Berry phase and in its relation to trans- 
port phenomena, coupled with many developments in the 
growth of novel complex magnetic systems with strong 
spin-orbit coupling (notably the manganites, pyrochlores 
and spinels) led to a strong resurgence of interest in the 
AHE and eventually to deeper understanding. 

Since 2003 many systematics studies, both theoretical 
and experimental, have led to a better understanding of 
the AHE in the metallic regime, and to the recognition of 
new unexplored regimes that present challenges to future 
researchers. As it is often the case in condensed matter 
physics, attempts to understand this complex and fasci- 
nating phenomenon have motivated researchers to couple 
fundamental and sophisticated mathematical concepts to 
real- world materials issues. The aim of this review is to 
survey recent experimental progress in the field, and to 
present the theories in a systematic fashion. Researchers 
are now able to understand the links between different 
views on the AHE previously thought to be in conflict. 
Despite the progress in recent years, understanding is 
still incomplete. We highlight some intriguing questions 
that remain and speculate on the most promising avenues 
for future exploration. In this review we focus, in par- 
ticular, on reports that have contributed significantly to 
the modern view of the AHE. For previous reviews, the 
reader may consult Pugh (Pugh and Rostoker 1953) and 
Hurd (Hurd 19721. For more recent short overviews fo- 



cused on the topological aspects of the AHE, we point 



the induced M = xH, with x the susceptibihty ( |Fert| the reader to the rev iews by Nagaosa (|Nagaosa| |2006| ), 
However, in zero iJ, time-reversal in- and by Sinova et aZ.(|Sinova efa!] |2004[). A review of 



and Jaoul 



1972). 



variance (TRI) is spontaneously broken in the former, 
but not in the latter. Throughout the period 1960-1989, 
the two Hall effects were often regarded as a common 
phenomenon that should be understood microscopically 
on the same terms. It now seems clear that this view 
impeded progress. 

By the mid-1980s, interest in the AHE problem had 
waned significantly. The large body of Hall data gar- 
nered from experiments on dilute Kondo systems in the 
previous two decades showed that p^-y p and therefore 
appeared to favor the skew-scattering mechanism. The 
points of controversy remained unsettled, however, and 
the topic was still mired in confusion. 

Since the 1980's, the quantum Hall effect in two- 
dimensional (2D) electron systems in semiconductor het- 
erostructures has become a major field of research in 
physics (Prange and Girvin 1987). The accurate quan- 



the modern semiclassical treament of AHE was recently 
written by Sinitsyn (lSinitsyn| |2008| . The present review 
has been informed by ideas explained in the earlier works. 
Readers who are not familair with Berry phase concepts 
may find it u seful to consult the elementary review by 
Ong and Lee 



Ong and Lee 



mentary by MacDonald and Niu (MacDonald and Niu 



2006) and the popular com- 



20041 



tization of the Hall conductance is the hallmark of this 
phenomenon. Both the integer ( jThouless et al. 1982) 
and fractional quantum Hall effects can be explained in 
terms of the topological properties of the electronic wave- 
functions. For the case of electrons in a two-dimensional 
crystal, it has been found that the Hall conductance is 
connected to the topological integer (Chern number) de- 
fined for the Bloch wavcfunction over the first Brillouin 
zone ( |Thouless et al.\ 1982[ ). This way of thinking about 
the quantum Hall effect began to have a deep impact 



Some of the recent advances in the understanding of 
the AHE that will be covered in this review are: 

1. When (j-^^ is independent of (Txx, the AHE can 
often be understood in terms of the geometric con- 
cepts of Berry phase and Berry curvature in mo- 
mentum space. This AHE mechanism is respon- 
sible for the intrinsic AHE. In this regime, the 
anomalous Hall current can be thought of as the 
unquantized version of the quantum Hall effect. In 
2D systems the intrinsic AHE is quantized in units 
of e^//i at temperature T — Q when the Fermi level 
lies between Bloch state bands. 

2. Three broad regimes have been identified when sur- 
veying a large body of experimental data for di- 
verse materials: (i) A high conductivity regime 
{(^xx > 10^ (f2cm)^^) in which a linear contribution 



to a. 



AH 

xy 



Oxx due to skew scattering dominates 
In this regime the normal Hall conductiv- 
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ity contribution can be significant and even domi- 
nate axy, (ii) An intrinsic or scattering-independent 
regime in which is roughly independent of a^x 
(10" (r2cm)-i < < 10^ (f^cm)-i), (iii) A bad- 
metal regime {axx < lO'^ (ficm)"^) in which ax^ 
decreases with decreasing axx at a rate faster than 
linear. 

3. The relevance of the intrinsic mechanisms can be 
studied in-depth in magnetic materials with strong 
spin-orbit coupling, such as oxides and diluted mag- 
netic semiconductors (DMS). In these systems a 
systematic non-trivial comparison between the ob- 
served properties of systems with well controlled 
materials properties and theoretical model calcula- 
tions can be achieved. 

4. The role of band (anti)-crossings near the Fermi en- 
ergy has been identified using first-principles Berry 
curvature calculations as a mechanism which can 
lead to a large intrinsic AHE. 

5. Semiclassical treatment by a generalized Boltz- 
mann equation taking into account the Berry cur- 
vature and coherent inter-band mixing effects due 
to band structure and disorder has been formu- 
lated. This theory provides a clearer physical pic- 
ture of the AHE than early theories by identify- 
ing correctly all the semiclassically defined mech- 
anisms. This generalized semiclassical picture has 
been verified by comparison with controlled micro- 
scopic linear response treatments for identical mod- 
els. 

6. The relevance of non-coplanar spin structures with 
associated spin chirality and real-space Berry cur- 
vature to the AHE has been established both the- 
oretically and experimentally in several materials. 

7. Theoretical frameworks based on the Kubo formal- 
ism and the Keldysh formalism have been devel- 
oped which are capable of treating transport phe- 
nomena in systems with multiple bands. 

The review is aimed at experimentalists and theorists 
interested in the AHE. We have structured the review 
as follows. In the remainder of this section, we provide 
the minimal theoretical background necessary to under- 
stand the different AHE mechanisms. In particular we 
explain the scattering-independent Berry phase mech- 
anism which is more important for the AHE than for 
any other commonly measured transport coefficient. In 
Sec. [n] we review recent experimental results on a broad 
range of materials, and compare them with relevant cal- 
culations where available. In Sec. IIIII we discuss AHE 
theory from an historical perspective, explaining links 
between different ideas which are not always recognized, 
and discussing the physics behind some of the past con- 
fusion. The section may be skipped by readers who do 



the present understanding of the metallic theory based 
on a careful comparison of the different linear response 
theories which are now finally consistent. In Sec. [V| we 
present a summary and outlook. 



B. Parsing the AHE: 

The anomalous Hall effect is at its core a quantum phe- 
nomena which originates from quantum coherent band 
mixing effects by both the external electric field and 
the disorder potential. Like other coherent interference 
transport phenomena (e.g. weak localization), it cannot 
be satisfactorily explained using traditional semiclassical 
Boltzmann transport theory. Therefore, when parsing 
the different contributions to the AHE, they can be de- 
fined semiclassically only in a carefully elaborated theory. 

In this section we identify three distinct contributions 
which sum up to yield the full AHE: intrinsic, skew 
scattering, and side-jump contributions. We choose this 
nomenclature to reflect the modern literature without 
breaking completely from the established AHE lexicon 
(see Sec. 



III). However, unlike previous classifications. 



we base this parsing of the AHE on experimental and mi- 
croscopic transport theory considerations, rather than on 
the identification of one particular effect which could con- 
tribute to the AHE. The link to semiclassically defined 
processes is established after developing a fully general- 
ized Boltzmann transport theory which takes inter-band 
coherence effects into account and is fully equivalent to 



microscopic theories (Sec. IV. A I . In fact, much of the the 



oretical effort of the past few years has been expended 
in understanding this link between semiclassical and mi- 
croscopic theory which has escaped cohesion for a long 
time. 

A very natural classification of contributions to the 
AHE, which is guided by experiment and by microscopic 
theory of metals, is to separate them according to their 
dependence on the Bloch state transport lifetime r. In 
the theory, disorder is treated perturbatively and higher 
order terms vary with a higher power of the quasiparti- 
cle scattering rate r^^. As we will discuss, it is relatively 
easy to identify contributions to the anomalous Hall con- 
ductivity, ax^ , which vary as and as r". In experi- 
ment a similar separation can sometimes be achieved by 
plotting cTxy vs. the longitudinal conductivity axx oc r, 
when T is varied by altering disorder or varying tem- 
perature. More commonly (and equivalently) the Hall 
resistivity is separated into contributions proportional to 
Pxx and pIx. 

This partitioning seemingly gives only two contribu- 



tions to a^,^ , one ~ r and the other 



The first 



contribution we define as the skew-scattering contribu- 
Note that in this parsing of AHE con- 



tion, a 



AH- 

xy 



not wish to be burdened by history. Section IV discusses 



tributions it is the dependence on r (or axx) which de- 
fines it, not a particular mechanism linked to a micro- 
scopic or semiclassical theory. The second contribution 
proportional to r° (or independent of axx) we further 
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separate into two parts: intrinsic and side-jump. Al- 
though these two contributions cannot be separated ex- 
perimentally by dc measurements, they can be sepa- 
rated experimentally (as well as theoretically) by defin- 



ing the intrinsic contribution, a 



AH—int 
xy 



as the extrapo- 



lation of the ac-interband Hall conductivity to zero fre- 
quency in the limit of r — > oo, with 1/t faster 
than Lj — !■ 0. This then leaves a unique definition for 
the third and last contribution, termed side-jump, as 

AH—sj — lyAH _ AH — skew _ AH—int 
xy — xy xy xy 

We examine these three contributions below ( still at 
an introductory level). It is important to note that 
the above definitions have not relied on identifications 
of semiclassical processes such as side-jump scattering 



(Berger 19701 or skew-scattering from asymmetric con- 



tributions to the semiclassical scattering rates (Smit 
1955 ) identified in earlier theories. Not surprisingly, the 
contributions defined above contain these semiclassical 
processes. However, it is now understood (see Sec. IV I, 
that other contributions are present in the fully general- 
ized semiclassical theory which were not precisely identi- 
fied previously and which are necessary to be fully con- 
sitent with microscopic theories. 

The ideas explained briefly in this section are substan- 
tiated in Sec. [H] by analyses of tendencies in the AHE 
data of several different material classes, and in Sec. |III| 
and Sec. |IV| by an extensive technical discussion of AHE 
theory. We assume throughout that the ferromagnetic 
materials of interest are accurately described by a Stoner- 
like mean-field band theory. In applications to real ma- 
terials we imagine that the band theory is based on spin- 
density-functional theory (jJones and Gunnarsson 1989 ) 



with a local-spin-density or similar approximation for the 
exchange-correlation energy functional. 



1. Intrinsic contribution to 



Among the three contributions, the easiest to eval- 
uate accurately is the intrinsic contribution. We have 
defined the intrinsic contribution microscopically as the 
dc limit of the interband conductivity, a quantity which 
is not zero in ferromagnets when SOI are included. 
There is however a direct link to semiclassical theory in 
which the induced interband coherence is captured by 
a momentum-space Berry-phase related contribution to 
the anomalous velocity. We show this equivalence below. 

This contribution to the AHE was first derived by KL 
(Karplus and Luttinger |1954) but its topological nature 



was not fully appreciated unt i l rece ntly (Jungwirth et al. 



2002b 


Onoda and Nagaosa 


2002 


1. 1 


wirth 


et al. (Jungwirth et al. 


2002b') 



The work of Jung- 



the experimental importance of the AHE in ferromag- 
netic semiconductors and also by the thorough earlier 
analysis of the relationship between momentum space 
Berry phases and anomalous ve locities in semiclassical 
transport theory by Niu e t al. (Chang and Niu 1996 



Sundaram and Niu 19991. The frequency-dependent 



inter-band Hall conductivity, which reduces to the intrin- 
sic anomalous Hall conductivity in the dc limit, had been 
evaluated earlier for a number of materials by Mainkar 
et al. (Mainkar et al. 19961 and Guo and Ebert (Quo 



and Ebert[ 1995 1 but the topological connection was not 
recognized. 

The intrinsic contribution to the conductivity is de- 
pendent only on the band structure of the perfect crys- 
tal, hence its name. It can be calculated directly from 
the simple Kubo formula for the Hall conductivity for an 
ideal lattice, given the eigenstates \n, k) and eigenvalues 
£„(fc) of a Bloch Hamiltonian H: 



_AH—int 



dk 



xlm 



n, k\vi{k)\n' , k){n' , k\vj{k)\n, k) 
(e„(fc)-e„,(fc))2 ■ 

(1.2) 



In Eq. (1.2) H is the fe-dependent Hamiltonian for the 
periodic part of the Bloch functions and the velocity op- 
erator is defined by 



v{k) = ^[r,H{k)]^lVkHik). 
in n 



(1.3) 



(1.21 



Note the restriction n' in Eq 

What makes this contribution quite unique is that, like 
the quantum Hall effect in a crystal, it is directly linked 
to t he top ological properties of the Bloch states. (See 
Sec. III.B[ ) Specifically it is proportional to the integra- 
tion over the Fermi sea of the Berry's curvature of each 
occupied band, or equivalently (Haldane 2004 Wang 



et al. 20071 to the integral of Berry phases over cuts 
of Fermi surface segments. This result can be derived by 
noting that 

, . 1^ , / , X {n,k\VkH(k)\n',k) ,^ 
e„'(fc) - e„(fe) 



Using this expression, Eq. (1.2 1 reduces to 



^AH—int 



dk 

(2^ 



/(£„(fc)) 6l(fc), (1.5) 



where e^-^ is the anti-symmetric tensor, a„(fc) is the 
Berry-phase connection a„{k) — i(n, fe| V^jn, fc), and 
b„(fc) the Berry-phase curvature 



b„(fc) = Vfc X a„(fe) 



(1.6) 



corresponding to the states {\n,k)}. 

This same linear response contribution to the AHE 
conductivity can be obtained from the semiclassical the- 



ory of wave-packets dynamics (Chang and Niu 1996 
Marder[ |2000[ [Sundaram and Niu[ |1999[ ). It can be 



shown that the wavepacket group velocity has an addi- 
tional contribution in the presence of an electric field: 
r^ = dEn{k)/hdk~{E/h)xbn{k). (See Sec. [TO]) The 
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intrinsic Hall conductivity formula, Eq. ( 1.5 1, is obtained 
simply by summing the second (anomalous) term over all 
occupied states. 

One of the motivations for identifying the intrinsic 
contribution (7 mi IS that it can be evaluated accu- 

J. y 

rately even for relatively complex materials using first- 
principles electronic structure theory techniques. In 
many materials which have strongly spin-orbit coupled 
bands, the intrinsic contribution seems to dominates the 
AHE. 



2. Skew scattering contribution to 

The skew scattering contribution to the AHE can be 
sharply defined; it is simply the contribution which is 
proportional to the Bloch state transport lifetime. It will 
therefore tend to dominate in nearly perfect crystals. It is 
the only contribution to the AHE which appears within 
the confines of traditional Boltzmann transport theory 
in which interband coherence effects are completely ne- 
glected. Skew scattering is due to chiral features which 
appear in the disorder scattering of spin-orbit coupled 
ferromagnets. This mechanism was first identified by 



Smit (Smit 1955 19581 



Treatments of semi-classical Boltzmann transport the- 
ory found in textbooks often appeal to the principle of 
detailed balance which states that the transition prob- 
ability Wn^m from n to m is identical to the transi- 
tion probabihty in the opposite direction (Wm^„). Al- 
though these two transition probabilities are identical 
in a Fermi's golden-rule approximation, since Wn—>n' — 
{2-K/h)\{n\V\n')\^5{En - E^'), where V is the perturba- 
tion inducing the transition, detailed balance in this mi- 
croscopic sense is not generic. In the presence of spin- 
orbit coupling, either in the Hamiltonian of the perfect 
crystal or in the disorder Hamiltonian, a transition which 
is right-handed with respect to the magnetization di- 
rection has a different transition probability than the 
corresponding left-handed transition. When the transi- 
tion rates are evaluated perturbatively, asymmetric chi- 
ral contributions appear first at third order. (See Sec. 
rV.Al. In simple models the asymmetric chiral contribu- 



tion to the transition probability is often assummed to 
have the form (see Sec. III.C.2.al: 



Wkk' = -Tl^fe X k' ■ M, 



(1.7) 



When this asymmetry is inserted into the Boltzmann 
equation it leads to a current proportional to the lon- 
gitudinal current driven by E and perpendicular to both 
E and Mg. When this mechanism dominates, both the 
Hall conductivity cth and the conductivity a are propor- 
tional to the transport lifetime t and the Hall resistivity 
pskew _ ^skew p2 therefore proportional to the longitu- 
dinal resistivity p. 

There are several specific mechanisms for skew scat- 
tering (see Sec III.C.2 and Sec. 



IV.Al. Evaluation of 



or resistivity requires simply that the conventional lin- 
earized Boltzmann equation be solved using a collision 
term with accurate transition probabilities, since these 
will generically include a chiral contribution. In prac- 
tice our ability to accurately estimate the skew scatter- 
ing contribution to the AHE of a real material is lim- 
ited only by typically imperfect characterization of its 
disorder. We emphasize that skew scattering contribu- 
tions to uh are present not only because of spin-orbit 
coupling in the disorder Hamiltonian, but also because 
of spin-orbit coupling in the perfect crystal Hamiltonain 
combined with purely scalar disorder. Either source of 
skew-scattering could dominate a^^~'''"^^ depending on 
the host material and also on the type of impurities. 

We end this subsection with a small note directed to 
the reader who is more versed in the latest development 
of the full semiclassical theory of the AHE and in its 



comparison to the microscopic theory (see Sec. IV. A and 
IV.B.2). We have been careful above not to define the 



skew-scattering contribution to the AHE as the sum of 
all the contributions arising from the asymmetric scatter- 
ing rate present in the collision term of the Boltzmann 
transport equation. We know from microscopic theory 
that this asymmetry also makes an AHE contribution or 
order r°. There exists a contribution from this asymme- 
try which is actually present in the microscopic theory 
treatment associated with the so called ladder diagram 
corrections to the conductivity, and therefore of order r° . 
In our experimentally practical parsing of AHE contri- 
butions we do not associate this contribution with skew- 
scattering but place it under the umbrella of side-jump 
scattering even though it does not physically originate 
from any side-step type of scattering. 



3. Side-jump contribution to a^jf^ 

Given the sharp defintions we have provided for the 
intrinsic and skew scattering contributions to the AHE 
conductivity, the equation 



AH 



AH-int 



^AH— skew 



AH-sj 
'xy 



(1.8) 



defines the side-jump contribution as the difference be- 
tween the full Hall conductivity and the two-simpler con- 
tributions. In using the term side-jump for the remaining 
contribution, we are appealing to the historically estab- 
lished taxonomy outlined in the previous section. Estab- 
lishing this connection mathematically has been the most 
controversial aspects of AHE theory, and the one which 
has taken the longest to clarify from a theory point of 
view. Although this classification of Hall conductivity 
contributions is often useful (see below), it is not gener- 
ically true that the only correction to the intrinsic and 
skew contributions can be physically identified with the 
side-jump process defined as in the earlier studies of the 
AHE (lBerger[p64l. 



the skew scattering contribution to the Hall conductivity 



The basic semiclassical argument for a side-jump 
contribution can be stated straight-forwardly: when 
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considering the scattering of a Gaussian wavepacket 
from a spherical impurity with SOI ( Hso = 
{l/2ni'^c^){r^^dV/dr)SzLz), a wavepacket with incident 
wave-vector k will suffer a displacement transverse to k 
equal to ^kfi^ /m^c^ . This type of contribution was first 



noticed, but discarded, by Smit (Smit 19581 and reintro- 
duced by Berger ( |Berger 1964[ ) who argued that it was 
the key contribution to the AHE. This kind of mecha- 
nism clearly lies outside the bounds of traditional Boltz- 
mann transport theory in which only the probabilities of 
transitions between Bloch states appears, and not micro- 
scopic details of the scattering processes. This contribu- 
tion to the conductivity ends up being independent of r 
and therefore contributes to the AHE at the same order 
as the intrinsic contribution in an expansion in powers 
of scattering rate. The separation between intrinsic and 
side-jump contributions, which cannot be distinquished 
by their dependence on t, has been perhaps the most 
argued aspect of AHE theory since they cannot be dis- 
tinquished by their dependence on scattering rate (see 
Sec.|IH.C.2.d|). 



As explained clearly in a recent review by Sinitsyn 



(Sinitsyn 20081, side-jump and intrinsic contributions 
have quite different dependences on more specific sys- 
tem parameters, particularly in systems with complex 
band structures. Some of the initial controversy which 
surrounded side jump theories was associated with phys- 
ical meaning ascribed to quantities which were plainly 
gauge dependent, like the Berry's connection which in 
early theories is typically identified as the definition of 
the side-step upon scattering. Studies of simple models, 
for example models of semiconductor conduction bands, 
also gave results in which the side-jump contribution 
seemed to be the same size but opposite in sign com- 



pared to the intrinsic contribution ( Nozieres and Lewiner 
1973[ ). We now understand ( [Sinitsyn et al.\ |2007p that 
these cancellations are unlikely, except in models with a 
very simple band structure, e.g. one with a constant 
Berry's curvature. It is only through comparison be- 
tween fully microscopic linear response theory calcula- 
tions, based on equivalently valid microscopic formalisms 
such as Keldysh (non-equilibrium Grenn's function) or 
Kubo formalisms, and the systematically developed semi- 
classical theory that the specific contribution due to the 
side-jump mechanism can be separately identified with 
confidence (see Sec. 



IV. A I 



Having said this, all the calculations comparing the 
intrinsic and side-jump contibutions to the AHE from a 
microscopic point of view have been performed for very 
simple models not immediately linked to real materials. 
A practical approch which is followed at present for ma- 
terials in which a"^^ seems to be independent of (t^x, is 
to first calculate the intrisic contribution to the AHE. If 
this explains the observation (and it appears that it usu- 
ally does), then it is deemed that the intrinsic mechanism 
dominates. If not, we can take some comfort from under- 
standing on the basis of simple model results, that there 
can be other contributions to which are also inde- 



pendent of axx and can for the most part be identified 
with the side jump mechanism. Unfortunately it seems 
extremelly challenging, if not impossible, to develop a 
predictive theory for these contributions, partly because 
they require many higher order terms in the perturbation 
theory that be summed, but more fundamentally because 
they depend sensitively on details of the disorder in a 
particular material which are normally unknown. 



II. EXPERIMENTAL AND THEORETICAL STUDIES ON 
SPECIFIC MATERIALS 



A. Transition-metals 



1. Early experiments 



Four decades after the discovery of the AHE, an empir- 
ical relation between magnetization and Hall resistivity 
was proposed independently by A. W. Smith and by E. 



M. Pugh ( 


Pugh 


and Sears 


1929 



1930 



Pugh and Lippert 1932 Smith 



AHE in Fe, Ni and Co and the alloys Co-Ni and Ni-Cu 
in magnetic fields up to 17 kG over large intervals in T 
(10-800 K in the case of Ni), and found that the Hall 
resistivity pn is comprised of 2 terms, viz. 



PH = RoH + RiM{T,H), 



(2.1) 



where M{T, H) is the magnetization averaged over the 
sample. Pugh defined R^y and i?i as the ordinary and 
extra-ordinary Hall coefficients, respectively. The latter 
Ri = Rs is now called the anomalous Hall coefficient (as 
in Eq.[rT|). 

On dividing Eq. (2.1 1 by p^, we see that it just ex- 



presses the additivity of the Hall currents: the total Hall 
conductivity ct*°* equals <t^J^ + , where is the 
ordinary Hall conductivity and is the AHE conduc- 
tivity. A second implication of Eq. (2.1 1 emerges when 
we consider the role of domains. The anomalous Hall co- 
efficient in Eq.(2.1) is proportional to the AHE in a sin- 
gle domain. As iJ — > 0, proliferation of domains rapidly 
reduces M{T,H) to zero (we ignore pinning). Cancella- 
tions of a^.^ between domains result in a zero net Hall 
current. Hence the observed AHE term mimics the field 
profile of M (T, H), as implied by Pugh's term RiM. The 
role of H is simply to align the AHE currents by rotat- 
ing the domains into alignment. The Lorentz force term 
c^xy^ is a "background" current with no bearing on the 
AHE problem. 

The most interesting implication of Eq. (2.1 1 is that. 



in the absence of if, a single domain engenders a spon- 
taneous Hall current transverse to both M and E. Un- 
derstanding the origin of this spontaneous off-diagonal 
current has been a fundamental problem of charge trans- 
port in solids for the past 60 years. The AHE is also 
called the spontaneous Hall effect and the extraordinary 
Hall effect in the older literature. 
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2. Recent experiments 

The resurgence of interest in the AHE motivated by 
the Berry-phase approach (Sec. I.B.I I has led to many 



new Hall experiments on 3c? transition metals and their 
oxides. Both the recent and the older literature on Fe 
and Fe304 are reviewed in this section. An important 
finding of these studies is the emergence of three distinct 
regimes roughly delimited by the conductivity u^x and 
characterized by the dependence of ax^ on Uxx- The 
three regimes are: 

i) A high conductivity regime for axx ^ 10^ (ficm)^^ 
in which a^^^'''^'^^ ~ dominates (Tx^ , 

ii) A good metal regime for cr^^ ~ 10'' — 10^ (ficm)"^ 
in which axy ~ cr^^,, 

iii) A bad metal/hopping regime for axx < 10* 
(f]cm)-i in which a^y^ - crl-^'^-^. 

We discuss each of these regimes below. 

High conductivity regime - The Hall conductivity in the 
high-purity regime, axx > 0.5 x 10® (ilcm)^^, is dom- 
inated by the skew scattering contribution cr^^'^™. The 
high-purity regime is one of the least studied experimen- 
tally. This regime is very challenging to investigate ex- 
perimentally because the field H required for saturating 
M also yields a very large ordinary Hall effect (OHE) 
and Rq tends to be of the order of Rs ( Schad et al.\ 
1998[). In the limit lUcT > 1, the OHE conductivity 



a"y^ may be nonlinear in H (lOc is the cyclotron fre- 
quency). Although axy™ increases as r, the OHE term 
axy^ increases as and therefore the latter ultimately 
dominates, and the AHE current may be unresolvable. 
Even though the anomalous Hall current can not always 
be cleanly separated from the normal Lorentz-force Hall 
effect in the high conductivity regime, the total Hall cur- 
rent invariably increases with axx in a way which provides 
compelling evidence for a skew-scattering contribution. 

In spite of these challenges several studies have man- 
aged to convincingly separate the competing contribu- 
tions and have identifed a dominant linear rel ation be- 
tween a^^ and axx for axx ^ 10® (ficm) ~ ^ ( Majum- 
dar and Berger 1973 Shiomi et al.\ 2009). In an early 



study Majumdar et al. (IMajumdar and Bergerl 119731) 



grew highly pure Fe doped with Co. The resulting ax" , 
obtained from Kohler plots extrapolation to zero field, 
show a clear dependence of cr^^^ ~ c^^x (Fig. |4]a). In a 
more recent study, a similar finding (linear dependence of 



^xx 

) was observed by Shiomi et al. 
20091) in Fe doped with Co, Mn, Cr, and S: 



'.. ( |Shiomi et aL^ 
i. In these stud- 



ies the high temperature contribution to axy (presumed 
to be intrinsic plus side jump) was substracted from axy 
and a linear dependence of the resulting cr^jf^ is observed 
(Fig. |4] a and b). In this recent study the conductivity 
is intentionally reduced by impurity doping to find the 
linear region and reliably exclude the Lorernz contribu- 
tion. The results of these authors show, in particular. 



that the slope of axx vs. ax^ depends on the species 
of the impurities as it is expected in the regime domi- 
nated by skew scattering. It is reassuring to note that 
the skewness parameters {Sskew — cr"^^ /o'xx) implied by 
the older and the more recent experiments are consistent, 
in spite of differences in the conductivity ranges stud- 
ied. Sskew is independent of axx (Majumdar and Berger 



1973 Shiomi et al. 2009 ) as it should be when the skew 
scattering mechanism dominates. Further experiments 
in this regime are desirable to fully investiage the dif- 
ferent dependence on doping, temperature, and impurity 
type. Also, new approaches to reliably disentangle the 
AHE and OHE currents will be needed to faicilitate such 
studies. 

Good metal regime - Experiments re-examining the AHE 
in Fe, Ni, and Co have b een performed by Miyasato et 
al. (Miyasato et al. 2007). These experiments indicate a 
regime of axy versus axx in which axy is insensitive to axx 
in the range axx ~ 10''-10® (f]cm)-i (see Fig. [5|. This 
suggests that the scattering independent mechanisms (in- 
trinsic and side-jump) dominates in this regime. How- 
ever, in comparing this phenomenology to the discussion 
of AHE mechanisms in Sec. |I.B| one must keep in mind 
that the temperature has been varied in the Hall data 
on Fe, Ni, and Co in order to change the resistivity, even 
though it is restricted to the range well below Tc (Fi g. [5 
upper panel). In the mechanisms discussed in Sec. |I.B 
only elastic scattering was taken into account. Earlier 
tests of the dependence of pxy carried by varying T 
were treated as suspect in the early AHE period (see Fig. 
|2]) because the role of inelastic scattering was not fully un- 
derstood. The effect of inelastic scattering from phonons 
and spin waves remains open in AHE theory and is not 
addressed in this review. 

Bad metal/hopping regime- Several groups have mea- 



sured axy in Fe and Fe304 thin-film ferromagnets (Feng 



et aL| |1975l JFernandez-Pacheco al.\ |2008| [Miyasato 



et a?.[ |2007| [Sangiao et al.\ 2009f Venkateshvaran et al._ 
2008[ ). (See Fig. [6|) Sangiao et al. (ISang iao et aZ.|[2009 ) 
studied epitaxial thin-films of Fe deposited by pulsed- 
laser deposition (PLD) on single-crystal MgO (001) sub- 
strates at pressures < 5 x 10~^ Torr. To vary axx over 
a broad range, they varied the film thickness t from 1 to 
10 nm. The p vs. T profile for the film with t = 1.8 nm 
displays a resistance minimum near 50 K, below which p 
shows an upturn which has been ascribed to localization 
or electron interaction effects (Fig. [t]). The magnetiza- 
tion Al is nominally unchanged from the bulk value (ex- 
cept possibly in the 1.3 nm film). AHE experiments were 
carried out from 2-300 K on these films and displayed as 
axy vs. axx plots together with previously published re- 
sults (Fig. [6|. In the plot, the AHE data from films with 
t < 2 nm fall in the weakly localized regime. The com- 
bined plot shows that Sangio et al.'s data are coUinear (on 
a logarithmic scale) with those measured on 1 ^m-thick 
films by Miyasato et al. ( Miyasato et al. 2007 ) . For the 
three samples with t = 1.3-2 nm, the inferred exponent in 
the dirty regime is on average 1.66. A concern is that 
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♦ Fe-Co alloy (Majumdar 1973) 
O Fe-Co alloy (Shiomi 09) 
« Fe-Cr alloy (Shiomi 09) 
M Fe-Mn alloy (Shiomi 09) 
k- Fe-Si alloy (Shiomi 09) 
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FIG. 4 cr^jf for pure Fe film doped with Cr, Co, Mn and Si 
vs. a^x at low temperatures {T = 4.2 K and T = 5.5 K ) (a). 
In many of the alloys, particularly in the Co doped system, 
the linear scaling in the higher conductivy sector, a^x > 10^ 
(Slcm)"^, implies that skew scattering dominates . After 
Ref. jMajumdar and Berger|[T9 73 and Ref. jShiomi et aL\^m9\ 



The data from ( Shiomi et al.^ ,2009) , shown also at a larger 



scale in (b), is obtained by substracting the high temperature 
contribuiton to axy In the data shown the ordinary Hall 
contributi on has been identifie d and substracted. [Panel (b) 
From Ref. 



Shiomi et al. 2009 



the data from the 1.3 nm film was obtained by subtract- 
ing a InT term from p (the subtraction procedure was 
not described). How locahzation affects the scahng plot 
is an open issue at present. In Sec. |II.E| we discuss re- 
cent AHE measurements in disordered polycrystalline Fe 



films with i < 10 nm by Mitra et al. ( [Mitra et aL[|2007 l 

Recent progress in understanding weak-localization cor- 
rections to the AHE is also reviewed there. 



In magnetite, Fe304, scaling of ax 



V 



a XX with 

a 1.6 — 1.8 was already apparent in early experi- 
ments on polycrystalline samples (Feng et al. 19751. 






10 10 
^'xx cm"'] 
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FIG. 5 Measurements of the Hall conductivity and resistivity 
in single-crystal Fe and in thin foils of Fe, Co, and Ni. The 
top and lower right panels show the T dependence of a^y and 
Pxx, respectively. The lower left panel plots \a-xy\ against axx- 
After Ref. jMiyasato eTal] [20071 
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FIG. 6 Combined plot of the AHE conductivity \axy\ versus 
the conductivity axx in epitaxial films of Fe grown on MgO 
with thickness t — 2.5, 2.0, 1.8 and 1. 3 nm ([Sangiao et al. 
2009[), in polycrystalline Fe3_2,Zn^04 ([Feng et a/.| I1975^, m 



thin-film Fe3-a;Zna;04 between 90 and 350 K ( Venkatcshvaran 
et al.\ | 2008|l a nd above the Verwey transition ( [Fernandez- 



Pacheco et 



20081 



Recently, two groups have re-investigated the AHE in 
epitaxial thin fil ms (data incl uded in Fig. [6| . Fernandez- 



Pacheco et al. ( Fernandez-Pachcco et al. [2008| mea- 



sured a series of thin-film samples of Fe304 grown by 
PLD on MgO (001) substrates in ultra-high vacuum, 
whereas Venkatcshvaran et al. ( Venkateshvaran et al. 



|2008| ) studied both pure Fe304 and Zn-doped magnetite 
Fe3-x'Zinx04 deposited on MgO and AI2O3 substrates 
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FIG. 7 The T dependence of p in epitaxial thin-films 
MgO(001)/Fe(t)/MgO with t = 1.8 nm and 2.5 nm. The 
inset shows how To, the temperature of the resistivity mini- 
mum, varies with t. [From Ref. Sangiao et al. 2009 ] 



FIG. 8 Longitudinal resistivity pxx vs. T for epitaxial 
Fe3_^Zn^04 films. The (001), (110) and (111) oriented films 
were grow n on MgO(OOl), MgO(llO) a nd AI2O3 substrates. 
[From Ref. Venkateshvaran et al. 2008 ] 



grown by laser molecular-beam epitaxy under pure Ar 
or Ar/0 mixture. In both studies, p increases monotoni- 
cally by a factor of ~10 as T decreases from 300 K to the 
Verwey transition temperature Ty = 120 K. Below Ty, 
p further increases by a factor of 10 to 100. The results 
for p vs. T from Venkateshvaran et al. (Venkateshvaran 
et al.\ |2008[ ) is shown in Fig. [sj 



The large values of p and its insulating trend imply 
that magnetite falls in the strongly localized regime, in 
contrast to thin-film Fe which lies partly in the weak- 
localization (or incoherent) regime. 

Both groups find good scaling fits extending over sev- 
eral decades of a^x with a ^ 1.6 — 1.8 when varying 



T. Fernandez-Pacheco et al. ( Fernandez-Pacheco et al. 



2008[ ) plot a^y vs. a^x in the range 150< T <300 K 
for several thicknesses t and infer an exponent a = 1.6. 
Venkateshvaran et al. ( Venkateshvaran et al. 2008 1 plot 



axy vs. axx from 90 to 350 K and obtain power-law fits 
with a — 1.69, in both pure and Zn-doped magnetite 
(data shown in Fig. [6]). Significantly, the 2 groups find 
that a is unchanged below Ty. There is presently no the- 
ory in the poorly conducting regime which predicts the 
observed scaling {axx < 10~^ {Qcm)~^). 



3. Comparison to theories 

Detailed first-principles calculations of the intrinsic 
contribution to the AHE conductivity have been per- 
formed for bcc Fe ( |Wang et a{] '2006': Y ao et a?.[|2004| , 
fee Ni, and hpc Co ( [Wang et al , 2007J . In Fe and Co, 
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FIG. 9 First-principles calculation of the band dispersions 
and the Berry-phase curvature summed over occupied bands. 
[From Ref. |Yao et al.\ |2004|] 



enhanced near avoided crossings, opposite direction for 
the upper and lower bands. A large contribution to cr™* 
when the crossing is at the Fermi energy so that only 
one of the two bands is occupied, e.g. near the point H 
in Fig. [9] A map showing the contributions of different 



regions of the FS to b^{k) is shown in Fig. [lOl The SOI 
can lift an accidental degeneracy at certain wavevectors 
k. These points act as a magnetic monopole for the Berry 

2003) . In the parame- 



curvature in fc-space (Fang et al. 



ter space of spin-orbit coupling, axy is nonpertur bative in 
nature. The effect of these "parity anomalies" (Jackiw 



of topological enhancement in the intrinsic AHE was re- 
ported for a tight-binding calculation with the 2 orbitals 



dzx and dy^ on a square lattice (Kontani et al. 20071. 



19841 on the Hall conductivity was first discussed by Hal- 
the values of axy inferred from the Berry curvature ri^(fc) dane ( |Haldane[ |1988[ ). A different conclusion on the role 
are 7.5 x 10^ and 4.8 x 10^ (ficm)"^, respectively, in rea- 
sonable agreement with experiment. In Ni, however, the 
calculated value —2.2 x 10^ (ficm)^^ is only 30% of the 
experimental value. 

These calculations uncovered the crucial role played 
by avoided-crossings of band dispersions near the Fermi 
energy ep. The Berry curvature is always strongly 



Motivated by the enhancement at the crossing points 
discussed above, Onoda oZ. (Onodae<aZ. 2006a 20081 



proposed a minimal model that focuses on the topologi- 
cal and resonantly enhanced nature of the intrinsic AHE 
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FIG. 10 First-principles calculation of the FS in the (010) 
plane (solid lines) and the Berry curvature in atomic units 
(color map) . [From Ref. |Yao et al] |2004[] 
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arising from the sharp peak in 6^(fe) near avoided cross- 
ings. The minimal model is essentially a 2D Rashba 



model (Bychkov and Rashba 19841 with an exchange 



oc r, it is suppressed 



an in- 



field which breaks symmetries and accounts for the mag- 
netic order and a random impurity potential to account 
for disorder. The model is discussed in Sec. IIV.DI In 
the clean limit, a^y is dominated by the extrinsic skew- 
scattering contribution cr^^'^^, which almost masks the in- 
trinsic contribution cr"^ . Since cr^^°^ 
by increased impurity scattering, whereas crjfjj 
terband effect - is unaffected. In the moderately dirty 
regime where the quasiparticle damping is larger than 
the energy splitting at the avoided crossing (typically, 
the SOI energy) but less than the bandwidth, a^^^ domi- 
nates 0"!^°™ . As a result , one expects a crossover from the 
extrinsic to the intrinsic regime. When skew scattering 
is due to a spin-dependent scattering potential instead of 
spin-orbit coupling in the Bloch states, the skew to intrin- 
sic crossover could be controlled by a different condition. 
In this minimal model, a well-defined plateau is not well- 
reproduced in the intrinsic regime unless a weak impurity 
potential is assumed ( Onoda et al. 2008 1 . This crossover 



may be seen when the skew-scattering term shares the 



20091 



same sign as the intrinsic one (Kovalev et al 

With further increase in the scattering strength, spec 
tral broadening leads to t he scaling relat i onship (7xv cx: 



a^-^, as discussed above (Kovalev et al. 2009 



Onoda 



et al. 2006b 2008). In the strong-disorder regime, a^x 
is no longer linear in the scattering lifetime r. A dif- 
ferent scaling, axy oc ct^^, attributed to broadening of 
the electronic spectrum in the intrinsic regime, has been 
proposed by Kontani et al. (Kontani et al. 2007). 



As discussed above, there is some experimental evi- 
dence that this scaling prevails not only in the dirty 
metallic regime, but also deep into the hopping regime. 



FIG. 11 Anomalous HaU effect in SrRuOa. (A) The mag- 
netization M, (B) longitudinal resistivity pxx, and (C) trans- 
verse resistivity pxy as functions of the temperature T for the 
single crystal and thin- film SrRuOa, as well as for Ca-doped 
Sro.8Cao.2Ru03 thin film. fiB is the Bohr magneton. [From 
Ref. 



Fang et al. 2003 



Sangio et al. (jSangiao et al. 
exponent a 
regime 



2009), for e.g., obtained the 



1.7 in epitaxial thin- film Fe in the dirty 
In manganite, scaling seems to hold, with the 
same nominal value of a, even below the Verwey tran- 
sition where charge transport is deep in the hopping 



regime ( Fernandez-Pacheco et al. 2008 Venkateshvaran 



et al. 2008[ ) . These regimes are well beyond the purview 
of either the minimal model, which considers only elas- 
tic scattering, or the theoretical approximations used to 



2006b 20081 



model its properties (Onoda et al 

Nonetheless, the experimental reports have uncovered 
a robust scaling relationship with a near 1.6, which ex- 
tends over a remarkably large range of axx- The origin 
of this scaling is an open issue at present. 



B. Complex oxide ferromagnets 

1. First-principles calculations and experiments on SrRuOs 

The perovskite oxide SrRuOs is an itinerant ferromag- 
net with a critical temperature Tc of 165 K. The electrons 
occupying the 4d t2g orbitals in Ru*'^ have a SOI en- 
ergy much larger than that for 3d electrons. Early trans- 
port investigations of this material were reported in Refs. 
Allen et al.\ 119961 and llzumi et al. \ 119971 The latter au- 



thors also reported results on thin- film SrTiOa. Recently, 
the Berry-phase theory has been applied to account for its 
AHE (Fang et al.\ |2003| [Mathieu et'al] |2004a|b[ ) , which 
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FIG. 12 The calculated transverse conductivity a^y as a 
function of the chemical potential ^ for SrRuOa. The chaotic 
behavior is the fingerprint of the Berry curvature distribution 
illustrated in Fig. [l3] [From Ref. |Fang et al.\ [2003)] 





FIG. 14 (Upper panel) Combined plots of the Hall conduc- 
tivity (Txy vs. magnetization M in 5 samples of the ruthen- 
ate Sri-xCaxRuOa (0 < a; < 0.4). The inset compares data 
Oxy at a; = (triangles) with calculated values (solid curve). 
(Lower panel) First-principles calculations of a^y vs. M for 
cubic and orthorhombic structures. The effect of broaden- 
ing on the curves is shown for the orthorhombic case. [From 
ReL iMathieu et al.\ |2004b|l 




FIG. 13 The Berry curvature b^(k) for a band as a function 
of fex = (kx,ky) with the fixed fez = 0. [From Ref. Fang 



et al. 2003 



is strongly T dependent (FigjTTJ:). Neither the KL the- 
ory nor the skew-scattering theory seemed adequate for 
explaining the T dependence of the inferred AHE con- 
ductivity axy (Fang et al. 



20031 



The experimental results motivated a detailed first- 
principles, band-structure calculation that fully incorpo- 
rated the SOI. The AHE conductivity a^y w as calculated 



directly using the Kubo formula Eq. (1.2 1 ( Fang et a/.) 
[2003). To handle numerical instabilities which arise near 
certain critical points, a fictitious energy broadening S 
— 70 meV was introduced in the energy denominator. 
Fig. 12 shows the dependence of axyifJ,) on the chemical 
potential /i. In sharp contrast to the diagonal conduc- 
tivity a^x, o'xyiti) fluctuates strongly, displaying sharp 
peaks and numerous changes in sign. The fluctuations 
may be understood if we map the momentum dependence 
of the Berry curvature bz{k) in the occupied band. For 
example. Fig. 13 displays &z(fc) plotted as a function of 
fc_L — (kx,ky), with kz flxed at 0. The prominent peak 
at fc_L = corresponds to the avoided crossing of the 
energy band dispersions, which are split by the SOI. As 
discussed in Sec. |I.B| variation of the exchange splitting 
caused by a change in the spontaneous magnetization M 
strongly affects Uxy in a nontrivial way. 



From the flrst-principles calculations, one may esti- 
mate the temperature dependence of Uxy by assuming 
that it is due to the temperature-dependence of the Bloch 
state exchange splitting and that this splitting is propor- 
tional to the temperature-dependent magnetization. The 
new insight is that the T-dependence of (Txy(T) simply 
reflects the M dependence of axy- at a finite tempera- 
ture T', the magnitude and sign of (Jxy may be deduced 
by using the value of M(T') in the zero-T curve. This 
proposal was tested against the results on both the pure 
material and the Ca-doped material Sri_a;Caa;Ru03. In 
the latter, Ca doping suppresses both Tc and M system- 
atically (Mathieu et al. 2004a|b ) . As shown in Fig. 14 
(upper panel), the measured values of M and axy, ob- 
tained from 5 samples with Ca content 0.4< a: < 0, fall on 
2 continuous curves. In the inset, the curve for the pure 
sample {x = 0) is compared with the calculation. The 



lower panel of Fig. 14 compares calculated curves of axy 
for cubic and orthorhombic lattice structures. The sen- 
sitivity of axy to the lattice symmetry reflects the dom- 
inant contribution of the avoided crossing near ep. The 
sensitivity to broadening is shown for the orthorhombic 



Kats et al. (Kats et al. 20041 also have studied the 



magnetic field dependence of pxy in an epitaxial film of 
SrRuOs. They have observed sign-changes in pxy near a 
magnetic field B = 3 T at T = 130 K and near B = 8 
T at 134 K. This seems to be qualitative consistent with 
the Berry-phase scenario. On the other hand the authors 
suggest that the intrinsic-dominated picture is likely in- 
complete (or incorrect) near Tc (Kats et al. 20041. 
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2. Spin chirality mechanism of the AHE in manganites 



In the manganites, e.g. Lai_xCaxMn03 (LCMO), the 
three t2g electrons on each Mn ion form a core local mo- 
ment of spin S — ^. A large Hund energy Jh aligns the 
core spin S with the s — ^ spin of an itinerant electron 
that momentarily occupies the eg orbital. Because this 
Hund coupling leads to an extraordinary magnetoresis- 
tance in weak iJ, the manganites are called colossal mag- 



netoresistance (CMR) materials (Tokura and Tomioka 



1999 



Eq. 



The double exchange theory summarized by 



( 2.2[ ) (below) is widely adopted to describe the onset 
of ferromagnetism in the CMR manganites. 

As T decreases below the Curie temperature Tc — 270 
K in LCMO, the resistivity p falls rapidly from ~15 mfi 
cm to metallic values <2 mficm. CMR is observed over a 
significant interval of temperatures above and below Tc, 
where charge transport occurs by hopping of electrons 
between adjacent Mn ions (Fig. 15 1). At each Mn site 
i, the Hund energy tends to align the carrier spin s with 
the core spin Si. 



Early theories of hopping conductivity (Holstein 1961 ) 
predicted the existence of a Hall current produced by the 
phase shift (Peierls factor) associated with the magnetic 
flux (f) piercing the area defined by 3 non-coUinear atoms. 
However, the hopping Hall current is weak. The observa- 
tion of a large p^y in LCMO that attains a broad maxi- 



mum in modest H (Fig. 15 d) led Matl et at. (Matl et at. 



1998 ) to propose that the phase shift is geometric in ori- 
gin, arising from the solid angle described by s as the 
electron visits each Mn site (s||<S^i at each site i as shown 
16|. To obtain the large pxy seen, one requires 

Since Si 



in Fig 

Si to cleEne a finite solid angle Q. 



gradually 

aligns with H with increasing field, this effect should dis- 
appear along with (f2), as observed in the experiment. 
This appears to be the first application of a geometric- 
phase mechanism to account for an AHE experiment. 
Subsequently, Ye et al. (Ye et al. 1999) considered 



the Berry phase due to the thermal excitations of the 
Skyrmion (and anti-Skyrmion) . They argued that the 
SOI gives rise to a coupling between the uniform magne- 
tization M and the gauge field b by the term XM ■ b. In 
the ferromagnetic state, the spontaneous uniform magne- 
tization M leads to a finite and uniform 6, which acts as 
a uniform magnetic field. Lyanda-Geller et al. (Lyanda- 



Geller et al. 20011 also considered the AHE due to the 
spin chirality fluctuation in the incoherent limit where 
the hopping is treated perturbatively. This approach, 
applicable to the high-T limit, complements the theory 



of Ye et al. (Ye et al. 19991 



The Berry phase associated with non-coplanar spin 
conflgurations, the scalar spin chirality, was flrst con- 
sidered in theories of high-temperature superconductors 
in the context of the flux distribution generated by the 
complex order parameter of the resonating valence bond 
(RVB) correlation deflned by Xiji which acts as the trans- 
fer integral of the "spinon" between the sites i and j ( Lee 
et al. 2006). The complex transfer integral also appears 




t'.H (T) 

FIG. 15 (a) The colossal magnetoresistance p vs. H in 
Lai_xCaxMn03 (Tc = 265 K) at selected T. (b) The Hall re- 
sistivity Pxy VS. H at temperatures 100 to 360 K. Above Tc, 
Pxy is strongly influenced by the MR and the susceptibility X- 
[From Ref. |Matl et al.\ [19981] 




FIG. 16 Schematic view of spin chirality. When circulates 
among three spins to which it is exchange-coupled, it feels a 
fictitious magnetic field b with flux given by the half of the 
solid angle Q subtended by the three spins. [From Ref. Lee 
eran|2006,l 



in the double-exchange model specified by 

H = - ^t,j{cl^Cja + h.C.)-JHj2^- 



claO-a/dC^p, (2.2) 



where Jh is the ferromagnetic Hund's coupling between 
the spin cr of the conduction electrons and the localized 
spins Si. 

In the manganese oxides. Si represents the local- 
ized spin in t2g-orbitals, while and c are the opera- 
tors for Cg-electrons. (The mean-field approximations of 
Hubbard-like theories of magnetism, the localized spin 
may also be regarded as the molecular field created by 
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the conduction electrons themselves, in which case Jh is 
replaced by the on-site Coulomb interaction energy U .) 
In the limit of large Jh , the conduction electron spin s is 
forced to align with Si at each site. The matrix element 
for hopping from i —^ j is then given by 



et al. 1999 ) derived in the continuum approximation the 



— i ?■ 7 6 



COS 



(2.3) 



where \xi) is the two-component spinor spin wave func- 
tion with quantization axis \\Si. The phase factor e^'J 
acts like a Peierls phase and can be viewed as originat- 
ing from a fictitious magnetic field which influences the 
orbital motion of the conduction electrons. 

We next discuss how the Peierls phase leads to a gauge 
field, i.e. flux, in the presence of non-coplanar spin con- 
figurations. Let Si, Sj, and Sk be the local spins at sites 
i, j, and k, respectively. The product of the three trans- 
fer integrals corresponding to the loop i j ~* k ^ i 



{ni\nk){nk\nj){nj\ni) 

(1 + rii ■ Tij + Tij ■ Til; + rik ■ rij) + irii ■ {rij x n^) 

(2.4) 



where |ni) is the two-component spinor wavefunction of 
the spin state polarized along rii = Si/\Si\. Its imag- 
inary part is proportional to Si ■ {Sj x Sk), which cor- 
responds to the solid angle H. subtended by the three 
spins on the unit sphere, and is called the scalar spin 
chirality (Fig. 16). The phase acquired by the electron's 
wave function around the loop is e*^/^, which leads to 
the Aharonov-Bohm (AB) effect and, as a consequence, 
to a large Hall response. 

In the continuum approximation, this phase factor is 
given by the flux of the "effective" magnetic field b ■ dS = 
V X a • dS where dS is the elemental directed surface 
area defined by the three sites. The discussion implies 
that a large Hall current requires the unit vector n{x) — 
S{x) /\S{x)\ to fluctuate strongly as a function of x, 
the position coordinate in the sample. An insightful way 
to quantify this fluctuation is to regard n(a;) as a map 
from the x-y plane to the surface of the unit sphere (we 
take a 2D sample for simplicity). An important defect 
in a ferromagnet - the Skyrmion (Sondhi et al. 19931 



occurs when n{x) points down at a point x' in a region 
A of the x-y plane, but gradually relaxes back to up at 
the boundary of A. The map of this spin texture wraps 
around the sphere once as x' roams over A. The number 
of Skyrmions in the sample is given by the topological 
index 



Ns — dxdy n 
J A 



dn dn 

dx dy 



= [ dxdyh, (2.5) 

J A 



where the first integrand is the directed area of the image 
on the unit sphere. Ns counts the number of times the 
map covers the sphere as A extends over the sample. The 
gauge field b produces a Hall conductivity. Ye et al. (Ye 



coupling between the field b and the spontaneuos magne- 
tization M through the SOI. The SOI coupling produces 
an excess of thermally excited positive Skyrmions over 
negative ones. This imbalance leads to a net uniform 
"magnetic field" b (anti)parallel to M, and the AHE. 
In this scenario, p^y is predicted to attain a maximum 
slightly below Tc, before falling exponentially to zero as 

The Hall effect in the hopping regime has been dis- 
cussed by Holstein (Holstein 19611 in the context of 
impurity conduction in semiconductors. Since the en- 
ergies Sj and Ek of adjacent impurity sites may differ sig- 
nificantly, charge conduction must proceed by phonon- 
assisted hopping. To obtain a Hall effect, we consider 3 
non-coUinear sites (labelled as i = 1, 2 and 3). In a field 
H, the magnetic flux (j) piercing the area enclosed by the 
3 sites plays the key role in the Hall response. Accord- 
ing to Holstein, the Hall current arises from interference 
between the direct hopping path 1 — s- 2 and the path 
1 — > 3 — + 2 going via 3 as an intermediate step. Taking 
into account the changes in the phonon number in each 
process, we have 



{2,NxTl,NyTl), 
{2,NxTl,Nx,Tl), 
(2.6) 



where N\ , Ny are the phonon numbers for the modes 
A, A', respectively. 

In a field H, the hopping matrix element from Ri to 
Rj includes the Peierls phase factor exp[— i(e/c) dr ■ 
A{r)]. When we consider the interference between the 



two processes in Eq. (2.6 1, the Peierls phase factors 



combine to produce the phase shift exp(i27r0/(/)o) where 
(j)o — hc/e is the fiux quantum. By the Aharonov-Bohm 
effect, this leads to a Hall response. 

As discussed, this idea was generalized for the mangan- 
ites by replacing the Peierls phase factor with th e Berr y 
phase factor in Eq. (2.4 1 ( jLyanda-Geller et al. 20011. 
The calculated Hall conductivity is 

an = G({£}) cos(%/2) cos(0,fc/2) cos(0h/2) sin(r!/2), 

(2.7) 

where {e} is the set of the energy levels Sa (a = k), 
and Oij is the angle between rii and Uj. When the aver- 
age of an over all directions of Ua is taken, it vanishes 
even for finite spontaneous magnetization m. To obtain 
a finite a^y, it is necessary to incorporate SOI. 

Assuming the form of hopping integral with the SOI 
given by 



(2., 



the Hall conductivity is proportional to the average of 



[fifjfc • {nj X nk)\[ni ■ n2 x na]. 



(2.9) 
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FIG. 17 Comparison between experiment and the theoreti- 
cal prediction Eq.(2.10l. Scaling behavior between the HaU 

The soUd 



resistivity pu and the magnetization M is shown, 
line is a fit to Eq.| |2.10 1; the dashed line is the numerator of 
Eq.(2.10l only. There are no fitting parameters except the 
overall scale. [From Ref. Chun et al. 2000 ] 
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Taking the average of n's with m = M /Msat where Msat 
is the saturated magnetization, we finally obtain 



Pxy 



mix 

rxy 



(1 + m?Y 



(2.10) 



This prediction has been tested by the experiment of 
Chun et al. 



(Chun et al. 



scaling law for the anoma 



2000) shown in Fig. 17 



The 

ious Pxy as a function of |M| 
obtained near Tc is in good agreement with the experi- 
ment . 

Similar ideas have been used by Burkov and Ba- 



lents (Burkov and Balents 20031 to analyze the vari 



able range hopping region in (Ga,Mn)As. The spin- 
chirality mechanism for the AHE has also been applied 
to Cr02 ( Yanagihara and Salamon]|2002[ 2007), and the 
element Gd (Baily and Salamon 2005). In the former 



case, the comparison between p^y and the specific heat 
supports the claim that the critical properties of pxy are 
governed by the Skyrmion density. 

The theories described above assume large Hund cou- 
pling. In the weak- Hund coupling limit, a perturba- 
tive treatment in Jh has been developed to relate the 
AHE conductivity to the scalar spin chirality (Tatara 



and Kawamura 2002 1 . This theory has been applied to 



metallic spin-glass systems (Kawamura 20071. 



FIG. 18 Temperature dependence of the magnetization M 
(a). Hall resistivity p^y (b), Hall conductivity a^y (c), in four 
crystals of Lai-xSrxCoOs (0.17 < x < 0.30) (all measured 
in a field = 1 T). In Panels a, b and c, the data for x 
— 0.17 and 0.20 were multiplied by a factor of 200 and 50, 
respectively, (d) The Hall conductivity a^y at 1 T in the four 
crystals plotted against M. Results for x = 0.17 and 0.20 
were multiplied by factor s 50 and 5, respectively. [From Ref. 
Onose and Tokural 1200611 



at low T, and proposed the relevance of spin-ordered 
clusters and orbital disorder scattering to the AHE in 
the low-T limit. Subsequently, a more detailed inves- 
tigation of Lai-xSrxCoOs was rep orted by Onose and 
Tokura fOnose and Tokura 



20061. Figs. 18 



a. 



b and 
re- 



and axy, 
0.30. The 



c summarize the T dependence of M, p^y 
spectively in four crystals with 0.17 < x < 
variation of pxx vs. x suggests that a metal-insulator 
transition occurs between 0.17 and 0.19. Whereas the 
samples with x > 0.2 have a metallic Pxx-T profile, the 
sample with x — 0.17 is non-metallic (hopping conduc- 
tion). Moreover, it displays a very large MR at low T 
and large hysteresis in curves of M vs. H, features that 
are consistent with a ferromagnetic cluster-glass state. 



3. Lanthanum cobaltite 

The subtleties and complications involved in analyz- 
ing the Hall conductivity of tunable ferromagnetic ox- 
ides are well illustrated by the cobaltites. Samoilov et 



al. (Samoilov et al. 



and Salamon 2003 



1998 1, and Baily and Salamon (Baily 
I investigated the AHE in Ca-doped 
lanthanum cobaltite Lai-xCaxCoOs, which displays a 
number of unusual magnetic and transport properties. 
They found an unusually large AHE near Tc as well as 



When the Hall conductivity is plotted vs. M (Fig. 18 
d), (Txy shows a linear dependence on M for the most 
metallic sample {x = 0.30). However, for x — 0.17 and 
0.20, there is a pronounced downturn suggestive of the 
appearance of a different Hall term that is electron-like 
in sign. This is most apparent in the trend of the curves 
b 



of Pxy VS. T in Fig. 18 b. Onose and Tokura (Onose and 
Tokura 2006 1 propose that the negative term may arise 
from hopping of carriers between local moments which 
define a chirality that is finite, as discussed above. 
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4. Spin chirality mechanism in pyroclilore ferromagnets 

In the examples discussed in the previous subsection, 
the spin-chirahty mechanism leads to a large AHE at 
finite temperatures. An interesting question is whether 
or not there exist ferromagnets in which the spin chirality 
is finite in the ground state. 

Ohgushi et al. (Ohgushi et al. 2000) considered the 



ground state of the non-coplanar spin configuration in the 
Kagome lattice, which may be obtained as a projection 
of the pyrochlore lattice onto the plane normal to (1,1,1) 
axis. Considering the double exchange model Eq. (2.2 1, 



they obtained the band structure of the conduction elec- 
trons and the Berry phase distributi on. Quite similar 
to the Haldane model (Haldane 1988) or the model dis- 



cussed in Eq. (3.21 1, the Chern number of each band be- 
comes nonzero, and a quantized Hall effect results when 
the chemical potential is in the energy gap. 

Turning to real materials, the pyrochlore ferromag- 
net Nd2Mo207 (NMO) provides a test-bed for exploring 
these issues. Its lattice structure consists of two inter- 
penetrating sublattices comprised of tetrahedrons of Nd 
and Mo atoms, res pectively (the sublattices are shifted 
along the c-axis) ( Taguchi et al.\ 2001 Yoshii et al. 
|2000). While the exchange between spins on either sub- 
lattice is ferromagnetic, the exchange coupling J^/ be- 
tween spins of the conducting d-electrons of Mo and lo- 
calized /-electron spins on Nd is antiferromagnetic. 
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FIG. 19 Anomalous Hall effect in Nd2Mo207. Magnetic 
field dependence of (A) the magnetization and (B) the trans- 
verse resistivity {p^y) for different temperatures. [From Ref. 
[Taguchi et oE] |200T| ] 



The dependences of the anomalous Hall resistivity pxy 
on H at selected temperatures are shown in Fig. [19] The 
spins of Nd begin to align antiparallel to those of Mo be- 
low the crossover temperature T* = 40 K. Each Nd spin 
is subject to a strong easy-axis anisotropy along the line 
from a vertex of the Nd tetrahedron to its center. The 
resulting noncoplanar spin configuration induces a trans- 
verse component of the Mo spins. Spin chirality is ex- 



pected to be produced by the coupling J^j , which leads to 
the AHE of d-electrons. An analysis of the neutron scat- 
tering experiment has determined the magnetic structure 
(Taguchi et al. 2001). The tilt angle of the Nd spins is 



close to that expected from the strong limit of the spin 
anisotropy, and the exchange coupling Jfd is estimated 
as Jfd ~ 5 K. This leads to a tilt angle of the Mo spins of 
~5°. From these estimates, a calculation of the anoma- 
lous Hall conductivity in a tight-binding Hamiltonian of 
triply degenerate t2g bands leads to an ~ 20 (ficm)"^, 
consistent with the value measured at low T. In a strong 
H, this tilt angle is expected to be reduced along with 
Pxy. This is in agreement with the traces displayed in 
Fig. [191 

The T dependence of the Hall conductivity has also 
been analyzed in the spin-chirality scenario by incorpo- 
rating spin fluctuations ( Onoda and Nagaosa 2003 1 . The 



result is that frustration of the Ising Nd spins leads to 
large fluctuations, which accounts for the large p ob- 
served. The recent observation of a sign-c hange in a^y in 



a field H applied in the [1,1,1] direction (Taguchi et al 



2003 1 is consistent with the sign-change of the spin chi- 



rality. 

In the system GdsMosOr (GMO), in which Gd^+ {<f) 
has no spin anisotropy, the low-T AHE is an order- 
of-magnitude smaller than that in NMO. This is con- 



sistent with the spin-chirality scenario (Taguchi et al. 



20041. The effect of the spin chirality mechanism on 



the finite-frequency conductivity (jf{{uj) has been inves- 



tigated (Kezsmarki et al. 2005) 



In another work, Yasui et al. ( Yasui et al.\ 2006 



2007) performed neutron scattering experiments over a 
large region in the (i/, T)-plane with H along the [0,1,1] 
and [0,0,1] directions. By fitting the magnetization 
M^d{H, T) of Nd, the magnetic specific heat Cniag{H, T), 
and the magnetic scattering intensity ImngiQ, H, T), they 
estimated Jdj — 0.5 K, which was considerably smaller 
than estimated previously (Taguchi et al.\ 2001). Fur- 



thermore, they calculated the thermal average of the spin 
chirality {Si ■ Sj x Sk) and compared its value with that 
inferred from the AHE resistivity pxy . They have empha- 
sized that, when a 3-Tesla field is applied in the [0,0, 1] 
direction (along this direction H cancels the exchange 
field from the Mo spins), no appreciable reduction of p^y 
is observed. These recent conclusions have cast doubt on 
the spin-chirality scenario for NMO. 

A further puzzling feature is that, with H applied in 
the [1, 1, 1] direction, one expects a discontinuous transi- 
tion from the two-in, two-out structure (i.e. 2 of the Nd 
spins point towards the tetrahedron center while 2 point 
away) to the three-in, one-out structure for the Nd spins. 
However, no Hall features that might be identified with 
this cancellation have been observed down to very low T. 
This seems to suggest that quantum fiuctuations of the 
Nd spins may play an important role, despite the large 
spin quantum number {S — |). 

Machida has discussed the possible relevance of 
spin chirality to the AHE in the pyrochlore Pr2lr207 
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FIG. 20 Plot of AHE conductivity cjahe vs. conduc- 
tivity o" for anatase Tii_xCox02-5 (triangles) and rutile 
Tii_xC0x02-5 (diamonds). Grey symbols are data taken by 
other groups. The inset shows the expanded view of data for 
anatase with x= 0.05 (the open and closed triangles a re for 
T > 150 K and T < 100 K, respectively. [From Ref. [Ueno 
et al.\ 120071 



(Machida et al. 2007a|b I. In this system, a novel "Kondo 
effect" is observed even though the Pr'^+ ions with S = 1 
are subject to a large magnetic anisotropy. The magnetic 
and transport properties of i?2Mo207 near the phase 
boundary between the spin glass Mott insulator and fer- 
romagnetic metal by changing the rare earth ion R has 
been studied ( jKatsufuji et al. 



20001 



5. Anatase and Rutile Tii_xCox02- 



In thin-film samples of the ferromagnetic semiconduc 

Ueno et al. ( |Ueno et al. 



tor anatase Tii_xCOxO 



2-5, 



2008|) have reported scahng between the AHE resistance 

The AHE conductivity 
scales with the conductivity axx as a^,^ cx a^-^ 



and the magnetization M. 



XV -XX (Fig.|20l. 
A similar scaling relation was observed in another poly- 
morph rutile. See also Ref. |Ramaneti et al.\ |2007| for 
related work on Co-doped Ti02. 



C. Ferromagnetic semiconductors 

Ferromagnetic semiconductors combine semiconduc- 
tor tunability and collective ferromagnetic properties in 
a single material. The most widely studied ferromag- 
netic semiconductors are diluted magnetic semiconduc- 
tors (DMS) created by doping a host semiconductor with 
a transition metal which provides a localized large mo- 
ment (formed by the d-electrons) and by introducing 
carriers which can mediate a ferromagnetic coupling be- 
tween these local moments. The most extensively studied 
are the Mn based (III,Mn)V DMSs, in which substitut- 
ing Mn for the cations in a (III,V) semiconductor can 
dope the system with hole carriers; (Ga,Mn)As becomes 



ferromagnetic beyond a concentration of 1%. 

The simplicity of this basic but generally correct model 
hides within it a cornucopia of physical and materials 
science effects present in these materials. Among the 
phenomena which have been studied are metal-insulator 
transitions, carrier mediated ferromagnetism, disorder 
physics, magneto-resistance effects, magneto-optical ef- 
fects, coupled magnetization dynamics, post-growth de- 
pendent properties, etc. A more in-depth discussion of 
these materials, both from the experimental and theoret- 
ical point of view, can be found in the recent review by 
Jungwirth et al (Jungwirth et al. 2006). 



The AHE has been one of the most fundamental char- 
acterization tools in DMSs, allowing, for example, direct 
electrical measurement of transition temperatures. The 
reliability of electrical measurement of magnetic prop- 
erties in these materials has been verified by compari- 
son with remnant magnetization measurements using a 
SQUID magnetometer (Ohno et al. 19921. The relative 



simplicity of the effective band structure of the carriers 
in metallic DMSs, has made them a playing ground to 
understand AHE of ferromagnetic systems with strong 
spin-orbit coupling. 

Experimentally, it has been established that the AHE 
in the archetypical DMS system (Ga,Mn)As is in the 



mechanism 


i.e. p^^ cx pI^ 


Chun et al. 


2007 


Edmonds 


et al. 


2002 


Pu et al. 


2008 


kuzmetov et al. V 


i004 


1. I'he 


studies of . 


Edmonds et al. 


2002 and Chun et al. 


12007 



have established this relationship in the non-insulating 

materials by extrapolating the low temperature Pxy{B) 
to zero field and zero temperature. This is illustrated in 



Fig. 21 where metallic samples, which span a larger range 
than the ones studied by [Edmonds et al.\ |2002[ show a 
clear Rs ~ p^^ dependence. 

DMSs grown requires non-equilibrium (low tempera- 
ture) conditions and the as-grown (often insulating) ma- 
terials and post-grown annealed metallic materials show 
typically different behavior in the AHE response. A 
similar extrapolating procedure performed on insulating 
(Ga,Mn)As seems to exhibit a somewhat linear depen- 
dence of Rs on pxx- On the other hand, considerable 
uncertainty is introduced by the extrapolation to low 
temperatures because p^^ diverges and the complicated 
magetoresistance of pxx is a priori not known in the low 
T range. 



A more recent study by Pu et al. (Pu et al. 2008 lof 
(Ga,Mn)As grown on InAs, such that the tensile strain 
creates a perpendicular anisotropic ferromagnet, has es- 
tablished the dominance of the intrinsic mechanism in 
metallic (Ga,Mn)As samples beyond any doubt. Measur- 
ing the longitudinal thermo-electric transport coefficients 
( Pxx, Pxy, axy, and Uxx where J = aE + ai-VT)), one 

can show that given the Mott relation a — ^ 3^'^ {§%) Ef 
and the empirical relation Pxy{B = 0) = XM^p^^, the re- 
lation between the four separately measured transport 
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FIG. 21 (a) Gai_i:Mni:As samples that show insulating and 
metallic behavior defined by dpxx/dT near T — Q. (b) Rs 
vs. pxx extrapolated from pxy{B) data to zero field and low 
temperatures. [From Ref. [Chun et al.\ |2007| ] 
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(2.11) 



The fit to the ^ and n parameters are shown in Fig 
Fixing n = 1 does not produce any good fit to the data 
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as indicated by the dashed lines in Fig. 22 These data 
excludes the possibility of a n = 1 type contribution to 
the AHE in metallic (Ga,Mn)As and further verify the 
validity the Mott relation in these materials. 

Having established that the main contribution to the 
AHE in metallic (Ga,Mn)As is scattering-independent 
contributions, rather than skew-scattering contributions, 
DMSs are an ideal system to test our understanding of 
AHE. In the regime where the largest ferromagnetic crit- 
ical temperatures are achieved ( for doping levels above 
1.5% ), semi-phenomenological models that are built on 
Bloch states for the band quasiparticles, rather than lo- 
calized basis states appropriate for the localized regime 



(Berciu and Bhatt 



2001), provide the natural starting 



amiltonian which reproduces many 



point for a model 

of th e observed experimental effect s ([Jungwirth et al. 



2006 Sinova and Jungwirth 20051. Recognizing that 



the length scales associated with holes in the DMS com- 
pounds are still long enough, a. k ■ p envelope function 
description of the semiconductor valence bands is appro- 
priate. To understand the AHE and magnetic anisotropy, 
it is necessary to incorporate intrinsic spin-orbit coupling 
in a realistic way. 

A successful model for (Ga,Mn)As is specified by the 
effective Hamiltonian 

n = nKL + Jpd J2 ■ ^(^)'^(^ - + (2-12) 
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FIG. 22 (Top) Zero B field p^y and pxx for four samples 
grown on InAs substrates (i.e. perpendicular to plane easy 
axis). Annealed samples, which produce perpendicular to 
plane easy axes, are marked by a *. The inset indicates the 
B dependence of the 7% sample at 10 K. (Bottom) Zero-field 
Nerst coefficient axy for the four samp les. The solid red curves 
indicate the best fit using Eq. ( |2.11 1 and the dashed curves 
the best fit setting n=l. [From Ref. 



Pu et al. 2008 ] 



where Hkl is the six-band Kohn-Luttinger (KL) k ■ p 
Hamiltonian (Dietl et al. 20011, the second term is the 



short-range antiferromagnetic kinetic-exchange interac- 
tion between local spin Sj at site Rj and the itinerant 
hole spin (a finite range can be incorporated in more real- 
istic models), and Tidis is the scalar scattering potential 
representing the difference between a valence band elec- 
tron on a host site and a valence band electron on a Mn 
site and the screened Coulomb interaction of the itinerant 
electrons with the ionized impurities. 

Several approximations can be used to vastly simplify 
the above model, namely, the virtual crystal approxima- 
tion (replacing the spatial dependence of the local Mn 
moments by a constant average) and mean field theory 
in which quantum and thermal fiuctuations of the local 



moment spin-orientations are ignored (Dietl et al. 2001 



Jungwirth et al. 20061. In the metallic regime, disorder 
can be treated by a Born approximation or by more so- 
phisticated, exact-diagonalization or Monte-Carlo meth- 
2002a^ 'Sc hliemann and MacDon- 
2002; Yang et al.\\2003^. 



ods ([Jungwirth et al. 
Hdl 



2002 



Sinova et al] 



Given the above simple model Hamiltonian the AHE 
m be computed if one assumes that the intrinsic 
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Berry's phase (or Karplus-Luttinger) contribution will 
most likely be dominan, because of the large SOI of the 
carriers and the experimental evidence showing the domi- 
nance of scattering-independent mechanisms. For practi- 
cal calculations it is useful to use, as in the intrinsic AHE 



studies i n the oxides (Fang et at. 2003 Mathieu et al.\ 
2004a|b I, the Kubo formalism given in Eq. ( 1.2 ) with dis- 
order induced broadening, F, of the band-structure (but 
no side-jump contribution). The broadening is achieved 
by substituting one of the (£„(fc) — £„'(fc)) factors in the 
denominator by (£„(fc) — £n'{k) + iT). Applying this 
theory to metallic (III,Mn)V materials using both the 4- 
band and 6-band k ■ p description of the valence band 
electronic structure one obtains results in quantitative 
agreement with experimental data in (Ga,Mn)As and 
(In,Mn)As DMS ( |Jungwirth ei! aL| |2002bD. In a foll ow 
up calculation Jungwirth et al. ( Jungwirth et aL[[2003 l, a 
more quantitative comparison of the theory with experi- 
ments was made in order to account for finite quasiparti- 
cle lifetime effects in these strongly disorder systems. The 
effective lifetime for transitions between bands n and n' , 
Tn^n' = l/rn,n'i Can be Calculated by averaging quasipar- 
ticle scattering rates obtained from Fermi's golden rule 
including both screened Coulomb and exchange poten- 
tials of randomly distributed substitutional Mn and com- 
pensating defects as done in the dc Boltzman transport 
studies (Jungwirth et al. 2002a Sinova et al. 20021. A 



systematic comparison between theoretical and experi- 
mental AHE data is shown in Fig. 23 (Jungwirth et al. 
2003 ). The results are plotted vs. nominal Mn concentra- 



tion X while other parameters of the seven samples stud- 
ied are listed in the figure legend. The measured aAH 
values are indicated by filled squares; triangles are theo- 
retical results obtained for a disordered system assuming 
Mn-interstitial compensationg defects. The valence band 
hole eigenenergies e„fe and eigenvectors \nk) are obtained 
by solving the six-band Kohn-Luttinger Hamiltonian in 
the presence of the exchange field, h ~ NMnSJpdZ (Jung- 



wirth et al. 



2006). Here N, 



Mn 



Ax /a 



|)jv/s is the Mn" 



density in the Mna;Gai_2;As epilayer with a lattice con- 
stant a^MSi the local Mn spin S — 5/2, and the exchange 
coupling constant Jpd = 55 mcV nm^"^. 

In general, when disorder is accounted for, the the- 
ory is in a good agreement with experimental data over 
the full range of Mn densities studied from x ~ 1.5% to 
X = 8%. The effect of disorder, especially when assum- 
ing Mn-interstitial compensation, is particularly strong 
in the x = 8% sample shifting the theoretical aAH much 
closer to experiment, compared to the clean limit the- 
ory. The remaining quantitative discrepancies between 
theory and experiment have been attributed to the reso- 
lution in measuring experimental hole and Mn densities 



( Jungwirth et al. 2003 1 



We conclude this section by mentioning the anoma- 
lous Hall effect in the non-metallic or insulating/hopping 
regimes. Experimental studies of (Ga,Mn)As digital fer- 
romagnetic heterostructures, which consist of submono- 
layers of MnAs separated by spacer layers of GaAs, 
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FIG. 23 Comparison between experimental and theoretical 
anomalous Hall conductivities. After Ref. [Jungwirth et aL\ 

mm 



have shown longitudinal and Hall resistances of the hop- 
ping conduction type, R^x oc T" exp[(To/r)''], and have 
shown that the anomalous Hall resistivity is dominated 
by h opping with a sub l inear dependence of Rah on 
Rxx ( [Allen et al] [2004| [Shen et al] [2008[ ), similar to 

Studies in 



experimental observations on other materials, 
this regime are still not as systematic as their metallic 
counterparts which clearly indicate a scaling power of 2. 
Experiments find a sublinear dependence of Rah ^ Rxx' 
with P ~ 0.2 — 1.0 depending on the sample studies (Shen] 



et al. 2008 1. A theoretical understanding for this hop- 
ping regime remains to be worked out still. In previous 
theoretical calculations based on the hopping conduction 
with the Berry phase (Burkov and Balents 2003 1 showed 
the insulating behavior Rah ^ oo as Rxx oo but 
failed to explain the scaling dependence of Rah on Rxx ■ 



D. Other classes of materials 



1. Spinel CuCr2Se4 

As mentioned in earlier sections, a key prediction of 
the KL theory (and its modern generalization, based on 
the Berry-phase approach) is that the AHE conductivity 
is independent of the carrier lifetime r (dissipation- 
less Hall current) in materials with moderate conductiv- 
ity. This implies that the anomalous Hall resistivity 
varies as p^. Moreover, tend to be proportional to 



the observed magnetization M^. We write 



AH 



ShM, 



(2.13) 



with Sh a constant. 

Previously, most tests were performed by comparing 
Pyx vs. p measured on the same sample over an ex- 
tended temperature range. However, because the skew- 
scattering model can also lead to the same prediction 
Pyx ~ when inelastic scattering predominates, tests 
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FIG. 24 Log-log plot of the quantity \p'xy\/nh vs. p in twelve 
crystals of Br-doped spinel CuCr2Se4_xBrx with Uh the hole 
density (p is measured at 5 K; pxy is measured at 2 and 5 K). 
Samples in which pyx is electron-like (hole-like) are shown as 
open (closed) circles. The straight-line fit implies Ipxyl/nn = 
Ap" with a = 1.95 ± 0.08 and A = 2.24x10"^^ (SI units). 



[From Ref. Lee et al. 2004 



at finite T are inconclusive. The proper test requires a 
system in which p at 4 K can be varied over a very large 
range without degrading the exchange energy and mag- 
netization M. 

In the spinel CuCr2Se4, the ferromagnetic state is sta- 
bilized by 90° superexchange between the local moments 
of adjacent Cr ions. The charge carriers play only a weak 
role in the superexchange energy. The experimental proof 
of this is that when the carrier density n is varied by a 
factor of 20 (by substituting Se by Br), the Curie tem- 
perature Tc decreases by only 100 K from 380 K. Signif- 
icantly, M at 4 K changes negligibly. The resistivity p at 
4 K may be varied by a factor of 10'^ without weakening 
M. Detailed Hall and resistivity measurements were car- 



ried out by Lee et al. (Lee et al. 
of CuCr2Se4, 



2004 ) on twelve crystals 
measured at 5 



Brx. They found that Py. 
K changes sign (negative to positive) when x exceeds 0.4. 
At a; = 1, Pyx attains very large values (~ 700 /ificm at 

5K). 

Lee et al. (Lee et al. 20041 showed that the magni- 
tude Ipyx 



FIG. 25 Combined plots of a^y in the Heusler alloy Co2CrAl 
versus the magnetization M at selected temperatures from 36 
to 278 K. The inset shows the inferred values of a^y ^ = Oxy/M 
at each T. [From Ref. [Husmann and Singh[ |2006 ] 



eral groups ( [Block et 11] |2004[ [Galanakis et al.\ |2002 1 



have argued that the conduction electrons are fully spin 
polarized ( "half metal" ) . The absence of minority carrier 
spins is expected to simplify the analysis of the Hall con- 
ductivity. Hence this system is potentially an important 
system to test theories of the AHE. 

The AHE has been investigated on single crystals with 
stoichiometry Co2.o6Cri.o4Alo gg (Husmann and Singh 
j2006|. Below the Curie temperature Tc = 333 K, the 
magnetization M increases rapidly, eventually saturating 
to a low-T value that corresponds to 1.65 /is (Bohr mag- 
neton) per formula unit. Husmann and Singh show that, 
below -310 K, M{T) fits well to the form [1 - [T/Tcf] 5 . 
Assuming that the ordinary coefficient i?o is negligible, 
they found that the Hall conductivity a^y — Pyx/p^ is 
strictly linear in M (expressed as axy = (^hM) over the 
T interval 36-278 K (Fig. |25]). They interpret the hn- 
ear variation as consistent with the intrinsic AHE theory. 
The value of aj^ = 0.383 G/(47rilcm) inferred is similar 
to values derived from measurements on the dilute Ni 
alloys, half Heuslers and silicides. 



n varies as over 3 decades in p (Fig. 24 1, 3. Fei_yCoj,Si 



consistent with the prediction of the KL theory. 

The AHE in the related materials CuCr2S4, 
CuxZnxCr2Se4 {x = ^) and Cu3Te4 has been investi- 



gated by Oda et al. (Oda et al. 20011 



2. Heusler Alloy 

The full Heusler alloy Co2CrAl has the L2i lattice 
structure and orders ferromagnetically below 333 K. Sev- 



The silicide FeSi is a non-magnetic Kondo insulator. 
Doping with Co leads to a metallic state with a low den- 
sity p of holes. Over a range of Co doping (0.05 < y < 
0.8) the ground state is a helical magnetic state with 
a peak Curie temperature Tc — 50 K. The magnetiza- 
tion corresponds to I pB (Bohr magneton) per Co ion. 
The tunability allows investigation of transport in a mag- 



netic system with low p. Manyala et al. (Manyala et al. 



20041 observe that the Hall resistivity pn increases to 
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FIG. 26 Hall conductivity a^y of ferromagnetic metals and 
heavy fermion materials (collected at 1 kG and 5 K unless 
otherwise noted). Small solid circles represent Fei_yCoySi 
for 5 < T < 75 K and 500 G< H < 50 kG with y = 0.1 
(filled circles), 0.15 (filled triangles), y = 0.2 (+), and 0.3 
(filled diamonds). MnSi data (5-35 K) are large, solid squares 
connected by black line. Small asterisks are (GaMn)As data 



for 5< T < 120 K. The rising line 



M is consistent 



with the intrinsic AHE in itinerant ferromagnets while the 
falling line axy ~ M" *^ applies to heavy fermions. [From Ref. 



Manyala et al. 2004 



^1.5 /zficm (at 5 K) at the doping y= 0.1. By plotting 
the observed Hall conductivity a^y against M , they find 
Oxy = ShM with Sh ~ 0.22 G/(47rricm). In contrast, in 
heavy fermion systems (which include FeSi) <Jxy ^ M^^ . 
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FIG. 27 (a) HaU resistivity p^^ vs. H in MnSi at selected T 
from 5 to 200 K. At high T, pyx is linear in H, but gradually 
acquires an anomalous component p'y^ — Pyx — RoB with a 
prominent "knee" feature below Tc = 30 K. (b) Matching of 
the field profiles of the anomalous Hall resistivity p'yx to the 
profiles of p^M, treating Sh and Ro as adjustable parameters. 
Not e the positive cur vature of the high-field segments. [From 
Ref. ILee et al.\ 1200711 



4. MnSi 

MnSi grows in the non-centrosymmetric B20 lattice 
structure which lacks inversion symmetry. Competi- 
tion between the exchange energy J and Dzyaloshinsky- 
Moriya term D leads to a helical magnetic state with a 
long pitch A (^--^180 A). At ambient pressure, the heli- 
cal state forms at the critical temperature Tc = 30 K. 
Under moderate hydrostatic pressure P, Tc decreases 
monotonically, reaching zero at the critical pressure Pc 
= 14 kbar. Although MnSi has been investigated for 
several decades, interest has been revived recently by a 
neutron scattering experiment which shows that, above 
Pc, MnSi displays an unusual magnetic phase in which 
the sharp magnetic Bragg spots at P < Pc are replaced 
by a Bragg sphere ( |Pfleiderer et al. 20041. Non- Fermi 



liquid exponents in the resistivity p vs. T are observed 
above Pc. 

Among ferromagnets, MnSi at 4 K has a low resistivity 
(p ~ 2-5 /zficm). The unusually long carrier mean-free- 
path £ implies that the ordinary term a^J^ ~ is greatly 



enhanced. In addition, the small p renders the total Hall 
voltage difficult to resolve. Both factors greatly compli- 
cate the task of separating o-^y^ from the AHE conductiv- 
ity. However, the long £ in MnSi presents an opportunity 
to explore the AHE in the high-purity limit of ferromag- 
nets. Using high-resolution measurements of the Hall 
resistivity py,j: (Fig 
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a), Lee et al. (Lee et al. 20071 



recently accomplished this separation by exploiting the 
large lon gitudi nal m agnetoresistance (MR) p{H). From 
Eq. ^ and ( |2.13[ ), we have Py^iH) = Shp{HYM{H), 

yx — Pyx = Pyx — RqB 



where p. 



At each T, the field 
profiles of p'y,-^{H) and M{H) are matched by adjusting 
the two 77-independent parameters Sh and Rq (Fig. 27 
b). The inferred parameters Sh and Rq are found to be 



T independent below Tc (Fig.[28|a). 

The Hall effect of MnSi under hydrostatic pressure (5- 
11.4 kbar) was measured recently (Lee et al. 20081. In 



addition to the AHE and OHE terms and a"y , Lee 
et al. observed a well-defined Hall term (t£, with an un- 

jy y 



usual profile. As shown in Fig. 29 (note that in the figure 
and a^y^ are labeled a^y and respectively) , the 
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FIG. 28 



Comparison of the anomalous Hall conductivity 
measured in a 1-Tesla field 



and the ordinary term a^y 

a^„, inferred from the measured M (solid curve 



'xy 

in MnSi. 

and Eq. (2.131, is strictly independent of £. Its T dependence 
reflects that of M(T). a^y ~ is calculated from Rq. The 



^xy 

conductivity at zero H, a 
[From Ref. |Lee et al.\ |2007 



is shown as a dashed curve. 



new term appears abruptly at 0.1 T, rapidly rises to a 
large plateau value, and then vanishes at 0.45 T (curves 
at 5, 7 and 10 K). From the large magnitude of cr^j,, and 
its restriction to the field interval in which the cone angle 
is non-zero, the authors argue that it arises from the cou- 
pling of the carrier spin to the chiral spin textures in the 



helical magnetization, as discussed in Sec. II. B The au- 
thors note that MnSi under pressure provides a very rare 
example in which the three Hall conductivities co-exist 
at the same T. 



5. MnsGea 

The AHE in thin-film samples of MnsGea was investi- 
gated by Zeng et al. (Zeng et al. 20061. They express 



the AHE resistivity pah, which is strongly T depen- 
dent (Fig. 30 i) , as the sum of the skew-scattering term 



a{M)pxx and the intrinsic term b(M)p'^^, viz. 

PAH = a{M)pxx + b{M)pl^. (2.14) 
the intrinsic AHE conductivity 



The quantity 6(M) is 

O'AH-int- 

To separate the 2 terms, Zeng et al. plotted the quan- 
tity PAH /{M{T)pxx) against pxx with T as a parameter. 
For T < 0.8T(7, the plot falls on a straight line with 
a small negative intercept (solid squares in Fig. 
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The intercept yields the skew-scattering term a{ATj/M 
whereas the slope gives the intrinsic term 6(M) /Al . From 
the constant slope, they derive their main conclusion that 
aAH-int varies linearly with M. 

To account for the linear-M dependence, Zeng et 



al. ( Zeng et al.\ 2006 1 identify the role of long- wavelength 
spin waves which cause fiuctuations in the local direction 
of M{x) (and hence of fJ). They calculate the reduc- 
tion in O'AH-int Jind show that it varies linearly with M 
(Fig.lSOlc). 
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FIG. 29 (Main panel) ~Pyx vs. H in MnSi under hydrostatic 
pressure P = 11.4 kbar at several T < Tc, with H nominally 
along (111). The large Hall anomaly observed (electron-like 
in sign) arises from a new chiral contribution a^y to the total 
Hall conductivity. In the phase diagram (inset) the shaded 
region is where axy is resolved. The non-Fermi liquid region 
is shaded blue. [From Ref. |Lee et al.\ [2008] ] 
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FIG. 30 (a) The T dependence of pah in Mns Gea be- 
fore (solid squares) and after (solid circles) subtraction of the 
skew-scattering contribution. The thin curve is the resistivity 
p. (b) Plots of Pah / (M pxx) vs. pxx before (solid squares) 
and after (solid circles) subtraction of skew-scattering term, 
(c) Comparison of calculated aAH-int vs. (open circles) 
with experimental values before (squares) and after (solid ci 
cles) subtraction of skew-scattering term. [From Ref. Zeng 
[eFat] [2Q06| ] 
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6. Layered dichalcogenides 

The layered transition-metal dichalcogenides are com- 
prised of layers weakly botmd by the van der Waals force. 



Parkin and Friend (Parkin and Friend 1980 1 have shown 



that a large number of interesting magnetic systems may 
be synthesized by intercalating 3d magnetic ions between 
the layers. 

The dichalcogenide FexTaS2 typically displays proper- 
ties suggestive of a ferromagnetic cluster-glass state at 
low T for a range of Fe content x. However, at the com- 
position X — ^, the magnetic state is homogeneous. In 
single crystals of Fe i TaS2, the easy axis of M is parallel 
to c (normal to t he TaS2 layers). Morosan et al. (Ma 
rosan et al.\ 2007[ ) observed that the curves of M vs. H 
display very sharp switching at the coercive field at all 
T <Tc (160 K). In this system, the large ordinary term 



NH 

xy 



complicates the extraction of the AHE term 



Converting the py x-H curves to axy-H cu rves (Fig. 31 1 
Checkelsky et al. ( Checkelsky et al. 2008 ) infer that the 

2a^" by assuming that 
This method provides a direct mea- 

As shown 



jum p mag nitude ^Oxy equals ^^^y 
Eq. (plSl) is vahd ^' 

AH 
>xy 

in Fig. |3l|b, both the inferred cTx^ and measured M 
are nearly T-independent below 50 K, but the former 
deviates sharply downwards above 50 K. Checkelsky et 
al. (Checkelsky et al. 



surcmcnt of crt^,, without knowledge of i?o 




50 100 150 
T(K) 



FIG. 31 (a) Hystersis loops of <Jxy vs. H in FeiTaS2 calcu- 

4 

20081) propose that the deviation lated from the measured p and Py^ curves. The hnear portions 



represents a large, negative inelastic-scattering contribu- correspond t^o while the jump magnitude Act,,, equals 
tion cr™ that involves scattering from chiral textures of 



2a. 



the spins which increase rapidly as T approaches . 
In support, they show that the curve of a^J^/M{T) vs. 
T matches (within the resolution) that of Ap{T)^ with 
Ap(T) = p(T)-p(0). 



. (b) Comparison of Aa^y with the magnetization M 
measured at 0.1 T. Within the resolution, Aaxy seems to be 
proportional to M below 50 K, but deviates sharply from M 
at higher temperatures, reflecting the growing dominanc e of 
a negative, inelastic-scattering term a'J^. [From Ref. | Check-] 
elsky et al.\ 120081] 



theoretically the model 



E. Localization and AHE 

The role of localization in the anomalous Hall effect is 
an important issue, and there have been several works 
on this subject. ( For a review of theoretical works, see 
Ref. [Woelfle and Muttalib[[2006|) The weak localization 



effect on the normal Hall effect due to the external mag- 



netic field has been studied by Fukuyama (Fukuyama 
[l980( ), and the relation Sa^y^/a^y = 25afx'^/uxx has 
been obtained where 50^^ represents the correction of 
the physical quantity O by the weak localization effect. 
This means that the Hall coefficient is not subject to the 
change due to the weak localization, i.e. Sp^^if^ = 0. An 



early experiment by Bergmann and Ye (Bergmann and 



Ye 1991 1 on the anomalous Hall effect in the ferromag- 



netic amorphous metals showed almost no temperature 
dependence of the anomalous Hall conductivity, while 
the diagonal conductivity shows the logarith mic temper- 
ature dependence . Langenfeld and Woel fle ( Langenfeld 
and Wolflel 119911 iWoelfle and Muttalibl 120061) studied 



+ E E e*^'""'")-^-[F + z(fcxfc')- J,]4',.,Cfc,. 

(2.15) 

where Ji is proportional to the angular momentum of the 
impurity at Ri , and the term containing it describes the 
spin-orbit scattering. They discussed the logarithmic cor- 
rection to the anomalous Hall conductivity in this model, 
and found that Coulomb anomaly terms vanished iden- 
tically, and the weak localization correction was cut-off 
by the phase-breaking lifetime due to the skew scat- 
tering, explaining the experiment by Bergmann and Ye 
(Bergmann and Ye 1991) . Dugaev, Crepieux, and Bruno 
( Dugaev et al. 2002 1 found that for the two-dimensional 



case the correction to the Hall conductivity was loga- 
rithmic in the ratio max(Ttr/Tsoj ^tr/T^), with Ttr being 
the transport lifetime and tso the lifetime due to the 
spin-orbit scattering. More recently, Mitra et al. (Mi 
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FIG. 32 The resistance (7?o)-temperature dependence of the 
coefficients A_R and AyiH defined in Eq. (2.161. Different sym- 
bols correspond to different methods of sample preparation. 
[From Ref. |Mitra et al'] |2007[] 



tra et al. 20071 first observed the logarithmic tempera- 



ture dependence of the anomalous Hall conductivity in 
the ultrathin film of polycrystalline Fe of sheet resis- 
tance Rxx less than ^ 3 kfi. They defined the quantities 
A^((5) = {l/RoLoo){SQ/Q) for the physical quantity Q 
with respect to the reference temperature T = Tg = 5 K 



by 6Q = QiT)-Q{To), Rq = R^To), and Lqo 
Defining the coefficients Ar and Ah by 



Fig. |32] shows the i?o-dependence o f the coefficients Ar 
and Aah defined in the above Eq. (2.16). 



The change in the interpretation comes from the fact 
that the phase-breaking lifetime in the Fe film is 
mostly from scattering by the magnons and not from 
skew scattering, which allows a temperature regime 
where max(l/rs, l/rso, ^//) << T0 << l/^r ( r^: spin- 
flip scattering time, ujh- internal magnetic field in the fer- 
romagnet ). In this region, they found that the weak lo- 
calization effect can lead to coefficient of the logarithmic 
temperature dependence of cr^y proportional to the factor 
t^fy^'Vl^fy + where af^^*'^ is the contri- 

bution from the skew scattering (side jump) mechanism. 
Assuming that that ratio cr^y^ /a^y^ decreases as the 
sheet resistance Rq increases, they were able to explain 
the sample-dependence of the logarithmic correction to 



the anomalous Hall conductivity (Mitra et al. 20071 



The absence of the logarithmic term in ( Bergmann and 



Ye 1991 1 is interpreted as being due to a large ratio of 



^xy"^ l^xy i'^ their sample. It is interesting that the 
ratio cr^y^ /a^y'^ can be estimated from the coefficient 
of the logarithmic term. 

Up to now, we have discussed the weak localization 



effect. When the disorder strength increases, a metal- 
insulator transition will occur. For a normal metal under 
external magnetic field, the system belongs to the uni- 
tary class, and in 2D all the states are localized for any 
disorder (Lee and Ramakrishnan 19851. In the quan- 



tum Hall system, however, the extended states can sur- 
vive at discrete energies at the center of the broadened 
density of states at each Laudau level. This extended 
state carries the quantum Hall current. Field theoret- 
ical formulation of this localization problem has been 
developed (Prange and Girvin 1987). In the presence 



of the external magnetic field, a Chern-Simons term ap- 
pears in the non-linear sigma model whose coefficient is 
axy. Therefore, the scaling variables are axx and (Jxy, 
i.e., the two-parameter scaling theory should be applied 
instead of the single-parameter scaling. It has been dis- 
cussed that the scaling trajectory in the (Jxy — <^xx plane 
has the fixed point at (tTxy, fxa;) = (("■+ l/2)(e^//i),(To) 
where (Tq is some finite value, and axy scales to the 
quantized value n(e'^/h) when the initial value (given 
by the Boltzmann transport theory) lies in the range 
(n- l/2)(e2//i) < ai°J < {n+l/2){e'^/h). In contrast to 
this quantum Hall system, there is no external magnetic 
field or the Landau level formation in the anomalous Hall 
system, and it is not trivial that the same two-parameter 
scahng theory applies to this case. 



Onoda and Nagaosa ( Onoda and Nagaosal 2003 1 stud- 



/tt/i. i ed this problem using the generalized Haldane model 



(Haldane 1988) which shows the quantum Hall effect 



without the external magnetic field. They calculated the 
scaling function of the localization length in the finite- 
width stripe sample in terms of the iterative transfer 



matrix method by MacKinnon (MacKinnon and Kramer 



(2.16) 119831), and found that two-parameter scaling holds even 



without an external magnetic field or Landau level for- 
mation. For the experimental realization of this quan- 
tized anomalous Hall effect, the \a^J\ given by Boltz- 
mann transport theory (without the quantum correction 
) is larger than e^/(2/i), and the temperature is lower 



than Ti^oc ^ Spe 



where £p is the Fermi energy 



and c is a constant of the order of unity. Therefore, the 

Hall angle a^y ja^xx should not be so small, hopefully of 
the order of 0.1. However, in the usual case, the Hall an- 
gle is at most 0.01, which makes the quantized anomalous 
Hall effect rather difficult to realize. 



III. THEORETICAL ASPECTS OF THE AHE AND 
EARLY THEORIES 

In this section, we review recent theoretical develop- 
ments as well as the early theoretical studies of the AHE. 
First, we give a pedagogical discussion on the difference 
between the normal Hall effect due to the Lorentz force 



and the AHE (Sec. III.Al. In Sec. III.B we discuss the 
topological nature of the AHE. In this Sec. |III.C| we 
present a wide survey of the early theories from a modern 
viewpoint in order to bring them into the context of the 
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present linear transport theory formalisms now used as 
a framework for AHE theories. 



A. Symmetry considerations and analogies between normal 
Hall effect and AHE 

Before describing these recent developments, we pro- 
vide an elementary discussion that may facilitate under- 
standing of the following sections. 

The Hall effect is one of the fundamental transport 
phenomena in solid-state physics. Its occurrence is a 
direct consequence of broken time-reversal symmetry in 
the ferromagnetic state, T . The charge current J is T- 
odd,i.e., it changes sign under time reversal. On the 
other hand, the electric field E is T-even. Therefore, 
Ohm's law 



J = (jE 



(3.1) 



relates two quantities with different T-symmetries, which 
implies that the conductivity a must be associated with 
diispative irreversible processes, and indeed we know that 
the Joule heating Q = aE"^ /2 always accompanies the 
conductivity in Eq. (3.1 1. This irreversibility appears 



only when we consider macroscopic systems with contin- 
uous energy spectra. 

Next, we consider the other aspect of the T-symmetry, 
i.e., the consequences of the T-symmetry of the Hamil- 
tonian which governs the microscopic dynamics of the 
system. This important issue has been formulated by 
Onsager, who showed that the response functions satisfy 
the following relation (Landau et al.\ 



19841 



:B)=e^e0Kp^{uj;,r',r;-B), (3.2) 



where Kafj^uj; , r, r'; B) is the response of a physical quan- 
tity a at position r to the stimulus conjugate to the quan- 
tity P at position r' with frequency to. Here Sais/s) = ±1 
specifies the symmetry property of a (/3) with respect 
to the T-operation. B suggests a magnetic field, but 
represents any time-reversal breaking field. In the case 
of a ferromagnet B can be associated with the magne- 
tization A/, the spontaneously generated time-reversal 
symmetry breaking field of a ferromagnet. The conduc- 
tivity tensor aab at a given frequency is proportional to 
the current-current response function. We can there- 
fore make the identification a — *■ Ja, P Jb, where Ji 
{i — X, y, z) are the components of the current operator. 
Since, = = — 1, we can conlcude that 



<yab{^\ B) = cr6a(w; -B). 



(3.3) 



Hence, we conclude that Gab is symmetric with respect 
to a and h in systems with T-symmetry. The anti- 
symmetric part aab{^) — (^bai^) is finite only if T- 
symmetry is broken. 

Now we turn back to irreversibility for the general form 
of the conductivity Gab, for which the dissipation is given 



by ( [Landau ar||1984 | 

ah 

= \ Y.^Re{aab + <Jba)MEaE;) 



Im((T<jf, - o-ba)Im(K£'b)]. 



(3.4) 



The real part of the symmetric combination aab + o'ba 
and the imaginary part of the antisymmetric combination 
Cah ~ CTfca contribute to the dissipation, while the imagi- 
nary part of (Jab + <Jba and the real part of Gab — o'ba rep- 
resent the dispersive (dissipationless) responses. There- 
fore, He^Gab — <Jba), which corresponds to the Hall re- 
sponse, does not produce dissipation. This means that 
the physical processes contributing to this quantity can 
he reversible, i.e., dissipationless, even though they are 
not necessarily so. This point is directly related to the 
controversy between intrinsic and extrinsic mechanism 
for the AHE. 

The origin of the ordinary Hall effect is the Lorentz 
force due to the magnetic field H: 



-v X H, 



(3.5) 



which produces an acceleration of the electron perpen- 
dicular to its velocity v and H. For free electrons 
this leads to circular cyclotron motion with frequency 
Lj(. — eH/imc). The Lorentz force leads to charge accu- 
mulations of opposite signs on the two edges of the sam- 
ple. In the steady state, the Lorentz force is balanced by 
the resultant transverse electric field, which is observed as 
the Hall voltage Vh. In the standard Boltzmann-thcory 
approach, the equation for the electron distribution func- 
tion f{p,x) is given by (Ziman, 19671: 



dl 
dt 



dl 
dx 



F 



dl 
dp 



dl 
dt 



(3.6) 



coll. 



where p and x are the momentum and real space coor- 
dinates, respectively. Putting F = —e{E+ - x H), and 
using the relaxatin time approximation for the collision 
term, — ^(/ — fcq), with /^q being the distribution func- 
tion in thermal equilibrium, one can obtain the steady 
state solution to this Boltzmann equation as 



/(p)=/eq(p)+5(p), 



with 



9{P) 



dfoq 

-re „ V 



de 



E 



H xE 



(3.7) 



(3.8) 



where order terms have been neglected and the mag- 
netic field Hz is assumed to be small, i.e. oJcT <C 1. The 
current density J = —e J (^^^^/(p) is obtained from 



Eq. ([3^8]) to order 0{E) and 0{EH) as 
J ^ (jE + (JhH X E, 



(3.9) 
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where 



1 

is the conductivity, while 



d^P 2 2 

(2^" ^ de 



cq 



(27r)3 



2 2 5/( 

7j e T 



er 



de 



(3.10) 



(3.11) 



is the ordinary Hall conductivity. 

Now let us fix the direction of the electric and magnetic 
fields a,s E — EyEy and H = H^ez- Then the Hall 
current is along the x-direction and we write it as 



Jx — (7HEy, 



(3.12) 



with the Hall conductivity an = ne^ t"^ H ^ / vn? c with n 
being the electron density. The Hall coefficient Rh is 
defined as the ratio Vh/JxHz, and from Eqs. (3.101 and 
(3.11 1, one obtains 



Rh — — - 



(3.13) 
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FIG. 33 Tight-binding model on a square lattice under 
the magnetic flux </> for each plaquette. We choose the 
gauge where the phase factor of the transfer integral along 



a;-direction is given by e 



with iy being the y-component 



of the lattice point position i = {i^ 



This result is useful to determine the electron density 
n experimentally, since it does not contain the relax- 
ation ti me r . Note here that for fixed electron density 
n, Eq. (3.131 means an oc Ha^^ and that pn esc H, in- 
dependent of the relaxation time r (or equivalently axx)- 
In other words, the ratio {cthI/o'xx = <^ 1. This 
relation will be compared with the AHE case below. 

What happens in the opposite limit lOcT 3> 1? In 
this limit, the cyclotron motion is completed many times 
within the lifetime r. This implies the closed cyclotron 
motion is repeated many times before being disrupted 
by scattering. When treated quantum mechanically the 
periodic classical motion leads to kinetic energy quantiza- 
tion and Landau level formation, which leads in turn to 
the celebrated quantum Hall effect. As is well known, 
in the 2D electron gas realized in semiconductor het- 
erostructures. Landau-level q uantization leads to t he cel - 
ebrated quantum Hall effect (Prange and Girvin 19871. 



If we consider the free electrons, i.e. completely neglect- 
ing the potential (both periodic and random) and the 
interaction among electrons, one can show that an is 

2 

given hy j^v where v is the filling factor of the Landau 
levels. In combination with electron localization, gaps 
at integer filling factors lead to quantization of a^. For 
electrons in a two-dimensional cyrstal, the Hall quantiza- 
tion is still quantized, a property which can be traced to 
the topological properties of Bloch state wavefunctions 
discussed below. CPrange and GirvinI |1987p. 



Now let us consider the magnetic field effect for elec- 
trons under the influence of the periodic potential. For 
simplicity we study the tight-binding model 



H 



(3.14) 



The magnetic field adds the Peierls phase factor to the 
transfer integral Uj between sites i and j, viz. 



t. 



iijexp[mjj], 



with 



he 



dr ■ A{r). 



(3.15) 



(3.16) 



Note that the phase factor is periodic with the period 2tt 
with respect to its exponent. Since the gauge transfor- 
mation 



Ci c^exp[i6i 
together with the redefinition 



(3.17) 



(3.18) 



itself is not a phys- 

a-jk + a.ki + au 



keeps the Hamiltonian invariant, 
ical quantity. Instead, the flux (p = aij 
per square plaquette is the key quantity in the problem 
(see Fig 33 ) 

For example, one may choose a gauge in which Uj 
along the directions ±y is a real number t, while tn^^ — 
texp[—i(f>iy] to produce a uniform flux distribution (j) in 
each square plaquette. Here iy is the y-component of the 
lattice point position i = {ix,iy)). The problem is that 
the ij,-dependence breaks the periodicity of the Hamil- 
tonian along the y-axis, which invalidates the Bloch the- 
orem. This corresponds to the fact that the vector po- 
tential A{r) = {—Hy, 0) for the uniform magnetic fleld 
H is y-dependent, so that momentum component Py is 
no longer conserved. However, when 0/(27r) is a rational 
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number n/m, where n and m are integers that do not 
share a common factor, one can enlarge the unit cell by 
m times along the y-direction to recover the periodicity, 
because exp[— Z(/)(jy-|-TO)] = exp[— Therefore, there 
appear m sub-bands in the Ist-Brillouin zone, each of 
which is characterized by a Chern number related to the 
Hall response as will be described in section HI.B (Thou- 



less et al. 19821. The message here is that an external 



commensurate magnetic field leads to a multiband struc- 
ture with an enlarged unit cell, leading to the quantized 
Hall response when the chemical potential is within the 
gap between the sub-bands. When 0/(27r) is an irrational 
number, on the other hand, one cannot define the Bloch 
wavefunction and the ele ctronic structure is described by 
the Hofstadter butterfly fHofstadterl 11976 1 . 



Since flux 2n is equivalent to zero flux, a commensu- 
rate magnetic field can be equivalent to a magnetic field 
distribution whose average fiux is zero. Namely, the to- 
tal fiux penetrating the m-plaquettes is 27rrt, which is 
equivalent to zero. Interestingly, spatially non-uniform 
flux distributions can also lead to quantum Hall effects 
This possibility was flrst considered by Haldane (Hal- 



dane 19881, who studied a tight-binding model on the 



honeycomb lattice. He introduced complex transfer inte- 
grals between the next-nearest neighbor sites in addition 
to the real one between the nearest-neighbor sites. The 
resultant Hamiltonian has translational symmetry with 
respect to the lattice vectors, and the Bloch wavefunction 
can be defined as the function of the crystal momentum 
k in the first Brillouin zone. The honeycomb lattice has 
two sites in the unit cell, and the tight-binding model 
produces two bands, separated by the band gap. The 
wavefunction of each band is characterized by the Berry 
phase curvature, which acts like a "magnetic field" in 
fc-space, as will be discussed in sections |III.B| and |IV.A[ 

The quantum Hall effect results when the Fermi en- 
ergy is within this band gap. Intuitively, this can be 
interpreted as follows. Even though the total fiux pen- 
etrating the unit cell is zero, there are loops in the unit 
cells enclosing a nonzero fiux. Each band picks up the 
fiux distribution along the loops with different weight 
and contributes to the Hall response. The sum of the 
Hall responses from all the bands, however, is zero as ex- 
pected. Therefore, the "polarization" of Hall responses 
between bands in a multiband system is a general and 
fundamental mechanism of the Hall response, which is 
distinct from the classical picture of the Lorentz force. 

There are several ways to realize a fiux distribution 
within a unit cell in momentum space. One is the rela- 
tivistic spin-orbit interaction given by 

he 

where V is the potential, s the spin, and p the momentum 
of the electron. This Hamiltonian can be written as 



H. 



SOI 



= Asoi ■ P, 



(3.20) 



with Asoi 



2^2^2 {sxW) acting as the effective vector 
potential. Therefore, in the magnetically ordered state. 



i.e., when s is ordered and can be regarded as a c-number, 
Asoi plf^ys a role similar to the vector potential of an 
external magnetic field. Note that the Bloch theorem is 
valid even in the presence of the spin-orbit interaction, 
since it preserves the translational symmetry of the lat- 
tice. However, the unit cell may contain more than two 
atoms and each atom may have multiple orbit als. Hence, 
the situation is very similar to that described above for 
the commensurate magnetic field or the Haldane model. 

The following simple model is instructive. Let us con- 
sider the tight-binding model on a square lattice given 
by 



H 



i,cr,a—x,y 

i,a.a—x.y 

^ ^ ^sp^ ifjVi-^-a.a.cT h.C. 
i^a^a—x^y 



(3.21) 

On each site, we put the 3 orbitals s, Px and py, asso- 
ciated with the corresponding creation and annihilation 
operators. The first three terms represent the transfer 
of electrons between the neighboring sites as shown by 
Fig.^ 
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The signs in front of t's are determined by the 
relative sign of the two orbitals connected by the trans- 
fer integrals, and all fs are assumed to be positive. The 
last term is a simplified SOI in which the z-component 
of the spin moment is coupled to that of the orbital mo- 
ment. This 6-band model can be reduced to a 2-band 
model in the ferromagnetic state when only the a = +1 
component is retained. By the spin-orbit interaction, the 
p-orbitals are split into 



1 



72 



{\Px)±i\Py)) 



(3.22) 



at each site, with the energy separation 2A. There- 
fore, when only the lower energy state |p— ) and the 
s-orbital \s) are considered, the tight- binding Hamilto- 
nian becomes H = J^k''^^ (^)^(^)'^(^) with tp{k) — 
[sk,a=i,Pk-.a=i]'^ and 



h{k) 



Eg — 2ts{cos kx + COS ky) \/2tsp{i sin kx + sin ky) 
\/2tsp{—i sin kx + sin ky) Ep + tp{cos kx + cos ky) 

(3.23) 

Note that the complex orbital |p— ) is responsible for 
the complex off-diagonal matrix elements of h{k). This 
produces the Hall response as one can easily see from the 



formula Eq. (1.5 1 and Eq. (1.6) in section LB When the 



Fermi energy is within the band gap, the Chern number 
for each band is ±1 when e^ — 4ts < Ep + 2tp, and they are 
zero otherwise. This can be understood by considering 
the effective Hamiltonian near fc = 0, which is given by 



h{k) 



(3.24) 



29 





i 








tsp 5 






* 


— ^cX-/ 


is 


r 




— > 


> 





FIG. 34 Tight-binding model on a square lattice for a 3 band 
model made from s and p^, —ipy orbitals with polarized spins. 
The transfer integrals between s and — ipy orbitals become 
complex along the y-direction in a way that is equivalent to 
effective magnetic flux 



with e = [{es — its) + {sp + 2ip)]/2 and m = (Es — 
Ats) — (Cp + 2tp), which is essentially the Dirac Hamil- 
tonian in (2+l)D. One can calculate the Berry curvature 
b±{k) = Tbz{k±)ez defined in Eq. ( T6 1 for the upper(+) 
and lower (— ) bands, respectively as 



2[kl 



,213/2 ■ 



(3.25) 



The integral over the two dimensional wavevector fcj^ = 
{kx, ky) leads to Hall conductance (Jh — |sgn(m)e^//i in 
this continuum model when only the lower band is fully 
occupied (Jackiw 1984). The distribution Eq. (3.251 in- 



dicates that the Berry curvature is enhanced when the 
gap m is small near this avoided band crossing. Note 
that the continuum model Eq. (3.24) cannot describe 



the value of the Hall conductance in the original tight- 
binding model since the information is not retained over 
all the first Brillouin zone. Actually, it has been proved 
that the Hall conductance is an integer multiple of e^/ft. 
due to the single-valudeness of the Bloch wavefunction 
within the first Brillouin zone. However, the change of 
the Hall conductance between positive and negative m 

can be described correctly. 

Onoda and Nagaosa (lOnoda and Nagaosal l2002|) 



demonstrated that a similar scenario emerges also for 
a 6-band tight binding model of t2g-orbitals with SOI. 
They showed that the Chern number of each band can 
be nonzero leading to a nonzero Hall response arising 
from the spontaneous magnetization. This suggests that 
the anomalous Hall effect can be of topological origin, 
a reinterpretation o f the intrinsic contribut i on fou nd by 
Karplus-Luttinger (Karplus and Luttinger 1954). An- 



other mechanism which can produce a flux is a non- 
coplanar spin structure with an associated spin chirality 
as has been discussed in section Hl.BI 

The considerations explained here have totally ne- 
glected the finite lifetime r due to the impurity and/or 



phonon scattering, a generalization of the parameter uJcT 
which appeared in the case of the external magnetic field. 
A full understanding of disorder effects usually requires 
the full power of the detailed microscopic theories re- 



viewed in Sec. IV The ideas explained in this section 
though are sufficient to understand why the intrinsic 
Berry phase contribution to the Hall effect can be so 
important. The SOI produces an effective vector poten- 
tial and a "magnetic flux distribution" within the unit 
cell. The strength of this flux density is usually much 
larger than the typical magnetic field strength available 
in the laboratory. In particular, the energy scale "/kJc" 
related to this flux is that of the spin-orbit interaction for 
Eq. (3.20), and hence of the order of 30 meV in 3d tran- 



sition metal elements. This corresponds to a magnetic 
field of the order of 300 T. Therefore, one can expect that 
the intrinsic AHE is easier to observe compared with the 
quantum Hall effect. 



B. Topological interpretation of the intrinsic mechanism: 
relation between Fermi sea and Fermi surface properties 

The discovery and subsequent investigation of the 
quantum Hall effect led to numerous important concep- 
tual advances (Prange and Girvin 19871. In particular. 



the utility of topological considerations in understanding 
electronic transport in solids was first appreciated. In 
this section, we provide an introduction to the applica- 
tion of topological notions such as the Berry phase and 
the Chern number to the AHE problem, and explain the 
long-unsuspected connections between these considera- 
tions and KL theory. A more elementary treatment is 
given in Ref . Ong a nd Lee[ |2006| 



The topological expression given in Eq. (1.5) was first 
obtained by Thouless, Kohmoto, Nightingale and Nijs 
(TKNN) (Thouless efai] [T982| for Bloch electrons in a 



2D insulating crystal lattice immersed in a strong mag- 
netic field H . These authors showed that each band is 
characterized by a topological integer called the Chern 
number 



' dky 



(3.26) 



According to Eq. (1.5), the Chern number C„ gives the 
quantized Hall conductance of an ideal 2D insulator, viz. 
Cxy = e^Cn/h. As mentioned in prior sections, these 
topological arguments were subsequently applied to semi- 



classical transport theory ( Sundaram and Niu 1999 1 and 



the AHE problem in itinerant ferromagnets (Jungwirth 



et al. 


2002b 


Onoda and Nagaosa 


2002 


1 , and the equiva- 


lence to the KL expression ( Karplus anc 


Luttinger 


19541 



for gh was established. (Note however, that in multiband 
systems, induced interband coherence among states with 
the same crystal momentum is a possibility.) 

Although the intrinsic Berry-phase effect involves the 
entire Fermi sea, as is clear from Eq. (1.5), it is usu- 



ally believed that transport properties of Fermi liquids 
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at low temperatures relative to the Fermi energy should 
be dependent only on Fermi-sruface properties. This is 
a simple consequence of the observation that at these 
temperatures only states near the Fermi surface can be 
excited to produce non-equillibirium transport. 

The apparent contradiction between conventional 
Fermi liquid theory ideas and the Berry phase theory 
of the AHE was resolved by Haldane (Haldanel 120041), 



showed that the Berry phase contribution to the AHE 
could be viewed in an alternate way. Because of its topo- 
logical nature, the intrinsic AHE, which is most naturally 
expressed as an integral over the occupied Fermi sea, can 
be rewritten as an integral over the Fermi surface. Let us 
start with the topological properties of the Berry phase. 
The Berry-phase curvature b„fc is gauge-invariant and 
divergence-free except at quantized monopole and anti- 
monopole sources with the quantum ±27r, i.e., 



qm6^{k-kni), = ±27r. (3.27) 



These monopoles and anti-monopoles appear at isolated 
k points in three dimensions. This is because in com- 
plex Hermitian eigenvalue problems, accidental degener- 
acy can occur by tuning the three parameters of fc. 

To gain further insight on the topological nature of 3D 
systems, it is useful to rewrite Eq. (1.5) as 



X 2^' 



(3.28) 



where 



K = J2k„ 



^ / d^kfienk) Kk- (3.30) 
Jf.b.z. 



We note that if the band dispersion e^k does not cross £p, 
is quantized in integer multiples (Cna) of a primitive 
reciprocal lattice vector Ga at T = 0. We have 



— / ^ CnaGa, 



(3.31) 



where the index a runs over the three independent prim- 



itive reciprocal lattice vectors (Kohmoto 1985) 



We next discuss the non-quantized part of the anoma- 
lous Hall conductivity. In a real material with multiple 
Fermi surface (FS) sheets indexed by a, we can describe 
each sheet by k'"p\s), where s — (s^,s^) is a parame- 
terization of the surface. It is convenient to redefine the 
Berry-phase connection and curvature as 



a\s) = a{k{s)) ■ dsik{s), 



bis) 



(3.32) 
(3.33) 



with eij the rank-2 antisymmetric tensor. 



By integrating Eq. (3.30) by parts, eliminating the in- 



tegration over the Brillouin zone boundary for each band. 



and dropping the band indices, we may write (at T = 0) 
the vector K in Eq. ( |3.29| as 



a a 

1 



(3.34) 



2tt 



ds^ hds^ b{s)k^p \s) 



+ ^y.Go.^ f ds-dis). (3.35) 

Here, Ca is the sum of Cna over the fully occupied bands, 
a labels sheets of the Fermi surface Sa, and (55^ is the 
intersection of the Fermi surface Sa with the Brillouin 
zone boundary i where k^p \s) jumps by Gai- Note that 
the quantity 



1 

2tt 



ds^ A rfs^ b{s) 



(3.36) 



gives an integer Chern number. Hence gauge invariance 
requires that 



(3.37) 



The second term in Eq. ( 3.35[ ) guarantees that is un- 
changed by any continuous deformation of the Brillouin 
zone into another primitive reciprocal cell. 



C. Early theoretical studies of the AHE 



(3.29) 1. Karplus-Luttinger theory and the intrinsic mechanism 



The pioneering work of Karplus and Luttinger (KL) 
( Karplus and Luttinger[ 1954 1 was the first theory of the 
AHE fully based on Bloch states ipnk- As a matter of 
course, their calculations uncovered the important role 
played by the mere existence of bands and the associated 
overlap of Bloch states. The KL theory neglects all lattice 
disorder, so that the Hall effect it predicts is based on an 
intrinsic mechanism. 

In a ferromagnet, the orbital motion of the itinerant 
electrons couples to spin ordering via the SOI. Hence all 
theories of the AHE invoke the SOI term, which, as we 
have mentioned, is described by the Hamiltonian given 
by Eq. |3.19| where V{r) is the lattice potential. The SOI 
term preserves lattice translation symmetry and we can 
therefore define spinors which satisfy Bloch's theorem. 
When SO interactions are included the Bloch Hamilto- 
nian acts in a direct product of orbital and spin-space. 
The total Hamiltonian Ht can be separated into contri- 
butions as follows: 



Ht — Hn 



Hsoi + He 



(3.38) 



where i/o is the Hamiltonian in the absence of SOI and 
He the perturbation due to the applied electric field E. 
In simple models the consequences of magnetic order can 
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be represented by replacing the spin operator in Eq.(3.19) 
by the ordered magnetic moment Ms as s ^ | ^ , where 
Mq is the magnitude of the saturated moment. Th e ma- 
trix element of He = —eEhXf, can be written as 

(3.39) 

where the "overlap" integral 



Jr (fc) = / dr ul^ir)—uMr) (3.40) 



are regular functions of k. In hindsight, J^" (fc) may 
be recognized as the Berry- phase connection a^i^k) dis- 
cussed in Sec. II.B.ll 

We next divide He into H]^ + H%, where H% {H^^) 
corresponds to the first (second) term on the right hand 
side of Eq. (3.391. Absorbing i?^ into the unperturbed 
Hamiltonian 



KL treated 



i.e., Hp — Hq + Hsoi 
the remaining term iJ^ as a perturbation. Accordingly, 
the density matrix p was written as 



p = poiHp) + pi 



(3.41) 



where po{Hp) is the finite-temperature equilibrium den- 
sity matrix and pi is the correction. KL assumed that pi 
gives only the ordinary conductivity, whereas the AHE 
arises solely from po^Hp). Evaluating the average veloc- 
ity Va as Va = Tr[pQVa], they found that 



where \F\ = {e/c)Hsov, m* is the effective mass, S is 
the number of incompletely filled d-orbitals, and p is the 
resistivity. 

Note that the p^-dependence implies that the off- 
diagonal Hall conductivity is independent of the 
transport lifetime r, i.e., it is well-defined even in the 
absence of disorder, in striking contrast with the diag- 
onal conductivity a. The implication that pn = cthP^ 
varies as the square of p was immediately subjected to 
extensive experimental tests, as described in Sec. |I.A| 

An important finding in the KL theory is that inter- 
band matrix elements of the current operator contribute 
significantly to the transport currents. This contrasts 
with conventional Boltzmann transport theory, where 
the current arises solely from the group velocity Vn,k ~ 
dSnM/dk. 

In the Bloch basis 



{r\iJJnk) 



Ah-', 



{r\Unk) 



(3.45) 



the TV-orbital Hamiltonian (for a given k) may be decom- 
posed into the 2N x 2N matrix h{k), viz. 



h{k) = {nk\h{k)\mk) al^{k)am{k). 



(3.46) 



The corresponding current operator for a given k is 



Jf,{t,k) = ^ (nfc 



dh{k) 
dhkn 



mk) a'l{k)am{k). (3.47) 



Va — — iei?fc ^ Pq(£'^^) J""(fc), (3.42) According to HcUman-Fcynman's theorem 



n,k 



where p'q is the derivative with respect to the energy. As 
it is clear here, this current is the dissipationless current 
in thermal equilibrium under the influence of the external 
electric field, i.e., H^. The AHE contribution arises from 
the interband coherence induced by an electronic field, 
and not from the more complicated rearrangements of 
states within the partially occupied bands. 

A second assumption of KL is that, in 3d metals, the 
SOI energy Hso £f (the Fermi energy) and W (the 
bandwidth), so that it suffices to consider Hgo to lead- 
ing order. Using Eq. (3.19), the AHE response is then 
proportional to jiV/sl, consistent with the empirical rela- 
tionship Eq. (2.1 1. More explicitly, to first-order in Hsoi, 
KL obtained 



k,n 



(3.43) 



where £„fc is the energy of the Bloch state for i^o, A is the 
averaged value of interband energy separation, and F„k 
is the force i{nk\[Hsoi,p]\'nk). This gives the anomalous 
Hall coefficient 



Rs 



2e'^Hso ^ I m 
mA^ 



(3.44) 



(3.48) 



the diagonal contribution to the current is (with n = m) 
Jl^''-^{t)^-eY,Vg{k)N^{k), (3.49) 



nk 



where Nn{k) = ajj(fc)a„(fc) is the occupation number 
in the state \nk). Eq. (3.491 corresponds to the ex- 
pression in Boltzmann transport theory mentioned be- 



fore. The inter-band matrix element (nk 



dh(k) 



ik) with 

n ^ m corresponds to inter-band transitions. As shown 
by KL, the inter-band matrix elements have profound 
consequences for the Hall current, as has become appar- 
ent from the Berry-phase approach. 



2. Extrinsic meclianisms 

a. Skew scattering 

In a series of reports, Smit (Smit 1955 19581 mounted 
a serious challenge to the basic findings of KL. In 
the linear-response transport regime, the steady-state 
current balances the acceleration of electrons by E 
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against momentum relaxation by scattering from impu- 
rities and/or phonons. Smit pointed out that this bal- 
ancing was entirely absent from the KL theory, and pur- 
ported to show that the KL term vanishes exactly. His 
reasoning wasa that the anomalous velocity central to 
KL's theory is proportional to the acceleration k which 
must vanish on average at steady state because the force 
from E cancels that from the impurity potential. 

More significantly, Smit proposed the skew-scattering 
mechanism (Fig.|3]) as the source of the AHE (Smit 1955 
1958). As discussed in Sec. LB. 2 in the presence of SOI, 



the matrix element of the impurity scattering potential 
reads 



{k's\V\k,s)=Vk.k' [6, 



((s'|cr|s) X k')-k) 



(3.50) 

Microscopic detailed balance would require that the tran- 
sition probability Wn^m between states n and m is 
identical to that proceeding in the opposite direction 
{Wm^n)- It holds, for example, in the Fermi's golden- 
rule approximation 



27r 



\{n\V\m)\^SiE„-Em), 



(3.51) 



where V is the perturbation inducing the transition. 
However, microscopic detailed balance is not generic. In 
calculations of the Hall conductivity, which involve the 
second Born approximation (third order in V), detailed 
balance already fails. In a simple N=l model, skew scat- 
tering can be represented by an asymmetric part of the 
transition probability 



kk' 



-T. ^kxk' ■ M, 



(3.52) 



When the asymmetric scattering processes is included 
(dubbed skew scattering), the scattering probability 
W{k — > k') is distinct from W{k' k). 

Physically, scattering of a carrier from an impurity in- 
troduces a momentum perpendicular to both the incident 
momentum k and the magnetization M. This leads to a 
transverse current proportional to the longitudinal cur- 
rent driven by E. Consequently, the Hall conductivity 
an and the conductivity a are both proportional to the 
transport lifetime r. Equivalently, pH = ohP^ is propor- 
tional to the resistivity p. 

As mentioned in Sec. |l]and in SecjlTj this prediction - 
qualitatively distinct from that in the KL theory - is con- 
sistent with experiments, especially on dilute Kondo sys- 
tems. These are systems which are realized by dissolving 
magnetic impurities (Fe, Mn or Cr) in the non-magnetic 
hosts Au or Cu. (At higher concentrations, these systems 
become spin glasses.) Although these systems do not ex- 
hibit magnetic ordering even at T as low as 0.1 K, their 
Hall profiles pn vs. H display an anomalous component 
derived from polarization of the magnetic local moments 
(Fig.|35]a). 

Empirically, the Hall coefficient Rh has the form 
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FIG. 35 The dependences of the Hall resistivity in AuFe on 
magnetic field (Panel a) and temperature (Panel b). [After 



Ref. Hurd 1972 



where Rq is the nominally T-independent OHE coeffi- 
cient and A is a constant (see Fig. 35 b). Identifying 



the second term with the paramagnetic magnetization 
M = xH, where xi^) oc 1/T is the Curie susceptibility 
of the local moments, we see that Eq. (3.531 is of the 
Pugh form Eq. ^J\ . 

Many groups have explored the case in which p and 
Ph can be tuned over a large range by changing the 
magnetic-impurity concentration c^. Hall experiments 
on these systems in the period 1970-1985 by and large 
confirmed the skew-scattering prediction p^ cx p. This 
led to the conclusion - often repeated in reviews - that 
the KL theory had been "experimentally disproved" . As 
mentioned in Sec. |lj this invalid conclusion ignores the 
singular role of TRI-breaking in ferromagnets. While 



Eq. (3.531 indeed describes skew scattering, the dilute 
Kondo system respects TRI in zero H. This essential 
qualitative differences between ferromagnets and systems 
with easily aligned local moments implies essential differ- 
ences in the physics of their Hall effects. 



b. Kondo theory 

Kondo has proposed a finite-temperature skew- 
scattering model in which spin waves of local moments at 
finite T lead to asymmetric scattering (Kondo) 19621. In 



Rh = Ro + A/T 



(3.53) 



the KL theory, the moments of the ferromagnetic state 
are itinerant: the electrons carrying the transport current 
also produce the magnetization. Kondo has considered 
the opposite limit in which nonmagnetic s electrons scat- 
ter from spin- wave excitations of the ordered c?-band local 
moments (with the interaction term JSn ■ s). Retaining 
terms linear in the spin-orbit coupling A and cubic in J, 
Kondo derived an AHE current that arises from tran- 
sition probabilities containing the skew-scattering term 
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k X k' . The T dependence of ph matches well (espe- 



cially near Tc) the pn vs. T profiles measured in Ni ( Jan 



1952 



1952 
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FIG. 36 Comparison of the curve of pH vs. T measured in 
Fe by Jan with Kondo's calculation [From Ref. [Kondo , 1962] ] 

Kondo has noted that a problem with his model is 
that it predicts that pH should vanish when the d orbital 
angular momentum is quenched as in Gd (whereas pn 
is observed to be large). A second problem is that the 
overall scale for pH (fixed by the exchange energies Fq, 
Fi and F2) is too small by a factor of 100 compared 
with experiment. Kondo's model with the s-d spin-spin 
interaction replaced by a (i(spin)-s(orbital) interaction 



has also been applied to antifcrromagnets (Maranzana 
|T967| ). 



c. Resonant skew scattering 

Resonant skew scattering arises from scattering of car- 
riers from virtual bound states in magnetic ions dissolved 
in a metallic host. Examples are XM, where X = Cu, 
Ag and Au and M = Mn, Cr and Fe. The prototypical 
model is a 3d magnetic ion embedded in a broad s band, 
as described by the Anderson model (Hewson 19931. As 



shown in Fig. |37| a spin-up s electron transiently occu- 
pies the spin-up bound state which lies slightly below 
ep- However, a second (spin-down) s electron cannot be 
captured because the on-site repulsion energy U raises 
its energy above ep- As seen in Fig. [37] the SOI causes 
the energies of the d orbitals (labelled by m) to be 
individually resolved. Here, cr = ± is the spin of the 
bound electron. The applied H merely serves to align 
spin a = + at each impurity. 

The scattering of an incident wave e^''"^ is expressed 
by the phase shifts Sl^{E) in the partial- wave expansion 
of the scattered wave. The phase shift is given by 



E) 



A 



(3.54) 



where A is the half-width of each orbital. 

In the absence of SOI, 5]^ {E) is independent of m and 
there is no Hall current. The splitting caused by SOI 
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FIG. 37 Sketch of the broadened virtual bound states of a 3d 
magnetic impurity dissolved in a non-magnetic metallic host. 
The SOI lifts the degeneracy of the d levels indexed by m, 
the magnetic quantum number. [From Ref. Fert and Jaoul| 
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results in a larger density of states at Ef for the orbital 
m compared with — m in the case where the up spin elec- 
tron is more than half-filled as shown in Fig. [37] Because 
the phase shift is sensitive to occupancy of the impurity 
state (Friedel sum-rule), we have 7^ ■ This leads 
to a right-left asymmetry in the scattering, and a large 
Hall current ensues. Physically, a conduction electron in- 
cident with positive z-component of angular momentum 
m hybridizes more strongly with the virtual bound state 
than one with negative — m. This results in more elec- 
trons being scattered to the left than the right. When 
the up spin density of states are filled less than half, i.e., 
E^^ > 0, the direction becomes the opposite, leading a 
sign change of the AHE. 

Explicitly, the splitting of i?^^ is given by 



El 



d± 



E"d± ± IrnXi, 



(3.55) 



wh ere is the SOI energy for spin a. Using Eq. (3.551 
in (3.54), we have, to order (A^r/A)^, 



cm rO ^aTTT- ■ ; 
(>2a = ('20^^^ Sm 



O9, 



2a 



A|m2 
4A2 



,3 rO 



Ccos<5L. 



(3.56) 



Inserting the phase shifts into the Boltzmann transport 
equation, we find that the Hall angle 70- ~ pZy/Pxx foi' 
spin a is, to leading order in Acr/A, 



3A„ 



crg^sin(2(5^^ - (5i)sin(5i, 



(3.57) 

where 61 is the phase shift of the p-wave channel which is 



assumed to be independent of ra and s (Fert et al. 1981 



Fert and Jaoul il972). Thermal averaging over a, we can 



obtain an estimate of the observed Hall angle 7. Its sign 
is given by the position of the energy level relative to ep- 
With the rough estimate Acr/A ~ 0.1, and sm&i ^ 0.1, 
we have I7I ~ 10"^. Without the resonant scattering 
enhancement, the typical value of 7 is ~ 10~^. 
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FIG. 38 The Hall coefficient Rh in heavy electron systems. 
Left panel: The T dependence of Rh in CeAla (open triangle), 
UPta (solid square), UAI2 (solid circle), and a single crystal 
of CeRu2Si2 (circle dot symbol). The field H is along the c- 
axis. Right panel: Interpretation of the Rh-T curve based on 
skew scattering from local moments. After Ref. |Hadzicleroux 
et al. 1986| and Lapierre et al. 1987[ [From Ref. Fert and 
Le7^[l987l ] 



Heavy-electron systems are characterized by a very 
large resistivity p above a coherence temperature Tcoh 
caused by scattering of carriers from strong fluctuations 
of the local moments formed by /-electrons at each lat- 
tice site. Below Tcoh, the local-moment fluctuations de- 
crease rapidly with incipient band-formation involving 
the /-electrons. At low T, the electrons form a Fermi 
liquid with a greatly enhanced effective mass. As shown 
in Fig. 38 the Hall coefficient Rh increases to a broad 
maximum at Tcoh before decreasing sharply to a small 
value in the low-T coherent-band regime. Resonant skew 
scattering has also been applied to account for the strong 
T depende nce of the Hall coefficient in CeCu2Si2, UBei3, 
andUPts (Coleman et al. 



and Levy 19871. 



1985 Fert et al. 1981 Fert 



d. Side-jump 

An extrinsic mechanism distinct from skew scattering 
is side-jump (Berger 1970[ ) (Fig. |3|. Berger considered 
the scattering of a Gaussian wavepacket from a spherical 
potential well of radius R given by 



V{r) 
V{r) 



2m 




(fc2 - kl) (r < R) 



in the presence of the SOI term Hsq = 
{l/2m^c^){r-^dV/dr)S,L,, where 5"^ (L,) is the 
z-component of the spin (orbital) angular momentum. 
For a wavepacket incident with wavevector k, Berger 
found that the wavepacket suffers a displacement Ay 
transverse to k given by 



(3.59) 



with Ac = h/mc the Compton wavelength. For k ^ kp = 
^jj-^ typical metals). Ay = 3 x 10^^^ m is far too 
small to be observed. 



In solids, however, the effective SOI is enhanced by 



band-structure effects by the factor (Fivaz 19691 



2m c 



2^2 



-T„ ^ 3.4 X 10*, 



(3.60) 



with T„ 



- {m*/3h^)j:„^^Jx^pyAE„^)\{m\q x p\n)\\ 
Here, AEnm — 0.5 eV is the gap between adjacent 
d-bands, x ~ 0-3 is the overlap integral, p = 2.5 x 
j^Q-io ^ jg nearest-neighbor distance, and ^ — 
-{ff /2m^c^){{r-^dV/dr)) = 0.1 eV is the atomic SOI 
energy. The factor in Eq. (3.601 is essentially the ra- 
tio of the electron rest mass energy mc^ to the energy 
gap AEnm- With this enhancement, the transverse dis- 
placement is Ay ~ 0.8 X 10""'^^ m, which renders the 
contribution relevant to the AHE. 

However, because the side-jump contribution to an is 
independent of t, it is experimentally difficult to distin- 
guish from the KL mechanism. We return to this issue 
in Sec. ITOI 



3. Kohn-Luttinger theory formalism 

a. Luttinger theory 

Partly motivated by the objections raised by Smit, 



Luttinger 


19581 


Kohn-Lutti 


nger 



ory ( Kohn and Luttinger 1957 Luttinger and Kohnj 



1958 1 . Employing a systematic expansion in terms of the 



impurity potential (p, he solved the transport equation for 
the density matrix and listed several contributions to the 
AHE current, including the intrinsic KL term, the skew- 
scattering term, and other contributions. They found 
that to zero-th order in f>, the average velocity is ob- 
tained as the sum of the 3 terms v'^f^^\up and v^^^ deflned 
by (iLuttingerl 119581) 



,(11) 



-ieEa 



up = —icEa 




k=0 



dka dkj3 



k=0 



(r > i?), (3.58) and 



Vp'^ — ieEoc 



d4_dJi 

dk„ dka 



(3.61) 



(3.62) 



(3.63) 



with jii the impurity concent ration. (For the definition 



of r 



(0) 



see Eq. (4.27) of Ref. Luttinger 

(11 



a/3 

The first term v 



19581 



^ is the skew-scattering contribu- 
tion, while the second U/ ; is the velocity obtained by KL 



(Karplus and Luttinger 
another term ot the sarnie' 



19541. The third term 



order as up. 



(b) 



IS 



The issues raised by comparing intrinsic vs. extrinsic 
AHE mechanisms involve several fundamental issues in 
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the theory of transport and quantum systems away from 
equUibrium. Following Smit (Smit 1955[ 1958) and Lut- 



tinger (Luttinger 19581, we consider the wavefunction 
expanded in terms of Bloch waves ipe^s, viz. 



(3.64) 



where i = (n, k) stands for the band index n and 
wavevector k. The corresponding expectation value of 
the position operator x/s is given by 



X0 



where J^'^ was defined in Eq. (3.401. Taking the 
quantum-statistical average using the density matrix 
{pT)e'.i — {0'i'0'*i), the expectation value can be written 
as 



dot. 



-a, + I 







(3.65) 



diPrk 
dk/B 



\k=k' + i '^{PT)i'iJp^ ■ (3.66) 



It may be seen that the second term in Eq. ( 3.66 1 does not 



contribute to the current because it is a regular function 
of fc, and the expectation value of its time-derivative is 
zero. Finally, one obtains the expression 



{V0) = I 



E 



dkg 



\k=k' 



by takin g the t ime-derivati ve of Eq.(3.66l 



Smit (Smit 



1955 1958) assumed o^) 



and obtained the diagonal part as 

dSnk 



i{pT)e. 



dki3 



(3.67) 



(3.68) 



while the off-diagonal part averages to zero be cause of 
the o scilla tory factor e^*'^^*"'^*')*. Inserting Eq. (3.68) in 
Eq. (3.67) and using {pt)u = one obtains the 



usual expression for the velocity, i.e., 

dSnk 



E(p^ 



dkp 



(3.69) 



which involves the group velocity, but not the anomalous 
velocity. A subtlety in the expression Eq. (3.67) is clar- 



ified by writing p — pq + fe where e is the adiabatic 
factor, and px — sfe^*. This leads to the expression 



h.nh' 



dkt 



(3.70) 



which approaches a finite value in the limit s ^ 0. This 
means that ^'''^gX^'"''' \k=k' ~ 1/s, i.e., a singular function 
of s. This is in sharp contrast to / itself, which is a 
regular function as s — > 0, which KL estimated. They 



considered the current or velocity operator instead of the 
position operator, the latter of which is unbounded and 
even ill-define d for periodic bou ndary conditions. 

Luttinger ( LuttingeT) 1958) argued that the ratio 
Vp'^'' /up equals ep/i-iniif), which may become less than 
1 for large impurity concentration Ui. However, if we as- 
sume if is comparable to ep, the expansion parameter is 
Eft/H, i.e., v'p^'^ /up = epr/h. 

The "metallicity" parameter spr/h plays a key role 
in modern quantum-transport theory, especially in the 
weak localization and interaction theory (Lee and Ra- 



makrishnan 1985). Metallic conduction corresponds to 
EpT/h ^ 1. More generally, if one assumes that the 
anomalous Hall conductivity <th is first order in the spin- 
orbit energy A, it can be written in a scaling form as 



e2 A J h 

(^H = T / 

ha £p \£pT 



(3.71) 



where the scahng function can be expanded as 

oo 

f{x) = (3.72) 

Here, c„ 's are constants of the order of unity. The leading 
order term c_ix~^ corresponds to the skew-scattering 
contribution oc r, while the second constant term cq is 
the contribution found by KL. If this expansion is valid, 
the intrinsic contribution by KL is always smaller than 
the skew scattering contribution in the metallic region 
with h/spT ^ 1. One should note however that the right 
hand side of this inequality should really have a smaller 
value to account for the weakness of the skew scattering 
amplitude. When both are comparable, i.e., h/epT ^ 1, 
one needs to worry about the localization effect and the 
system is nearly insulating, with a conduction that is of 
the hopping type. This issue is discussed in more detail 



in Sees. II.A lI.E and IV.D 



b. Adams-BJoujit formalism 
Adams and Blount (Adams and Blount 1959) ex- 



pressed the KL theory in a way that anticipates the mod- 
ern Berry-phase treatment by introducing the concept of 
"field-modified energy bands" and the "intracell" coordi- 
nate. They considered the diagonal part of the coordi- 
nate matrix in the band n as 



I— 



(3.73) 



which is analogous to Eq. (3.39) obtained by KL. The 



first term is the Wannier coordinate identifying the lat- 
tice site, while the second term - the "intracell" coor- 
dinate - locates the wavepacket centroid inside a unit 
cell. Significantly, the intracell nature of X™ implies 
that it involves virtual interband transitions. Although 
the motion of the wavepacket is confined to the conduc- 
tion band, its position inside a unit cell involves virtual 
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occupation of higher bands whose effects appear as a ge- 
ometric phase ( |Ong and Lee 2006 1 . They found that the 
curl Vfe X X"" acted Hke an effective magnetic field that 
lives in k space. 

The application of an electric field E leads to an 
anomalous (Luttinger) velocity, which gives a Hall cur- 
rent that is manifestly dissipationless. This is seen by 
evaluating the commutation relationship of cc^ and x'^. 
We have 



[x^ , x^] 



ih 



(3.74) 



where we have defined the field -B„(p) = Vp x 
This is analogous to the commutation relationship among 
the components of tt = p-\- {e/c)A in the presence of the 
vector potential, i.e., [tt^jTt^] — ihe^^\\S/r x ^('")]a = 
ihB\{r). The anomalous velocity arises from the ficti- 
tious "magnetic field" B„{p) (which lives in momentum 
space), and non-commutation of the gauge-covariant co- 
ordinates x^'s. This insight anticipated the modern idea 
of Berry phase curvature (see Sees, 
details) . 

Taking the commutator between and the Hamilto- 
nian ff„„ = En{p) — F'^x'^, we obtain 

t;„„ - -i[x„ i/„„] = ^^1^ -Fx S„(p). (3.75) 

The second term on the right hand side is called the 
anomalous velocity. Adams and Blount reproduced the 



III.B and IV.A for 



results of Luttinger ( Luttinger 1958 1 and demonstrated 
it for the simple case of a uniform field Bn{p) = D. 



As recognized by Smit (Smit 1955 19581, the currents 



associated with the anomalous velocities driven by E and 
by the impurity potential mutually cancel in steady state. 
However, D introduces corrections to the "driving term" 
and the "scattering term" in the transport equation. As 
a consequence, the current J is 



3 dE I 



Ft 

m 



FxD ( t'^P \ 



\ 3mTA J 



{F X D) 



.(3.76) 



Note that the current arises entirely from the average 
of the "normal current" . However, the second term in 
Eq. (3.761 is similar to the anomalous velocity and is 



consistent with the conclusions of KL ( Karplus and Lut 



tinger, 19541 and Luttinger (Luttinger 



19581 



The con- 
sensus now is that the KL contribution exists. However, 
in the clean limit r — > oo, the leading contribution to 
the AHE conductivity comes from the skew-scattering 
term. A more complete discussion of the semiclassical 
treatment is in Sec. IIV.AI 



Nozieres-Lewiner theory 



Adopting the premises of Luttinger's theory (Lut- 
1958[), Nozieres and Lewiner (Nozieres and 



tmger 



Lewiner 



1973 1 investigated a simplified model comprised 



of one conduction and one valence band to derive all the 
possible contributions to the AHE. The Fermi level ep 
was assumed to lie near the bottom of the conduction 
band. Integrating over the valence band states, Nozieres 
and Lewiner derived an effective Hamiltonian for a state 
k in the conduction band. The derived position operator 
is rcff = r + p, where the new term p which involves the 
SOI parameter A is given by 



-Xk X S. 



(3.77) 



This polarization or effective shift modifies the transport 
equation to produce several contributions to ajj . In their 
simple model, the anomalous Hall current J ah e (besides 
that from the skew scattering) is 



J AHE = 2Ne^XE X S, 



(3.78) 



where N is the carrier concentration, and S is the aver- 
aged spin polarization. The sign is opposite to that of the 
intrinsic term Jintrinsic — ~2N \E x S obtained from 
the SOI in an ideal lattice. They identified all the terms 
as arising from the spin-orbit correction to the scattering 
potential, i.e., in their theory there is no contribution 
from the intrinsic mechanism. However, one should be 
cautious about this statement because these cancellations 
occur among the various contributions and such cancel- 
lations can be traced back to the fact that the Berry cur- 
vature is independent of k. In particular, it is often the 
case that the intrinsic contribution survives even when 
/i lies inside an energy gap, i.e., (the side-jump 

contribution is zero in this case). This issue lies at the 
heart of the discussion of the topological aspect of the 
AHE. 



IV. LINEAR TRANSPORT THEORIES OF THE AHE 

In this section we describe the three linear response 
theories now used to describe the AHE. All three theories 
are formally equivalent in the ei?r 3> 1 limit. The three 
theories make nearly identical predictions; the small dif- 
ferences between the revised semiclassical theory and the 
two formalism of the quantum theory are thoroughly un- 
derstood at least for several different toy model systems. 
The three different approaches have relative advantages 
and disadvantages. Considering all three provides a more 
nuanced picture of AHE physics. Their relative corre- 
spondence has been shown analytically in several simple 

[2007t . The gen- 



models (Sinitsyn 2008 Sinitsyn et al 



eralized semiclassical Boltzmann transport theory which 
takes the Berry phase into account is reviewed in Sec. 
|IV.A| this theory has the advantage of greater physi- 
cal transparency, but it lacks the systematic character of 
the microscopic quantum-mechanical theories whose ma- 
chinery deals automatically with the problems of inter- 
band coherence. Another limitation of the semiclassical 
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Boltzmann approach is that the "quantum correction" to 
the conductivities, i.e., the higher order terms in H/epT, 
which lead to the Anderson localization and other inter- 
esting phenomena (Lee and Ramakrishnan 1985[ ) cannot 
be treated systematically. 

The microscopic quantum-mechanical linear response 
theories based on the Kubo formalism and the Keldysh 
(non-equiliribrum Green's function) formalism are re- 
viewed in Sec. |IV.B.l| and Sec. |IV.C| respectively. These 
formulations of transport theory are organized differently 
but are essentially equivalent in the linear regime. In 
the Kubo formalism, which is formulated in terms of the 
equilibrium Green's functions, the intrinsic contribution 
is more readily calculated, especially when combined with 
first-principles electronic structure theory in applications 
to complex materials. The Keldysh formalism can more 
easily account for finite lifetime quantum scattering ef- 
fects and take account of the broadened quasiparticle 
spectral features due to the self-energy, while at the same 
time maintaining a structure more similar to that of semi- 
classical transport theory (Onoda et al. 2006a I. AWe 



contrast these techniques by comparing their applications 
to a common model, the ferromagnetic Rashba model in 
two-dimensions (2D) which has been studied intensively 
in recent years (Sec |IV.D I. 



A. Semiclassical Boltzmann approach 

As should be clear from the complex phenomenology of 
the AHE in various materials classes discussed in Sec. |TT] 
it is not easy to establish a one- size- fits- all theory for 
this phenomenom. From a microscopic point of view the 
AHE is a formidable beast. In subsequent sections we 
outline systematic theories of the AHE in metallic sys- 
tems which employ Keldysh and Kubo linear response 
theory formalisms and are organized around an expan- 
sion in disorder strength, characterized by the dimen- 
sionlcss quantity hjepT. A small value for this parame- 
ter may be taken as a definition of the good metal. We 
start with semiclassical theory, however, because of its 
greater physical transparency. In this section we out- 



line the modern version (Sinitsyn 2008 Sinitsyn et al. 



2007) of the semiclassical transport theory of the AHE. 



This theory augments standard semiclassical transport 
theory by accounting for coherent band-mixing by the 
external electric field (which leads to the anomalous ve- 
locity contribution) and by a random disorder potential 
(which leads to side jump). For simple models in which 
a comparison is possible, the two theories (i.e. (i) the 
Boltzmann theory with all side-jump effects (formulated 
now in a proper gauge-invariant way) and the anoma- 
lous velocity contribution included and (ii) the metallic 
limit of the more systematic Keldysh and Kubo formal- 
ism treatments) give identical results for the AHE. This 
section provides a more compact version of the material 



presented in the excellent review by Sinitsyn (Sinitsyn 



tail. Below we take fi. = 1 to simplify notation. 

In semiclassical transport theory one retreats from a 
microscopic description in terms of delocalized Bloch 
states to a formulation of transport in terms of the dy- 
namics of wavepackcts of Bloch states with a well defined 
band momentum and position {n,kc,rc) and scattering 
between Bloch states due to disorder. The dynamics of 
the wavepackets between collisions can be treated by an 
effective Lagrangian formalism. The wavepacket distri- 
bution function is assummed to obey a classical Boltz- 
mann equation which in a spatially uniform system takes 
the form ( Sundaram and NTu[ |1999 I : 



dt 



hkr 



djj_ 
dkr 



= - '^{^I'jfl - ^1,1' II')- 



(4.1) 



Here the label I is a composition of band and momenta 
{n,kc) labels and fli'^i, the disorder averaged scattering 
rate between wavepackets defined by states I and I' , is to 
be evaluated fully quantum mechanically. The semiclas- 
sical description is useful in clarifying the physical mean- 
ing and origin of the different mechanism contributing to 
the AHE. However, as explained in this review, contribu- 
tions to the AHE which are important in a relative sense 
often arise from inter-band coherence effects which are 
neglected in conventional transport theory. The tradi- 
tional Boltzmann equation therefore requires elaboration 
in order to achieve a successful desc ription of the AHE. 

There is a substantial li t eratu re ( Berger 1970 Jung 
|wirth et al. 2002b Smit 19551 on the application of 
Boltzmann equation concepts to AHE theory (see Sec 



III.Cl. However stress was often placed only on one of 



the several possible mechanisms, creating a lot of confu- 
sion. A cohesive picture has been lacking until recently 



( |Sinitsyn[ [2008, .Sinitsyn et aZ.[ [2006, ,2005, ). In particu- 
lar, a key problem with some prior theory was that it in- 
correctly ascribed physical meaning to gauge dependent 
quantities. In order to build the correct gauge-invariant 
semiclassical theory of AHE we must take the following 
steps (|Sinitsyn[[2008|: 



i) obtain the equations of motion for a wavepacket 
constructed from spin-orbit coupled Bloch elec- 
trons, 

ii) derive the effect of scattering of a wave-packet 
from a smooth impurity, yielding the correct gauge- 
invariant expression for the corresponding side- 
jump, 

iii) use the equations of motion and the scattering 
rates in Eq. (4.1 1 and solve for the non-equillibirum 



2008 ) , to which we refer readers interested in further de- 



distribution function, carefully accounting for the 
points at which modifications are required to ac- 
count for side-jump, 

iv) utilize the non-equilibrium distribution function to 
calculate the dc anomalous Hall currents, again ac- 
counting for the contribution of side-jump to the 
macroscopic current. 
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The validity of this approach is partially established in 
the following sections by direct comparison with fully mi- 
croscopic calculations for simple model systems in which 
we are able to identify each semiclassically defined mech- 
anisms with a specific paxt of the microscopic calcula- 
tions. 



1. Equation of motion of Bloch states wave-packets 

We begin by defining a wavepacket centered at position 
Tc with average momentum kc'. 



k 

A key aspect of this wavepacket is that the complex func- 
tion, sharply peaked around kc, must have a very specific 
phase factor in order to have the wavepacket centered 
around r^- This can b e shown to be ( jMarder) |2000( |Sun-| 
daram and Niul 119991): 



Wk^,rA^) = |wfc^,^^(fc)| exp[i(fc - fee) • a„], (4.3) 

where a„ = {unk\idkUnk) is the Berry's connection of the 
Bloch state (see Sec. LB I. We can generate dynamics for 



the wavepacket parameters fc^, and by constructing 
a semiclassical Lagrangian from the quantum wavcfunc- 
tions: 

d 

hkc ■ Tc + hkc ■ an{kc) - E{kc) + eF(rc).(4.4) 



C 



All the terms in the above Lagrangian are common to 
conventional semiclassical theory except for the second 
term, which is a geometric term in phase space depending 
only on the path of the trajectory in this space. This term 
is the origin of the momentum-space Berry phase ( (Berry 



1984) effects in anomalous transport in the semiclassica 



formalism. The corresponding Euler-Lagrange equations 
of motion are: 



hkr 



-eE 
dkc 



hkc X bji(^kc 



(4.5) 
(4.6) 



where 6„(fcc) = V x a„ is the Berry's curvature of 
the Bloch state. Compared to the usual dynamic equa- 
tions for wave-packets formed by free electrons, a new 
term emerges due to the non-zero Berry's curvature of 
the Bloch states. This term, which is already linear 
in electric field E, is of the Hall type and as such will 
give rise in the linear transport regime to a Hall cur- 
rent contribution from the entire Fermi sea, i.e., j^^^n — 
-e^Ex ^j:^fo{Enk)bn{k). 

2. Scattering and the side-jump 

From the theory of elastic scattering we know that the 
transition rate w/^c in Eq. (4.1 1 is given by the T-matrix 



element of the disorder potential: 

u;vi = 27r\Tn\''6{ei. -ei). 



(4.7) 



The scattering T-matrix is defined by T/'j = {l'\V\^pl), 
where V is the impurity potential operator and [ipi) is 
the eigenstate of the full Hamiltonian H — Ho + V 
that satisfies the Lippman-Schwinger equation 1-0;) = 
|/)-|-(e; — Hq + ir|)~^V\^lJl). j?/';) is the state which evolves 
adiabatically from \l) when the disorder potential is 
turned on slowly. For weak disorder one can approxi- 
mate the scattering state \^i) by a truncated series in 
powers of V;/' = {l\V\l'): 



l^;>«IO+E 



Vi"i 
' ei" H 



r 



(4.8) 



Using this expression in the above definition of the T- 
matrix and substituting it into Eq. (4.7 1, one can expand 



the scattering rate in powers of the disorder strength 



ujii, +ujii, ••• , (4.9) 



,(3) 



,(4) 



where uj 



(2) 



(3) 



2^^{\Vu,\^)d^sSiel-el,), 



ui;,,' = 27r I > ' '" '' '.^"'^ -I- c.c.j 6{ei -ev), 

(4.10) 

and so on. 

We can always decompose the scattering rate into 
components that are symmetric and anti-symmetric in 



the state indices: 



,(s/a) 



{lohi ± LJiii)/2. In conven- 



tional Boltzmann theory the AHE is due solely to the 



anti-symmetric contribution to the scattering rate ( Smit 
T955|. 



The physics of this contribution to the AHE is quite 
similar to that of the longitudinal conductivity. In par- 
ticular, the Hall conductivity it leads to is proportional 
to the Bloch state lifetime r. Since W;,; is symmetric, the 

leading contribution to uj^f^ appears at order . Partly 
for this reason the skew scattering AHE conductivity con- 
tributions is always much smaller than the longitudinal 
conductivity. (It is this property which motivates the 
identification below of additional transport mechanism 
which contribute to the AHE and can be analyzed in 
semiclassical terms.) The symmetric part of tjj\f} is not 
essential since it only renormalizes the second order re- 

(2) 

suit for Lolii and the antisymmetric is given by 

4'"^ = -{2TTfY.5i.^i - tl")lIn{Vu'Vvl"Vlnl)d^s5{^l - ev). 

(4.11) 

This term is proportional to the density of scatterers, n^. 
Skew scattering has usually been associated directly with 
^(3a)^ neglecting in particular the higher order term 
which is proportional to nf and should not be disregarded 
because it gives a contribution to the AHE which is of 
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the same order as the side-jump contribution considered 
below. This is a common mistake in the semiclassical 



analyses of the anomalous Hall effect ( Sinitsyn 2008 1 



Now we come to the interesting side-jump story. Be- 
cause the skew scattering conductivity is small, we have 
to include effects which are absent in conventional Boltz- 
mann transport theory. Remarkably it is possible to pro- 
vide a successful analysis of one of the main additional 
effects, the side-jump correction, by means of a careful 
semiclassical analysis. As we have mentioned previously 
side jump refers to the microscopic displacement Sri^ii 
experienced by a wave-packets formed from a spin-orbit 
coupled Bloch states, when scattering from state I to 
state I' under the influence of a disorder potential. In 
the presence of an external electric field side jump leads 
to an energy shift AUij' = —eE ■ Sri jr. Since we are 
assuming only elastic scattering, an upward shift in po- 
tential energy requires a downward shift in band energy 
and vice-versa. We therefore need to adjust Eq. (4.1 1 by 

9/o(ei), 



adding J^u ' 



-eE ■ Sri II to the r.h.s. 



'''' . . . 

But what about the side-jump itself? An expression for 

the side-jump Sri^i associated with a particular transi- 
tion can be derived by integrating through a transition 
(ISinitsynl 120081) . We can write 



dk'w{k)w*{k')e- {re''' '') <fc,(r) 

w„fc(r)e'(^('=')"^('=»*/''| = + (4.12) 

J dkr 



d_ 
~dt 



dr I dkw* {k)unk{r) 
I J 



d 

i^w{k)unk{.r) 

Ok 



This expression is equivalent to the equations of motion 
derived within the Lagrangian formalism; this form has 
the advantage of making it apparent that in scattering 
from state I to a state a shift in the center of mass 
coordinate will accompany the velocity deflection. From 
Eq. (4.121 it appears that the scattering shift will go ap- 



proximately as: 
Sriii « 



(4.13) 



This quantity has usually been associated with the side- 
jump, although it is gauge-dependent and therefore arbi- 
trary in value. The correct expression for the side jump 
is similar to this one, at least for the smooth impurity po- 
tentials situation, but was derived only recently by Sinit- 



syn et al. (Sinitsyn et al. 2006 20051: 



Srm = {ui'\i-^ui>) 



d_ 
dk' 



{uiV^m: 



- Dk' ,k&T:g[{uv\ui)], 
(4.14) 

where arg[a] is the phase of the complex number a and 
Dk' ,k — + -gk,- The last term is essential and makes 
the expression for the resulting side-jump gauge invari- 



ant. Note that the side-jump is independent of the de- 
tails of the impurity potential or of the scattering pro- 
cess. As this discussion shows, the side jump contribu- 
tion to motion during a scattering event is analogous to 
the anomalous velocity contribution to wave-packet evo- 
lution between collisions, with the role of the disorder 
potential in the former case taken over in the latter case 
by the external electric field. 



3. Kinetic equation for the semiclasscial Boltzmann distribution 



Equations (4.7 1 and (4.14) contain the quantum me- 



chanical information necessary to write down a semiclass- 
cial Boltzmann equation that takes into account both 
the change of momentum and the coordinate shift dur- 
ing scattering in the presence of a driving electric field E. 
Keeping only terms up to linear order in the electric field 



the Boltzmann equation reads ( Sinitsyn et al. 2006 1 



9h 

dt 



eE ■ Vol 



= - E If I - fi' - 



eE ■ Sri', 



(4.15) 

where v^i is the usual group velocity Vqi = dei/dk. Note 
that for elastic scattering we do not need to take account 
of the Pauli blocking which yields factors like /;(1 — /;/) 
on the r.h.s. of Eq. (4.151, and that the collision terms 
are linear in // as a consequence. (For further discus- 



sion of this point see Appendix B in Ref. Luttinger and 
Kohn 1955| ) This Boltzmann equation has the standard 



form except for the coordinate shift contribution to the 
collision integral explained above. Because of the side- 
jump effect, the collision term does not vanish when the 
occupation probabilities // are replaced by their ther- 
mal equilibrium values when an external electric field is 
present: Uei)- foiei-eE-Srw) « -^^eE-Sri'i + 0. 
Note that the term containing uji^'fi should be written 



1,1 



If I. In making this simplification we are 
imagining a typical simple model in which the scattering 
rate depends only on the angle between k and k'. In 

that case, tOi'f, = and we can ignore a complication 
which is primarily notational. 

The next step in the Boltzmann theory is to solve for 
the non-equilibrium distribution function // to leading 
order in the external electric field. We linearize by writ- 
ing // as the sum of the equilibrium distribution /o(e;) 
and non-equilibrium corrections: 



fi^Mei)+gi + g'i 



.dist 



(4.16) 



where we split the non-equilibrium contribution into two 
terms gi and g^'^'^^ in order to capture the skew scatter- 
ing effect, gi and g^'^^^^ solve independent self-consistent 



time-independent equations ( (Sinitsyn et al. 2006j): 



eE 



(4.17) 
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and 

E (sr*^* - gf'-' - ^-^^E ■ <5n,) = 0. 

(4.18) 

In Eg. (4.141 we have noted that Srin — Srw. To solve 
Eq. (4.17) we further decompose gi = gf + gf^ + gf^ so 
that 

Y.u^\l''\gt - gt) + Y,J;^){gf - gf) = 0, (4.19) 



The first term is the intrinsic contribution which should 
be by now familiar to the reader: 



(4.24) 



Next are the effects due to coordinate shifts during scat- 
tering events (for E along the y-axis): 



(4.25) 



Y.4r{gt - at) + Y.^^){9f gf) = 0. (4.20) 



Here is the usual diagonal non-equilibrium distribu- 
tions function wh ich ca n be shown to be proportional to 
n^^. From Eq. (4.111 and Lolf,"''' ^ ni, it foll ows t hat 

nf and Eq. (4.201, it 



9i 



3a 



Finally, from (jj[ff' 



follows that gf ~ n°; illustrating the dangers of ignoring 



the 



,(4a) 



contribution to ui /' . One can also show from 



Eq. (4.181 that g^ 



4. Anomalous velocities, anomalous Hall currents, and 
anomalous Hall mechanisms 

We are now at the final stage where we use 
the non-equilibrium distribution function derived from 



Eqs. (4.17 1-(4.20 1 to compute the anomalous Hall cur- 
rent. To do so we need first to account for all contri- 
butions to the velocity of semiclassical particles that are 
consistent with this generalized semiclassical Boltzmann 
analysis. In addition to the band state group velocity 
Vqi = defdk, we must also take into account the veloc- 
ity contribution due to the accumulations of coordinate 
shifts after many scattering events, another way in which 
the side-jump effect enters the theory, and the veloc- 
ity contribution from coherent band mixing by the elec- 



trie field (the anomalous velocity effect) (|Nozieres and 


Lewiner| 1973 Sinitsyn 


2008 


Sinitsyn et al. 


2006 


): 



Vl = 



dk 



bix eE+y ^uJinSrvi. 



I' 



Combining Eqs. (4.161 and (4.211 we obtain the total 
current 



follows from the distribution function correction due to 
side jumps while 



eY^{gi/ Ey)^ujvi{5ri,i)a 
I v 



(4.26) 



is the current due to the side-jump velocity, i.e., due to 
the accumulation of coordinate shifts after many scatter- 
ing events. Since coordinate shifts are responsible both 
for cr^y*'*' and for cr^-^, there is, unsurprisingly, an inti- 
mate relationship between those two contributions. In 
most of the literature, cr°?** is usually considered to be 



idist 
xy 



part of the side-jump contribution, i.e., a\ 
a^^. We distinguish between the two because they are 
physically distinct and appear as separate contributions 
in the microscopic formulation of the AHE theory. 
Finally, a^y^ and cr^^^"*-' are contributions ari sing from 



the asymmetric part of the collision integral (Sinitsyn 
200l: 



T;^'^^eJ2igf/E,){voi), 



7 



-eY,igf/E,){voi)y 



(4.27) 



(4.28) 



According to the old definition of skew scattering both 
could be viewed as skew scattering contributions because 
they originate from the asymmetric part of the collision 
term (Smit 19551. However, if instead we define skew- 
scattering as the contribution proportional to n~^, i.e. 
linear in r, as in Sec. |I.B| it is only the first contribu- 



(4.21) tion, Eq. (4.27), which is the skew scattering (Leroux- 



Hugon and Ghazali 1972 Luttinger 19581. The second 



J = e E fivi = e E(/o(e;) + + gf + gf + gf '"^) give 2D Dirac band ^[sinitsyn et~al 



contribution, Eq. (4.281, was generally discarded in prior 
semiclassical theories, although it is parametrically of the 
same size as the side-jump conductivity. Explicit quanti- 

so far exist only for the mas- 
In parsing 



tative estimates of cr^^^ 



20071 



(^ eE + ^ uJviSrm) 



(4.22) 



This gives five non-zero contributions to the AHE up to 
linear order in E: 



total „int 
fj = (T 

xy xy 



adist 



ski 



sk2~sj 
'xy 



(4.23) 



this AHE contribution we will incorporate it within the 
family of side-jump effects due to the fact that it pro- 
portional to n^, i.e. independent of a^x- However, it is 
important to note that this "side-jump" contribution has 
no phyiscal link to the side-step experienced by a semi- 
classcial quasiparticle upon scattering. An alternative 
terminology for this contribution is intrinsic skew scat- 
tering to distinguish its physical origin from side-jump 
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deflections (Sinitsyn 20081, but, to avoid further confu- 



sion, we simply include it as a contribution to the Hall 
conductivity which is of the order of and originates 
from scattering. 

As we will see below, when connecting the microscopic 
formalism to the semiclasscial one, a"y, can be directly 
identified with the single bubble (Kubo formalism) con- 



adist 
V 



-<f con- 



tribution to the conductivity, cr° 

stitute the usually termed ladder-diagram vertex correc 
tions to the conductivity due to scattering and there- 
fore it is natural to group them together although their 



Looking more closely tr^^ we notice that every term 
depends on products of retarded Green's functions only 
or products of advanced Green's functions only. It can 
be shown that only the disorder free part of cr^^ is impor- 
tant in the weak disorder limit, i.e. , this contribution 
is zeroth order in the parameter 1/kplsc. The only ef- 
fect of disorder on this contribution (for metals) is to 
broaden the Green's functions (see below) through the 
introduction of a finite lifetime (Sinitsyn et al. 



20071 



It can therefore be shown by a similar argument that in 
general a^y, is of order l/ kphc and can be neglected in 
physical origins are distinct. ct|^^ is identified directly the weak scattering limit (Mahan 19901. Thus, impor- 
with the three-scattering diagram used in the literature 
This comparison has been made specifica lly for two sim- 
ple models, the massive 2D Dirac band (Sinitsyn et al. 



2007 ) and the 2D Rashba with exchange model ( |Borunda 



et al. 2007 1 



B. Kubo formalism 



1. Kubo technique for the AHE 



The Kubo formalism relates the conductivity to the 



quilibrium current-current correlation function (Kubo 



1957). It provides a fully quantum mechanical formally 



exact expression for the conductivity in linear response 
theory (Mahan 1990). We do not review the formal ma- 



chinery for this approach here since can be found in many 
textbooks. Instead, emphasize the key issues in studying 
the AHE within this formalism and how it relates to the 
semiclassical formalism described in the previous section. 

For the purpose of studying the AHE it is best to re- 
formulate the current-current Kubo formula for the con- 
ductivity in the form of the Bastin formula (see appendix 
A in ( jCrepieux and Bruno 2001 1 ) which can be manipu- 
lated into the more familiar form for the conductivity of 
the Kubo-Streda formula for the T = Hall conductivity 



tant disorder effects beyond simple quasiparticle lifetime 
broadening are contained only in cr^". For these reasons, 
it is standard within the Kubo formalism to neglect tJ^J^ 
and evaluate the cr^^ contribution with a simple lifetime 
broadening approximation to the Green's function. 

Within this formalism the effect of disorder on the 
disorder-configuration averaged Green's function is cap- 
tured by the use of the T-matrix, defined by the integral 
equation T = W + WGqT, where W = Y.r Vo5{r - u) is 
a delta-scatterers potential and Gq are the Green's func- 
tion of the pure lattice. From this one obtains 

G Go + GqTGo = Go + GqSG. (4.32) 

Upon disorder averaging we obtain 

E = {W)c + {WGoW)c + {WGoWGqW)c + ... (4.33) 

To linear order in the impurity concentration, tii, this 
translates to 

E(z, fc) = n,Vk,k + ^ I] Vk,k'Go{.k',z)Vk'M + ■■■ , 

(4.34) 

with Vfc.fc' = V{k — k') being the Fourier transform of the 
single impurity potential, which in the case of delta scat- 
terers is simply Vq (see Fig. 
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lib) 



(T^y where: 



-iy^ = ^,Mi'.G%F)vyG^{ep)),, (4.29) 



2TrV 



Jib) 



Tr(«,G^(eF){)j,G^(eF)+i)»G^(eF)%G^(eF)), 

(4.30) 



2 /- + 00 C^(f:\ 

oil = At / de/(e)TrKG«(e)t;,'^ 



' xy 



oo 



de 



G^ie) 
de 



VyG^ie) + c.c] 



(4.31) 



for a graphical reprenta- 
tion). Note that G and Gq are diagonal in momentum 
but, due to the presence of spin-orbit coupling, non- 
diagonal in spin-index in the Pauli spin-basis. Hence, 
the blue lines depicted in Fig. 39 represent G and are in 
general matrices in band levels. 

One effect of disorder on the anomalous Hall conduc- 
tivity is taken into account by inserting the disorder av- 
eraged Green's function, G^ /^, d irectly int o the expres- 
'sions for a^"^ and cr^^, Eq. |4.29| and Eq. ( |4.3l| . This 
step captures the intrinsic contribution to the AHE and 
the effect of disorder on it, which is generally weak in 
metallic systems. This contribution is separately identi- 
fied in Fig. [39] a. 

The so-called ladder diagram vertex corrections, also 



Here the subscript c indicates a disorder configuration 
average. The last contribution, a^y, was first derived by 
Streda in the c ontext of studying the quantum Hall effect 
( |Streda|[l982l ). In these equations G^/^(ef) = (e^-iJib 
i5)~^ are the retarded and advanced Green's functions 
evaluated at the Fermi energy of the total Hamiltonian. 



separately identified in Fig. 39 contribute to the AHE 
at the same order in l/kpl as the intrinsic contribution. 
It is useful to define a ladder-diagram corrected velocty 
vertex Vai^p) = Va + Svaiep), where 



n.Vr 



V 



Y,G%f){vc. + 6v^{eF))G^{ep), 



(4.35) 
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(a) 



intrinsic 



"side-jump" 



"sl<ew-scattering" 



(b) 



.r>+g>+...= r>+r> 



FIG. 39 Graphical representation of the Kubo formalism ap- 
plication to the AHE. The solid blue lines are the disorder 
averaged Green's function, G, the red circles the bare velocity 
vertex = dHo / dhka , and the dashed lines with crosses rep- 
resent disorder scattering {riiVQ for the delta-scatter model), 
(b) Svy is the velocity vertex renormalized by vertex correc- 
tions. 



as depicted in Fig. 39 b. Note again that Vaiep) and Va ~ 
dHo/dhka are matrices in the spin-orbit coupled band 
basis. The skew scattering contributions are obtained by 
evaluating, without doing an infinite partial sum as in 
the case of the ladder diagrams, third order processes in 
the disorder scattering shown in Fig. [39] 

As may seem obvious from the above machinery, calcu- 
lating the intrinsic contribution is not very difficult, while 
calculating the full effects of the disorder in a systematic 
way (beyond calculating a few diagrams) is challenging 
for any disorder model beyond the simple delta-scattering 
model. 

Next we illustrate the full use of this formalism for the 
simplest nontrivial model, massive 2D Dirac fermions, 
with the goals of illustrating the complexities present 
in each contribution to the AHE and the equivalence 
of quantum and semiclassical approaches. This model 
is of course not directly linked to any real material re- 
viewed in Sec. |TT] and its main merit is the possibility of 
obtaining full simple analytical expressions for each of 
the contributions. A more realistic model of 2D fermions 
with Rashba spin-orbit coupling will be discussed in Sec. 
|IV.C[ The ferromagnetic Rashba model has been used to 
propose a minimal model of AHE for materials in which 



band crossing near the Fermi surface dominate the AHE 



physics (Onoda et al. 2008) 



The massive 2D Dirac fermion model is specified by: 

Ha^v{k^(Tj. + kyay) + /:!^(7^ + Vdis, (4.36) 

where Vdis — X]i^o'5(r ~ (^x and CFy are Pauli 

matrices and the impurity free spectrum is = 
±a/A^ -I- {vkY where k — |k| and the labels ± distin- 
guish bands with positive and negative energies. We 
ignore in this simple model spin-orbit coupled disorder 
contributions which can be directly incorporated through 
similar calculations as in Crepieux et al (Crepieux and 



Bruno 2001). Within this model the disordered aver- 



aged Green's function is 



{eF-t++iT+){eF-e- 



+iT-) ' 

(4.37) 

± - Foi 



where F = -KruVg / [Av^) , Fi ^ Fco s(6>) 
FiC0s(6'), and cos^ = A/^/{vk)' 

this disorder model r cx 1 /ui . Using the result in 



7 = -Lo 
A^. Note that within 



Eq. (4.37) one can calculate the ladder diagram correc- 



tion to the bare velocity vertex given by Eq. (4.35) 



(1 



cos 



^ A(l + 3cos2 6')2 



(1 



3 cos^ ( 
(4.38) 

where 6 is evaluated at the Fermi energy. The details of 
the calculation of this vertex correction is described in 
Appendix A of Ref . |Sinitsyn et al.\ |2007[ Incorporating 
this result in a^y we obtain the intrinsic and side-jump 
contributions to the conductivity for e^? > A 



' xy 



2-nhV 



J2k Tr [vcr^GvayG] 



X^fe Tr {va^G(>VyG\ 

;^ cos B I ; 



3 sin"' 9 I 4 sin"' 9 
1-1-3 cos^e) (1-1-3 cos2 e)2 



(4.39) 



(4.40) 



The direct calculation of the skew-scattering diagrams of 



Fig. 
resu 



39 



IS a. 



sk 



TTs given by 



' 2-iThnVo (iA^ + ivkFyy- 



and the final total 



'xy 



47rri^(DfeF)2 + A2 ' 
3(vkF)* 



(4A2 + (l,feF)2)2 



4{vkFr 
iA^ + {vkF)'^ 

(vkF)'^A 



2TThnVo (4A2 + (i.A;j^)2)2 



(4.41) 



2. Relation between tlie Kubo and the semiclassical formalisms 

When comparing the semiclassical formalism to the 
Kubo formalism one has to keep in mind that in the 



semiclassical formalism the natural basis is the one that 
diagonalize the spin-orbit coupled Hamiltonian. In the 
case of the 2D massive Dirac model this is sometimes 
called chiral basis in the literature. On the other hand 
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xy 



\ 



..■<0>v 



intrinsic contribution _ 



xy 



side jump 



_sk2—sj 

r 

xy 



skew scattering 



FIG. 40 Graphical representation of the AHE conductivity in the chiral (band eigenstate) basis. The two bands of the two- 
dimensional Dirac model are labeled "±". The subset of diagrams that correspond to specific terms in the semiclassical 
Boltzmann formalism are indicated. 



in the application of the Kubo formalism it is simplest to 
compute the different Green's functions and vertex cor- 
rections in the Pauli basis and take the trace at the end 
of the calculation. In the case of the above model one 
can apply the formalism of Sec. |IV.A| and obtain the 
following results for the five distinct contributions: (j"y, 



adist sj 
^ xy ' ^ xy 7 



^sfc2- 



^-^ and cr^^^ . 



Below we quote the results 
for the complicated semiclassical calculation of each term 
(see Sec. IV of Sinitsyn et al. (Sinitsyn et al. 2007) for 
details) : 



47rfiv/A2 + (wfcj.)2' 



SI adtst 
(J — (J 

^ xy ^ xy 



sk-sj 



e^Akl 



27rVfcf + A2(fc2 + 4A2) ' 



'xy 



47rfi^(wfcF)2-t- A2[4A2 + {vkpYY ' 



(4.42) 



(4.43) 



(4.44) 



proportional to the density of scatterers. However, recall 
that in Sec. |I.B| we have defined the side-jump contri- 
bution as the disorder contributions of zeroth order in 
rii, i.e., T*^, as opposed to being directly linked to a side- 
step in the scattering process in the semiclassical theory. 
Hence, it is the sum of the three physically distinct pro- 
cesses a^y^'^* + af^y 
identical to Eq. 
tion. Therefore 



(4.40 



crj^2 ^iiich can be shown to be 



the o 



after some algebraic manipula- 
d notion of associating the skew 
scattering directly with the asymmetric part of the col- 
lision integral and the side-jump with the side-step scat- 
tering alone leads to contradictions with their usual as- 
sociation with respect to the dependence on t (or equiv- 
alenty 1/nj). We also note, that unlike what happens in 
simple models where the Berry's curvature is a constant 
in momentum space, e.g., the standard model for elec- 



trons in a 3D semiconductor conduction band (Nozieres 



and Lewiner 19731, the dependence of the intrinsic and 



side-jump contributions are quite different with respect 
to parameters such as Fermi energy, exchange splitting, 
etc. 



(z;/cf)*A 



27rfin,t/o (4A2 + {vkpYf 



(4.45) 



The correspondance with the Kubo formalism results 
can be seen after a few algebra steps. The contribu- 
tions cr™ * and are equal in both cases (Eq. (4.39) 
and Eq. ( 4.41[ )). As expected, the intrinsic contribution. 



C. Keldysh formalism 

Keldysh has developed a Green's function formalism 
applicable even to the nonequilibrium quantum states, 
for which the diagram techniques based on Wick's the- 



cr*^ , is independent of disorder in the weak scattering 



and Baym 


1962 


and Smith 


1986 



|Keldysh[ |1965| |Mahan[TT990l [Rammer 

. Unlike with the thermal (Matsubara) 

limit and the skew scattering contribution is inversely Green's functions, the Keldysh Green's functions are de- 
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fined for any quantum state. The price for this flexibility 
is that one needs to introduce the path-ordered prod- 
uct for the contour from t = —oo — )■ i = oo and back 
again from t = oo — *■ — oo. Correspondingly, four kinds 
of the Green's functions G^, G^, and G^ need to be 



considered, although only three are independent (Baym 



and Kadanoffj 


1961 


Kadanoff and Baym 


|1962||Keldysh| 


1965 


Mahan 


1990 


Rammer and Smith 


1986 


). There- 



fore, the diagram technique and the Dyson equation for 
the Green's function have a matrix form. 

In linear response theory, one can use the usual thermal 
Green's function and Kubo formalism. Since approxima- 
tions are normally required in treating disordered sys- 
tems, it is important to make them in a way which at 
least satisfies gauge invariance. In both formalisms this 
is an important theoretical requirement which requires 
some care. Roughly speaking, in the Keldysh formal- 
ism, G^ and G^ describe the single particle states, while 
represents the non-equilibrium particle occupation 
distribution and contains vertex corrections. Therefore, 
the self-energy and vertex corrections can be treated in a 
unified way by solving the matrix Dyson equation. This 
facilitates the analysis of some models especially when 
multiple bands are involved. 

Another and more essential advantage of Keldysh for- 
malism over the semiclassical formalism is that one can go 
beyond a finite order perturbative treatment of impurity 
scattering strength by solving a self-consistent equation, 
as will be discussed in the next subsection. In essence, 
we are assigning a finite spectral width to the semiclassi- 
cal wave packet to account for an important consequence 
of quantum scattering effects. In the Keldysh formalism, 
the semiclassical limit corresponds to ignoring the history 
of scattering particles by keeping the two time labels in 
the Greens functions identical. 

We restrict ourselves below to the steady and uniform 
solution. For more generic cases of electromagnetic fields, 
see Ref. Sugimoto et al. 2007 Let x = {t,x) be the 



time-space coordinate. Green's functions depend on two 
space-time points xi and X2, and the matrix Dyson equa- 
tion for the translationally invariant system reads ( Ram- 
mer and Smith[|1986| ): 



(e-F(p)-S(e,p))®G(£,p) = 1, 
G{e,p)®{e-H{p)-±{e,p)) - 1, 



(4.46) 



where we have changed the set of variables {xi;x2) to 
the center-of-mass and the relative coordinates, and then 
proceeded to the Wigner representation {X;p) by means 
of the Fourier transformation of the relative coordinate; 



{X,x) = 



X\ -f X2 



,Xi~ X2 



dt / dx S<^t-P-^)/f' . 



(4.47) 

withp — {e,p). In this Dyson equation, the product (g) lis 
reserved for matrix products in band indices, like those 
that also appear in the Kubo formalism. 

In the presence of the external electromagnetic field 



A^, we must introduce the mechanical or kinetic energy- 
momentum variable 



Trf,{X;p)=Pf, + eA^{X). 



(4.48) 



replacing p as the argument of the Green's function, as 
shown by Onoda et al ( jOnoda et al. 2006b) . In this rep- 
resentation, G^(X; 7r)/27ri is the quantum-mechanical 
generalization of the semi-classical distribution function. 
When an external electric field E is present, the equa- 
tion of motion, or equivalently, the Dyson equation, re- 
tains the same form as Eq. (4.461 when the product 



is replaced by the so-called Moyal product (Moyal 1949 



Onoda et al. 



exp 



2006b I given by 

^£;.(Kv, 



(4.49) 



Henceforth, d and d denote the derivatives operating 
on the right-hand and left-hand sides, respectively, and 
the symbol p = {s,p) is used to represent the mechanical 
energy-momentum tt. In this formalism, only gauge in- 
variant quantities appear. For example, the electric field 
E appears instead of the vector potential A. 

Expanding Eq. (4.491 in E and inserting the result in 



Eq. (4.46 ), one obtains the Dyson equation to linear order 
in E, corresponding to linear response theory. The linear 
order terms G^ and in E are decomposed into two 
parts as 



G< 
E — 



jd,f{s) + (g^ 



^E 

si 



/(£), 

/(£)■ 



(4.50) 
(4.51) 



Here, j(e) represents the Fermi distribution function. In 
these decompositions, the first term on the r.h.s. cor- 
responds to the nonequilibrium deviation of the distri- 
bution function due to the electric field E. The second 
term, on the other hand, represents the change in quan- 
tum mechanical wavefunctions due to E, and arises due 
to the multiband effect (Haug and Jauho 1996) through 



the noncommutative nature of the matrices. 

The corresponding separation of conductivity contri- 



butions is: fj,,- ~ CT, 



^11 



with 



(T,, = e h / ,^ J I 1 . Tr 



P 



{2TTh) 



Tr 



v,{p)G%jip)\ dj{e), (4.52) 
\{P) {G%{p)-G%{p))]m. 



(4.53) 



This is in the same spirit as the Streda version (Stfeda 



19821 of the Kubo-Bastin formula (Bastin et al. 1971 



Kubo , 1957 ) . The advantage here is that we can use the 
diagrammatic technique to connect the self-energy and 
the Green's function. For dilute impurities, one can take 
the series of diagrams shown in Fig. [41] corresponding 
to the T-matrix approximation. In this approximation. 
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the self-consistent integral equation for the self-energy 
and Green's function can be solved and the solution used 



to evaluate the first and second terms in Eq. (4.501 and 
Eq. Kbll. 



X"™ 



£(e) 



5$ "„ 



S(e.P') G(e,p') G(e.p") G(e,p') G(E.p") G(E,p- ) 



where cr = ± labels the two eigenstates \p, a) at momen- 
tum p. 

With this we can then obtain the intrisic contribution 
to the AHE by integrating over occupied states at zero 
temperature (Culcer et al. 2003 Dugaev et al. 20051: 



AH-int 
'xy 



2h 



Ao 



9{^l~eAp^)), (4.57) 



FIG. 41 Diagrammatic representation of the self-energy in 
the self-consistent T-matrix approximation in the Keldysh 
space, which is composed of the infinite series of multiple 
Born scattering amplitudes. [From Ref. Onoda et al. 2008 ] 



together w ith the 
and (Onoda 



In general, Eqs. (4.521 and (4.531. 
self-consistent equations for G"',G^ 
\et al., 2008) defines a systematic diagrammatic methoc 
for ca lculating in the Streda dec omposition ( S tfeda 



1982) of the Kubo-Bastin formula (Bastin et at 1971 



Kubo 19571 



D. Two-dimensional ferromagnetic Rashba model - a 
minimal model 

A very useful model to study fundamental aspects 
of the AHE is the ferromagnetic two-dimensional (2D) 
Rashba model (Bychkov and Rashba 1984): 



where p± denotes the Fermi momentum for the band a = 
±. 

An important feature of cr^j^^*"* is its enhancement 
in the interval of eo.+ < ^J^ < £o,-, where it approaches a 
maximum value close to e^/2/i, no matter how small Aq 
is (provided that it is larger than h/r). Near p = the 
Berry curvatures of the two bands are large and opposite 
in sign. The large Berry curvatures translate into large 
intrinsic Hall conductivities only when the chemical po- 
tential lies between the local maximum of one band and 
the local minimum of the other. This enhancement of 
the intrinsic AHE near avoided band crossings is illus- 
trated in Fig. |42]b, where it is seen to survive moderate 
disorder broadening of several times Aq. This peaked fea- 
ture arises from the topological nature of cr^j^~™*. As 
a consequence it is important to note that the result is 
non-perturbative in SOI; only a perturbative expansion 
on h/spT is justified. 



Hip)t 



P 

2m 



\p X & ■ e'' - Aqo-^ + V{x). (4.54) by Dugaev et al. ( [Dugaev et al. 



Cucler et al. ( Culcer et all 2003 ) were the first to study 



this model, obtai ning Eq. ([4.57[ ). They were was followed 

2005[ ) where the intrinsic 



Here m is the electron mass, A is the Rashba spin-orbit 
interaction strength, Aq is the mean field exchange split- 
ting, <T = ((Ta;, (Tj,, (Tz) and (T° are the Pauli and identity 
matrices, is the unit vector in the z direction, and 
V{x) = Vq J2i S{r — ri) is a 5-scatterer impurity potential 
with impurity density rii . Quantum transport properties 
of this simple but non-trivial model have been intensively 
studied in order to understand fundamental properties of 
the AHE in itinerant metallic ferromagnets. The metal- 
lic Rashba is a simple, but its AHE has both intrinsic 
and extrinsic contributions and both minority and ma- 
jority spin Fermi surfaces. It therefore captures most of 
the features that are important in real materials with a 
minimum of complicating detail. The model has there- 



contribution was calculated within the Kubo formalism. 
Although these studies found a non-zero cr^j^~™*, they 
did not calculate all the contributions arising from disor 
der (some aspects of the disorder treatment in (Dugaev 
et al. \ [2005 ) where corrected by these authors in (Sinit 



syn et al. 



2007|). 

T 



The intrinsic AHE a 



AH—int 



comes 

form both a^y and a^y. The first part a^y contains the 



intra-band contribution a. 



1(a) 
xy 



which is sensitive to the 



fore received a lot of attention (Borunda et al. 


^OOT: 
,,2007j 


Dugaev et al. 


2005 [Inoue et al.\ 2006 Kato et al. 


Kovalev et al. 


20091 20081 INunner et al.\ 2007 


Onoda 


et al. 2006b[ 2008). 



impurity scattering vertex correction. 

The calculation of Ox^y^ incorporating the effects of dis- 
order using the Kubo formalism, i.e., incorporating the 
ladder vertex corrections ("side-jump") and the leading 
0{Vq) skew-scattering contributions (Sec. IV.B.l I, yields 
a vanishing for the case where ep '^s, above the gap 
at p = {i.e. both subbands are occupied), irrespec- 
tive of the strength of the spin-independent scattering 
ampUtude ( Borunda et al. 2007 Inoue et olj [2006 Nun 



ner et al. 



The bare Hamiltonian has the band dispersion: 



20081). On the other hand, when only the 
is occ upied, aj^y^ is g i ven b y the 



ea{p) 



p 

2m 



aAp, 



majority band a = 

skew-scattering contribution ( [Borunda et al. 20071, 



A2p2 



A2 

^0' 



(4.55) 



1 X^pXD+ifi)AoA, 



illustrated in Fig. 42 a, and Berry-phase curvature 

X^h'^Aoa 



h ^T'inip^imp 



(3A2- 



A2 )2 



(4.58) 



bUp)^h' [VpX{i{p,a\Vp\p,a))Y 



2A3 



(4.56) 



in the leading-order in (l/njmp). 

Some properties of the ferromagnetic Rashba model 
appear unphyiscal in the limit r — > cx) (Onoda et al 
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FIG. 42 (a) Band dispersion of the ferromagnetic 2D Rashba 
model in the clean limit, (b) The intrinsic anomalous Hall 
conductivity of the Hamiltonian Eq. ( |4.54[ ) as functions of 
the Fermi level ep measured from the bottom of the majority 
band and Ti/t = mrii^pVo / ti^ , which is the Born scattering 
amplitude for A = Aq — 0. This figure is for the parameter 
set Ao — 0.1, 2mA = 3.59, and 2mVo/ fi'^ = 0.6, with energy 
unit has been taken as £o,-. [From Ref. Onoda et al. 2006b 



2008 ) : Gxy vanishes discontinuously as the chemical po- 
tential ^ crosses the edge eo,- of the minority band 
which leads to a diverging anomalous Nernst effect at 
fi — Eq^-, irrespective of the scattering strength, if one as- 
sumes the Mott relation to be valid for anomalous trans- 
port ( Smrcka and Stf eda [1977 1 . However, this unphys- 
ical property does not really hold, in fact, axy does not 
vanish even when both subbands are occupied, as shown 
by including all higher-order Born scattering amplitudes 
as it is done automatically in the numerical Keldysh ap- 
proach (Onoda et al. 2006a I. In particular, the skew- 



scattering contribution arises from the odd-order Born 
scattering (i.e., even order in the impurity potential) be- 
yond the conventional level of approximation, 0{Vq), 
that gives rise to the normal skew scattering contribu- 

This yields the uncon- 
(X l/riimp independent of 



tions (Kovalev et al. 



ventional behavior a: 



20081 

— sKeu 



Vo ( jKovalev et al. 20081. The possible appearance of 
the AHE in the case where both subbands are occupied 
was also suggested in the nume rical diagonalizatio n cal- 
culation of the Kubo formula ( Kato et al. , |2007 1 . The 
influence of spin-dependent impurities has also been an- 



alyzed (Nunner et al. 20081 



A numerical calculation of based on the Keldysh 
formalism using the self-consistent T-matrix approxima- 




cyj[e/{2nfi)] 



FIG. 43 Total anomalous Hall conductivity vs. a^x for 
the Hamiltonian Eq. ( |4.54[ ) obtained in the self-consistent 



r-matrix approximation to the Keldysh approach (Kovalev 
et al.\ pOOO] [Onoda et al] |2006b[ [2TO8|. Curves are for a 
variety of disorder strengths. The same parameter values 
have been taken as in Fig. |42] with the chemical potential 
being located at the center of the two subbands. The dashed 
curves represent the corresponding semiclassical results. [Af- 
ter Ref. [Kovalev et al\ [20091 ] 



the AHE as a function of cr^x at low temperatures ( Onoda 
et al.\ 2006a 20081. In particular, it shows a crossover 
from the predominant skew-scattering region in the clean 
limit {cTxy cx axx) to an intrinsic-dominated metallic re- 
gion {(7xy ~ constant). In this simple model no well 
defined plateau is observed. These results also suggests 
another crossover to a regime, referred to as the inco- 

where <t„ 



2008), 



'xy 



de- 



herent regime by (Onoda et al. 
cays with the disorder following the scaling relation of 
Gxy oc Gxx with n « 1.6. This scaling arises in the cal- 
culation due to the influence of finite-lifetime disorder 
broadening on cr^^^"*"*, while the skew-scattering contri- 
bution is quickly dimished by disorder as expected. Ko- 
valev et al. ( [Kovalev et al. . 2009 1 revisited the Keldysh 
calculations for this model, studying them numerically 
and analytically. In particular, their study extended 
the calculations to include the dependence of the skew- 
scattering contribution on the chemical potential both 
for ji < eo,- and for fi > Eq,- (as shown in Fig. 43 ). These 
authors demonstrated that changing the sign of the im- 
purity potential changes the sign of the skew-scattering 
contribution. The data collapse illustrated in Fig. 43 



tion, shown in Fig. 43 suggests three distinct regimes for 



then fails, especially near the intrinsic-extrinsic crossover. 
Data collapse in a^^ vs. cTxx plots is therefore not a gen- 
eral property of 2D Rashba models and should not be 
expected in real materials. There is, however, a suffi- 
cient tendency in this direction to motivate analyzing 
experiments by plotting data in this way. 

Minimal model- The above results have the following 
implications for the generic nature of the AHE (Onoda 
et a?.[|2006b[|2(308| . The 2D ferromagnetic Rashba model 
can be viewed as a minimal model that takes into ac- 
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count both the "parity anomaly" (Jackiw 1984) asso- 



ciated with an avoided-crossing of dispersing bands, as 
well as impurity scattering in a system with two Fermi 
surface sheets. Consider then a general 3D ferromag- 
net. When the SOI is neglected, majority and minority 
spin Fermi surfaces will intersect along lines in 3D. For 
a particular projection of Bloch momentum along the 
magnetization direction, the Fermi surfaces will touch at 
points. SOI will generically lift the band degeneracy at 
these points. At the fc-point where the energy gap is a 
minimum the contribution to crf^,^~*"* will be a maxi- 
mum and therefore the region around this point should 
account for most of , as the first-prin c iples calcula- 



tions seem to indicate ( Fang et at. 2003 Wang et at. 
2007[ |2006[ |Yao et al.\ |2004[ [200 7') We emphasize that 
the Berry phase contributions from the two bands are 
nearly opposite, so that a large contribution from this 
region of k-space accrues only if the Fermi level lies be- 
tween the split bands. Expanding the Hamiltonian at 
this particular fc-point, one can then hope to obtain an 
effective Hamiltonian of the form of the ferromagnetic 
Rashba model. Note that the gap Aq in this effective 
model Hamiltonian, Eq. (4.541, comes physically from 



the SOI splitting at this particular fe-point, while the 
"Rashba SOI" A is proportional to the Fermi velocity 
near this crossing point. In 3D, the anomalous Hall con- 
ductivity is then given by the 2D contribution integrated 
over pz near this minimum gap region, and remains of 
order of /ha (a the lattice constant) if no accidental 
cancellation occurs. 

We point out that a key assumption made in the above 
reasoning is that the effective Hamiltonian obtained in 
this expansion is sufficiently similar to the ferromagnetic 
Rashba model. We note that in the Rashba model it 
is Aq and not the SOI which opens the gap, so this 
is already one key difference. The SO interactions in 
any effective Hamiltonian of this type should in general 
contain at least Rashba-like and Dresselhaus-like contri- 
butions. Further studies examining the crossing points 
more closely near these minimum gap regions will shed 
further light on the relationship between the AHE in real 
materials and the AHE in simple models for which de- 
tailed perturbative studies are feasible. 

The minimal model outlined above suggests the pres- 
ence of three regimes: (i) the superclean regime domi- 
nated by the skew-scattering contribution over the intrin- 
sic one, (ii) the intrinsic metallic regime where Gxy be- 
comes more or less insensitive to the scattering strength 
and a^xi and (iii) the dirty regime with kp(. = 2epT /h < 
1 exhibiting a sublinear dependence of pxy o: p^^", or 
equivalently axy cx: cr"^ with n w 1.6. It is impor- 
tant to note that this minimal model is based on elas- 
tic scattering and cannot explain the scaling observed 
in the localized hopping conduction regime as (Txx is 
tuned by changing T. Nevertheless, if we multiply the 
2D anomalous Hall conductivity by a^^ with the lat- 
tice constant a ^ 5 A for comparison to the experi- 
mental results on three-dimensional bulk samples, then 



an enhanced of the order of the quantized value 

e'^/h in the intrinsic regime should be interpreted as 
~ /ha ~ 10'^ cm~^. This value compares well with 
the empirically observed cross over seen in the experimen- 



tal findings on Fe and Co, as discussed in Sec. II. A 



V. CONCLUSIONS: FUTURE PROBLEMS AND 
PERSPECTIVES ON AHE 

In this concluding section, we summarize what has 
been achieved by the recent studies of the AHE and what 
is not yet understood, pointing out possible directions for 
future research on this fascinating phenomenon. To keep 
this section brief, we exclude the historical summary pre- 



sented in Sec. LA and Sec. III.C which outline the early 
debate on the origin of the AHE. We avoid repeating all 
thel points highlighted already in Sec. |I.A| and focus on 
the most salient ones. Citations are kept to a minimum 
as well and we refer mostly to the sections in which the 
material was presented. 



a. Recent developments 

The renewed interest in the AHE, which has lead to 
a richer and more cohesive understanding of the prob- 
lem, began in 1998 and was fueled by other connected 
developments in solid state physics. These were: i) the 
development of geometrical and topological concepts use- 
ful in understanding electronic properties such as quan- 
tum phase interference and the quantum Hall effect. (Lee 
and Ramakrishnan 1985 Prange and Girvin 19871, ii) 



the demostration of the close relation between the Hall 
conductance and the toplogical Chern number revealed 



by the TKNN formula (Thouless et al. 19821, iii) the 



development of accurate first-principles band structure 
calculation which account realistically for SOI, and iv) 
the association of the Berry-phase concept (Berry 1984 1 



with the noncoplanar spin configuration proposed in the 
context of the resonating valence bond (RVB) theory of 



cuprate high temperature superconductors (Lee et 
20061. 



Intrinsic AHE- The concept of an intrinsic AHE, debated 
for a long time, was brought back to the forefront of the 
AHE problem because of studies which successfully con- 
nected the topological properties of the quantum states 
of m atter a nd the transport Hall response of a system. In 
Sec. 



and 



I.B.I 



we have defined cr^j^ 



both experimentally 
theoretically. From the latter, it is rather straight- 
forward to write cr^j^^~™* in terms of the Berry curvature 
in the fc-space, from which the topological nature of the 
intrinsic AHE can be easily recognized immediately. The 



topological non-perturbative quality of cr^ 



AH-int 



is high- 



lighted by the finding that for simple models with spon- 
taneous magnetization and SOI, bands can have nonzero 
Chern numbers even without an external magnetic field 
present. This means that expansion with respect to SOI 
strength is sometimes dangerous since it lifts the degen- 
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eracy between the up- and down-spin bands, leading to 
avoided band crossings which can invahdate such expan- 
sion. 

Even though the interpretations of the AHE in real 
systems are still subtle and complicated, the view that 
can be the dominant contribution to a^,^ in 
certain regimes has been strengthened by recent compar- 
isons of experiment and theory. The intrinsic AHE can 
be calculated from first-principles calculations or, in the 
case of semiconductors, using k ■ p theory. These calcu- 
lations have been compared to recent experimental mea- 
sureme nts fo r several materials such as Sri_a;Caa;Ru03 
(section II.Bl, Fe (section II.Al, CuCr2Se4_xBr2: (section 



TLDJ, and dilute magnetic semiconductors (section II. C I 
The calculations and experiments show semi-quantitative 
agreement. More importantly however, violations of the 
empirical relation an oc M have been established both 
theoretically and experimentally. This suggests that the 
intrinsic contribution has some relevance to the observed 
AHE. On the other hand, these studies do not always 
provide a compelling explanation for dominance of the 
intrinsic mechanism 

Fully consistent metallic linear response theories of the 
AHE - Important progress has been achieved in AHE 
theory. The semiclassical theory, appropriately modi- 
fied to account for interband coherence effects, has been 
shown to be consistent with fully microscopic theories 
based on Kubo and Keldysh formalisms. All three the- 
ories have been shown to be equivalent in the e^^r 3> 1 
limit, with each having their advantages and disadvan- 
tages (Sec. IV). Much of the debate and confusion in 
early AHE literature originated from discrepancies and 
farraginous results from earlier inconsistent application 
of these linear response theories. 

A semiclassical treatment based on the Boltzmann 
transport equation, but taking into account the Berry 
curvature and inter-band coherence effects, has been for- 



mulated (Sec. IV.Al. The physical picture for each pro- 



cess of AHE is now understood reasonably well in the 
case of elastic impurity scattering. 

More rigorous treatments taking into account the 
multi-band nature of the Green's functions in terms of 
Kubo and Keldysh formalism have been developed fully 
(section IV.Cl. These have been applied to a particu- 
lar model, i.e., the ferromagnetic Rashba model, with a 
static impurity potential which produces elastic scatter- 
ing. The ferromagnetic Rashba model has an avoided 
crossing which has been identified as a key player in the 
AHE of any material. These calculations have shown a 
region of disorder strength over which the anomalous Hall 
conductivity stays more or less constant as a function 
of axx, corresponding to the intrinsic-dominated regime. 
The emergence of this regime has been linked to the topo- 
logical nature of the intrinsic contribution, analogous to 
the topologically protexted quantized Hall efi'ect. 
Emergence of three empirical AHE regimes — Based on 
the large collections of experimental results and indi- 
cations from some theoretical calculations, it is now 



becoming clear that there are at least three different 
regimes for the behavior of AHE as a function of axx'- 
(i) {<^xx > 10^ (ricm)^^) A high conductivity regime in 
which ffx^ ~ cFxx, skew scattering dominates cr^^ , and 
the anomalous Hall angle an/'^xx is constant. In this 
regime however the normal Hall conductivity from the 
Lorentz force, proportional to a^^H, is large even for 
the small magnetic field H used to align ferromagnetic 
domains and separating cr^jf^ and cr^^ is therefore chal- 
lenging, (ii) (10"* (f]cm)-i < axx < 10^ (f^cm)"!) An in- 
trinsic or scattering-independent regime in which is 
roughly independent of (Txx ■ In this intermediate metallic 
region, where the comparison between the experiments 
and band structure calculations have been discussed, the 
intrinsic mechanism is assumed to be dominant as men- 
tioned above. The dominance of the intrinsic mechanism 
over side-jump is hinted at in some model calculations, 
but there is no firm understanding of the limits of this 
simplifying assumption, (iii) {(Txx < 10* (ficm)^-'^) A 
bad-metal regime in which ax^ decreases with decreas- 
ing axx at a rate faster than linear. In this strong dis- 
order region, a scaling an o: a^x with 1.6 < n < 1.7 
has been reported experimentally for a variety of materi- 
als discussed in Sec|TT] This scaling is primarily observed 
in insulating materials exhibiting variable range hopping 
transport and where axx is tuned by varying T. The ori- 
gin of this scaling is not yet understood and is a major 
challenge for AHE theory in the future. For metallic ul- 
trathin thin films exhibiting this approximate scaling, it 
is natural that a^ is suppressed by the strong disorder 
(excluding weak localization corrections). Simple con- 
siderations from the Kubo formula where the energy de- 
nominator includes a {H/tY is that an oc when this 
broadening is larger than the energy splitting between 
bands due to the SOI. Since in this large broadening 
regime axx is usually no longer linear in r, an upper limit 
oi fi — 2 for the scaling relation an oc a^x is expected. 
The numerical Keldysh studies of the ferromagnetic 2D 
Rashba model indicates that this power is /? ~ 1.6, close 
to what is observed in the limited dirty-metallic range 
considered in the experiments. It is a surprising feature 
that this scaling seems to hold for both the metallic and 
insulating samples. 



b. Future challenges and perspectives 

In the classical Boltzmann transport theory, the resis- 
tivity or conductivity at the lowest temperature is simply 
related to the strength of the disorder. However, quan- 
tum interference of the scattered waves gives rise to a 
quantum correction to the conductivity and eventually 
leads to the Anderson localization depending on the di- 
mensionality. At finite temperature, inelastic scattering 
by electron-electon and/or electron-phonon interactions 
give additional contributions to the resistivity, while sup- 
pressing localization effects through a reduction of the 
phase coherence length. In addition one needs to con- 
sider quantum correction due to the electron-electon in- 
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teraction in the presence of the disorder. These issues, 
revealed in the 80's, must be considered to scrutinize the 
microscopic mechanism of a^x or p^x before studying the 
AHE. This means that it seems unlikely that only axx(T) 
characterizes the AHE at each temperature. Wc can de- 
fine the Boltzmann transport only when the residual 
resistivity is well defined at low temperature before the 
weak localization effect sets in. Therefore, we need to 
understand first the microscopic origin of the resistivity. 
Separating the resistivity into elastic and inelastic con- 
tributions via Mathhiessen's rule is the first step in this 
direction. 

To advance understanding in this important issue, one 
needs to develop the theoretical understanding for the ef- 
fect of inelastic scattering on AHE at finite temperature. 
This issue has been partly treated in the hopping the- 



ory of the AHE described in section II. B where phonon 
assisted hopping was assumed. However, the effects of 
inelastic scattering on the intrinsic and extrinsic mech- 
anisms are not clear at the moment. Especially, spin 
fluctuation at finite T remains the most essential and dif- 
ficult problem in the theory of magnetism, and usually 
the mean field approximation breaks down there. The 
approximate treatment in terms of the temperature de- 
pendent exchange splitting, e.g., for SrRuOa (Sec. II. B I, 
needs to be reexamined by more elaborated method such 
as the dynamical mean field theory, taking into account 
the quantum/thermal fluctuation of the ferromagnetic 
moments. These type of studies of the AHE may shed 
some light on the nature of the spin fluctuation in fer- 
romagnets. Also the interplay between localization and 
the AHE should be pursued further in the intermediate 
and strong disorder regimes. These are all vital issues 
for quantum transport phenomena in solids in general, 
as well as for AHE specifically. 

Admitting that more work needs to be done, in Fig. 
[44| we propose a speculative and schematic crossover di- 
agram in the plane of diagonal conductivity axx of the 
Boltzmann transport theory (corresponding to the dis- 
order strength) and the temperature T. Note that a 
real system should move along the y-axis as tempera- 
ture is changed, although the observed axx changes with 
T. This phase diagram reflects the empirical fact that 
inelastic scattering kills off the extrinsic skew scattering 
contribution more effectively, leaving the intrinsic and 
side-jump contributions as dominant at finite tempera- 
ture. We want to help stress that the aim of this figure 
is to promote further studies of the AHE and to iden- 
tify the location of each region/system of interest. Of 
course, the generality of this diagram is not guaranteed 
and it is possible that the crossover boundaries and even 
the topology of the phase diagram might depend on the 
strength of the spin-orbit interaction and other details of 
the system. 

There still remain many other issues to be studied in 
the future. First- principles band structure calculations 
for AHE are still limited to a few number of materi- 
als, and should be extended to many other ferromag- 
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FIG. 44 A speculative and schematic phase diagram for the 
anomalous Hall effect in the plane of the diagonal conductivity 
a^x and the temperature T. 



nets. Especially, the heavy fermion systems are an im- 
portant class of materials to be studied in detail. Con- 
cerning this point, a more economical numerical tech- 



nique is now available (Marzari and Vanderbilt 19971. 



This method employs the maximally localized Wannier 
functions, which can give the best tight-binding model 
parameters in an energy window of several tens of eV's 
near the Fermi level. The algorithm for the calculation of 
axy in the Wannier interpolation scheme has also been de- 



veloped (Wang et at. 20061. Application of these newly 



developed methods to a large class of materials should 
be a high priority in the future. 

AHE in the dynamical regime is a related interest- 
ing problem. The magneto-optical effects such as the 
Kerr/Faraday rotation have been the standard experi- 
mental methods to detect the ferromagnetism. These 
techniques usually focus, however, on the high energy re- 
gion such as the visible light. In this case, an atomic or 
local picture is usually sufficient to interpret the data, 
and the spectra are not directly connected to the d.c. 
AHE. Recent studies have revealed that the small energy 
scale comparable to the spin-orbit interaction is relevant 
to AHE which is typically ~ 10 meV for 3d transition 
meV in DMSs CSinova et al. 



metal and 



100 



20031. 



'xy\ 



This means that the dynamical response, i.e 
in the THz and infrared region will provide important 
information on the AHE. 

A major challenge for experiments is to find exam- 
ples of a quantized anomalous Hall effect. There are 
two candidates at present: (i) a ferromagnetic insula- 



tor with a band gap (Liu et al. 20081, and (ii) a disorder 
induced Anderson insulator with a quantized Hall con- 
ductance (Onoda and Nagaosa 20031. Although theo- 



retically expected, it is an important issue to establish 
experimentally that the quantized Hall effect can be re- 
alized even without an external magentic field. Such a 
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finding would be the ultimate achievement in identify- 
ing an intrinsic AHE. The dissipationless nature of the 
anomalous Hall current will manifest itself in this quan- 
tized AHE; engineering systems using quantum wells or 
field effect transistors is a promising direction to realize 
this novel effect. 

There are many promising directions for extensions of 
ideas developed through studies of the AHE. For exam- 
ple, one can consider several kinds of "current" instead 
of the charge current. An example is the thermal/heat 
current, which can be also induced in a similar fashion 
by the anomalous velocity. The thermoelectric effect has 
been discussed briefly in Sec. H.C where combining all 



the measured thermoelectric transport coefficients helped 
settle the issue of the scaling relation a 



AH 

xy 



metallic DMSs ( |Pu et a?.[|2008| ). Recent studies in Fe al- 



loys doped with Si and Co discussed i n Sec. |ILA| followed 
a similar strategy ( Shiomi et al. 2009 ) . From the temper- 



ature dependence of the Lorentz number, they identified 
the crossover between the intrinsic and extrinsic mech- 
anisms. Further studies of thermal transport will shed 
some light on the essence of AHE from a different side. 

Spin current is also a quantity of recent great inter- 
est. A direct generalization of AHE to the spin current 
is the spin Hall effect (|Kato et al. 2004 Murakami et al. 



20031 Sinova et al. 2004 Wunderlich et al. 20051, which 



can be regarded as the two copies of AHE for up and 
down spins with the opposite sign of axy. In this effect, 
a spin current is produced perpendicular to the charge 
current. An interesting recent development in spin Hall 
effect is that the quantum spin Hall effect and topological 
insulators have been theoretically predicted and experi- 
mentally confirmed. We did not include this exciting 
new and still developing topic in this review article. In- 
terested readers are referred to the original papers and 



references therein ( 


Bernevig et al. 


2006 


Kane and Mele 


2005 


Koenig et al. 


2007 
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